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Обозначения и сокращения 


В – (англ.) Бше, синий спектральный канал (длина волны около 
450-515 нм). 

РЪР - (англ.) 412 а1 $12па] ргосеѕѕог, цифровой сигнальный 
процессор. 

ЕТМ+ - (англ.) епһапсеа ћетайс таррегріиѕ, расширенный 
тематический картограф плюс (камера КА «Г.ап4$а(-7»). 

С – (англ.) ргееп, зеленый спектральный канал (длина волны 
около 525...605нм). 

СРСРО - (англ.) репегаі-ригро ѕертарһіс ѕргосеѕѕіпе ипііѕ, СРО 
общего назначения технология использования графического 
процессора видеокарты для выполнения общих вычислений 
(вычислений, не предназначенных для решения задач вывода 
компьютерной графики). 

1СА - (англ.) шдерепдепе сотропепї апа|у$1$, анализ 
независимых компонент. 

[5КО — (англ.) пап Ѕрасе Кеѕеагсһћ Ограпіѕайоп, Космическое 
агентство Индии. 

ЈАХА - (англ.) Јарап Аегоѕрасе Ехр|огайоп Арепсу, Японское 
аэрокосмическое агентство. 

КАК! - (англ.) Когеа Аегоѕрасе Кеѕеагеһ Іпѕ0иќе, Корейский 
авиационно-космический научно-исследовательский институт. 

МІР ~ (англ.) ти@ауег регсерігоп, многослойный перцептрон. 

МАЅА – (англ.) МаНопа| Аегопаийс ѕапа Ѕрасе Айтіпіѕігабоп, 
Национальное управление по аэронавтике иисследованию 
космического пространства (США). 

МІК - (англ.) пеаг-іпігагей, «ближний» инфракрасный спект- 
ральный канал (длина волны около750...1400 нм.). 

МЕСММ – (англ.) попПпеаг се1шаг пеига! пеѓуогк, нелинейная 
ячеистая нейронная сеть. 

МОАА - (англ.) Майопа! Осеапісапа Аїтоѕрһегіс Айтіпіѕіга- 
поп, Национальное управление океанических и атмосферных 
исследований (США). 

РСА - (англ.) ргіпсіра! сотропепї апа1уѕ1$, метод главных 
компонентов. 

К - (англ.) гед, красный спектральный канал (длина волны 
около 630...690 нм). 
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ВВМ - (англ.) геѕігісіеа Во|тапп тасһіпе, ограниченная сто- 
хастическая машина Больцмана. 

КМІР - (англ.) геситепі ти|]ауег регсерігоп, многослойный 
перцептрон с обратными связями. 

МО - (англ.) зедиепиа! пита Іорітіғайоп, последовательная 
минимальная оптимизация. 

ЗУМ – (англ.) заррой уесїог тасһіпе, машина опорных 
векторов. 

МІК – (англ.) ѕһогі-уауеІіепе їһпігагей, коротковолновый инф- 
ракрасный спектральный канал (длина волны 1400...3000 нм.). 

050$ – (англ.) Опцед Ѕїаѓеѕ СеоІоріса! Зигуеу, Геологическая 
служба США. 

\О 5-84 – (англ.) могіа реойейс ѕуѕіет 1984, всемирная геоде- 
зическая система, мировой референц-эллипсоид. 

ҰТА - (англ.) уіппег ѓакеѕа1, принцип «победитель получает 
все». 

ОТМ - (англ.) ипіуегѕа! їгапѕуегѕе Мегсаюг, универсальная 
поперечная проекция Меркатора. 

ГИС – геоинформационная система. 

ГИСЧ - генератор псевдослучайных чисел. 

ДЗЗ – дистанционное зондирование Земли. 

ИНС - искусственная нейронная сеть. 

КА - космический аппарат. 

ЛА – летательный аппарат. 

НЖМД - накопитель на жестких магнитных дисках. 

ОС – операционная система. 

ПО – программное обеспечение. 

СКО – среднеквадратичное отклонение. 

ТВ – телевизионное (изображение). 

ТПВ – тепловизионное (изображение). 

ФНЧ - фильтр низких частот. 

ЦКМ - цифровая карта местности. 

ЭВМ - электронная вычислительная машина. 


Введение 


Совмещение данных дистанционного зондирования Земли 
(данных ДЗЗ) и цифровых карт местности (ЦКМ) в настоящее 
время является актуальной задачей, что подтверждается востребо- 
ванностью совмещения при решении следующих задач обработки 
данных ДЗЗ: 

– задачи автоматического и полуавтоматического уточнения 
топографических карт по актуальным спутниковым снимкам; 

— задачи оцифровки картографического материала, имеюще- 
гося в бумажном виде, поскольку, как правило, привязка данного 
картографического материала к той или иной системе координат 
требует дополнительного уточнения; 

– задачи первичной привязки спутниковых снимков по данным 
о положении и ориентации космического аппарата (КА), сделав- 
шего снимок. 

Разработка методов и алгоритмов автоматического совмеще- 
ния данных ДЗЗ и ЦКМ позволит повысить оперативность и уро- 
вень автоматизации решения перечисленных задач. 

Исследованию научных вопросов, связанных с обработкой 
изображений (в частности, с совмещением разнородных изображе- 
ний), посвящены работы известных отечественных и зарубежных 
ученых. Большой вклад в развитие научных исследований в этой 
области внесли: Алпатов Б.А., Визильтер Ю.В., Джанджгава Г.И., 
Евтушенко Ю.Г., Еремеев В.В., Желтов С.Ю., Злобин В.К., Серге- 
ев В.В., Сойфер В.А. и др. Значительное внимание этой проблеме 
уделяют и зарубежные ученые: Башков Е., Блейхут Р., Понс Ж., 
Прэтт У., Форсайт Д., Фукунага К. 

Большой вклад в развитие научных исследований в области 
исследования научных вопросов, связанных с теорией искусствен- 
ного интеллекта, теорией машинного обучения и теорией оптими- 
зации, внесли: Вапник В.Н., Демидова Л.А., Каширин И.Ю.., 
Кохонен Т., Корячко В.П., Ле Кун Я., Локтюхин В.Н., Осовский С., 
Платт Д., Розенблатт Ф., Ручкин В.Н., Скворцов С.В., Фукушима К., 
Хайкин С., Харалик Р.М., Червоненкис А.Я. и др. 

Анализ публикаций в области теории и практики совмещения 
данных ДЗЗ и ЦКМ позволяет говорить о существовании ряда 
проблем в данной области, к которым относятся: 
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— отсутствие полностью автоматического подхода к решению 
задачи совмещения, основанного исключительно на анализе дан- 
ных ДЗЗ и ЦКМ без использования дополнительной информации о 
положении КА, о существовании реперных точек, выделенных 
экспертом на спутниковом снимке и ЦКМ и других данных. 

Таким образом, существующие методы являются либо полу- 
автоматическими с существенным и, в ряде случаев, ключевым 
участием эксперта, либо требуют дополнительную информацию, 
часто отсутствующую в наборах исходных данных к совмещению; 

– существенные временные затраты на выполнение совмеще- 
ния, обусловленные полуавтоматизмом существующих методов; 

— отсутствие практических способов автоматического контроля 
за качеством совмещения, поскольку означенный контроль выпол- 
няется модифицированным алгоритмом совмещения, а значит 
перенимает от него все его недостатки, связанные с существен- 
ными временными затратами и низким уровнем автоматизации. 

Сложность задачи совмещения существенно возрастает в си- 
туациях, требующих оперативного совмещения, поскольку в дан- 
ном случае появляется четкое ограничение на время выполнения 
совмещения, выполнение которого не гарантируют существующие 
полуавтоматические методы и алгоритмы. 

В книге описываются системы искусственного интеллекта и 
разрабатываются алгоритмы моделирования этих систем и чис- 
ленные методы оптимизации, решающие некоторые задачи обуче- 
ния данных систем, реализующие в совокупности процесс совме- 
щения данных ДЗЗ и ЦКМ и свободные от перечисленных недос- 
татков. 


Глава 1 


РАЗРАБОТКА И АНАЛИЗ 
МАТЕМАТИЧЕСКОЙ МОДЕЛИ 
ПРОЦЕССА СОВМЕЩЕНИЯ 
РАЗНОРОДНЫХ ИЗОБРАЖЕНИЙ 


1.1. Математическая модель процесса 
совмещения цифровых карт местности 
с данными дистанционного зондирования Земли 


Данные ДЗЗ (спутниковые снимки), получаемые от сенсоров, 
установленных на КА, можно представить в виде изображения 
(спутникового снимка): 


В=6,}; (11) 


где: 6, – пиксель спутникового снимка В, находящийся в 1-Й 
строке пикселей, в / -м столбце; 6, — спектральная яркость 
пикселя 6, в К -м спектральном канале; Г — количество строк 


пикселей в спутниковом снимке В; Л — количество столбцов 
пикселей в спутниковом снимке В; К – количество спектральных 
каналов в спутниковом снимке ДВ. 

Спектральные яркости 6, должны быть масштабированы в 


диапазон [-1;1] для удобства дальнейшей работы с ними. Для 
масштабирования спектральных яркостей из первоначального 
диапазона уровней квантования может быть использовано линейное 
преобразование, которое для О и уровней имеет вид: 


^ 


би 
ние 
0/2 


) 
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Ьл Е [0;0); Би є7, 


где р — первоначальное (квантованное) значение спектральной 
яркости. При О = 256 формула (1.2) приобретает вид (1.3): 


Бн 
Е 15 
"10705 ин 
Пусть соота – функция преобразования координат (х, у) 


пикселя Б., в координаты (х,, у, ) системы координат спутнико- 
вого снимка В. Функция сооға определяется как: 
(х;, уу) = сооға((ї, 3)); (1.4) 
хуу; ЕК. 


ЦКМ М является множеством точечных, линейных и полиго- 
нальных объектов, объединенных в векторные слои в соответствии 
с атрибутивной информацией. Одной из форм описания ЦКМ М 
является ее представление в виде множества точек (х, у) таких, 


ЧТО: 
Р В р (1.5) 
х, У, Ха» Хь Уши» Уак Є, 
— соответственно нижняя и верхняя границы диапа- 


где. х х 


тіп? “тіп 
зона координат ЦКМ по оси ОХ рамки карты; ук, у — со0т- 
ветственно нижняя и верхняя границы диапазона координат ЦКМ 
по оси ОТ рамки карты. 

Каждой точке (х, у)є М ставится в соответствие некоторая 
атрибутивная информация (наличие или отсутствие объекта в 
точке, тип объекта и т.п.). 


Функция сотр , описывающая модель совмещения спутнико- 
вого снимка В с ЦКМ М |, задается как: 


сотр' = сотр( В, М ‚соога ); (1.6) 
р =  сотр'В); 
сотр' = 4 соот" = сотр (сооға); 


{07 7)) = сотр", Ј)); 
а 
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где: В' — преобразованное изображение размером /' строк на .Л' 
столбцов; соога' – функция преобразования координат пикселей 


изображения В' в координаты системы координат спутникового 
снимка В. 


Множество пикселей сотр ((1, /)) < В' имеет, в общем 


Е а 
случае, мощность |сотр (1, /))|—1 (то есть отображение 
> 


сотр': (1, 1) >> {(1,7')} не является взаимно однозначным). Таким 
образом, имеет место быть зависимость: 
Вр = 6'.. <> сотр'(Ь Ј). (1.7) 

Качество совмещения сотр оценивается с помощью количест- 
венной меры ошибки совмещения Ё. К мере Е предъявляются 
следующие требования: 

• Е > 0 при увеличении качества совмещения; 

• Е > при уменьшении качества совмещения; 

° Е должна быть определена для УВ, М ‚сотр; 

6 Е>0; УВ, М,сотр. 

Одной из мер Ё , удовлетворяющей перечисленным свойствам, 
является среднеквадратичное отклонение (СКО) Е: 


(х' р» Ув) = соога (2, Г)); 
Р={р}; р=((1, Ј)(х, у)); 


' —_ 2 ' _ У 
сотр '((1, у))| 


Хн) 
Р | 
где Р – множество пар реперных точек р, выделенных на спут- 
никовом снимке В и ЦКМ М вручную. Ручное составление мно- 
жества Р позволяет достичь высокой (практически идеальной) 


точности соответствия точек, выделенных на спутниковом снимке 
В , точкам, выделенным на ЦКМ М. 


Размерность Ё, — размерность х и у, что обеспечивается, в 


ғтѕе 


А(сотр', р) = 


Вы (сотр Н Р) сы 


> 


том числе, отсутствием нормировки Е (требование Е — оо при 
уменьшении качества совмещения). Таким образом, физический 
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смысл Е, заключается в оценке вариации в конкретных вели- 
чинах между идеальным смещением (математическое ожидание 
оцениваемой величины) и реальным результирующим смещением 
(оцениваемая величина) — ошибкой функции сотр. 

Функция совмещения сотр предполагает выполнение следую- 
щих действий: 

1) поиск Ипа, на спутниковом снимке В и карте М пар 


совпадающих точек р, : 


Р, м, = ИИ „(В.М ,сооға); (1.9) 
Рем = Рем) Рам (ам Леа 5 Ое Усы); (1.10) 
2) расчет для каждой пары р„„ точек вероятности совпадения 
уайайу( Ри ), 
3) составление множества Р, пар совпадающих точек: 
Р.м Є Рим < УАПАЙУ( Ри ) > Ти, (1.11) 
Тм Е 10,1), 


где Г, есть минимально допустимая для включения в множество 
Р. „ вероятность (порог вероятности) того, что в рассматриваемой 


е 


паре р,„ точки совпадают друг с другом; 
4) расчет еуа/ по множеству Р, преобразования сотр": 


е 


сотр' = еуаКР„). (1.12) 
Этап поиска пар совпадающих точек р,„ является самым 


сложным в ходе совмещения спутникового снимка В и ЦКМ М. 
Под «сложностью» в данном случае понимается как временная 
сложность, так и сложность по памяти. 

Действительно, как показано в [1], один из самых эффектив- 
ных по времени алгоритмов поиска совпадающих точек (контурное 
корреляционное совмещение) имеет сложность О (10°№). В работе 
[1] алгоритм контурного корреляционного совмещения применяет- 
ся для поиска совпадающих точек на ракурсе трехмерной модели 
местности и данных от камер, установленных на борту летательно- 
го аппарата (самолета или вертолета), но принципиального упро- 
щения данный факт не дает — нет принципиальной разницы в поис- 
ке совпадающих точек на фактически двухмерной сцене (спутни- 
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ковом снимке) и трехмерной сцене в случае, если нет необходи- 
мости в перегенерации трехмерного ракурса. 

Оптимизация алгоритма контурного корреляционного совме- 
щения приводит к существенным расходам памяти, необходимым 
для организации разнообразных кэшей, в том числе, кэшей потоков 
в случае, если используется параллельная версия алгоритма, опи- 
санного в [1]. 

Алгоритм контурного корреляционного совмещения дает низ- 
кую точность совмещения и нуждается в дополнительной настрой- 
ке при изменении погодных условий, характера местности и ряда 
других факторов. Таким образом, актуальной является задача 
разработки алгоритмов совмещения, имеющих низкую временную 
сложность, дающих высокое качество совмещения и, по возмож- 
ности, низкие затраты по памяти. 

Функция уаПАЙу оценки вероятности совпадения пары р, 
имеет сложную природу, зависящую в общем случае от вида функ- 
ции Йпа,„ и от значений параметров алгоритма, реализующего 
данную фукнцию. Наиболее простой реализацией уаПаЙу является 


оценка вероятности уаЛАЙу, по гистограмме отклонений д, и 5,: 
(Хеди ? Ум) = сооға (1а, э адн )); 
0, = Хеди — Хейди | д, = У гаи а Р ; 


(5, - М5, )/5, 2 +05, - М5,)/5 


х тах 
2 


где д ид — максимальные значения д, и д, соответст- 


х тах у тах 


РЕ); 
уапайу, (р,„) = 1- —_,(1.13) 


венно. Таким образом, уа/:2іѓу, оценивает расстояние от р,,„ до 


моды гистограммы в случае, если гистограмма является унимо- 
дальной. 


В качестве уайАЙу можно использовать уаЙаЙу, . 

Как правило, процесс совмещения носит итерационный харак- 
тер. Поиск оптимального преобразования сотр' происходит до тех 
пор, пока Ё(сотр',Р) > Т, где Т. є К – максимально допустимая 
ошибка совмещения сотр', при которой процесс совмещения бу- 
дет остановлен. На каждой итерации предполагается постепенное 
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уменьшение Т, что приводит к росту |Р, |, и в свою очередь, 
приводит к тому, что для расчета преобразования сотр’ 


используется больше точек из Р 


гам! ейи 


. Таким образом, следствием 


постепенного уменьшения Т, является постепенное уменьшение 
вклада в сотр' пар точек р„„, оценка вероятности совпадения ко- 
торых уаПайЙу( р.) выполнена не точно. 


Функция еуа[ реализует алгоритм выбора конкретного вида 
функции сотр’ и расчета ее параметров. Выбор происходит, в 


основном, в зависимости от | Р, |, а также от ряда других факторов, 
таких как, например, распределение уайайу на Р, . В качестве 
функции сотр' используются следующие преобразования В -> В' 


и сооға — соога': 

• аффинное преобразование, линейный сдвиг и масштабирова- 
ние (применяется чаще всего для задачи совмещения данных ДЗЗ и 
ЦКМ); 

• преобразование Гельмерта [2]; 

• полиномиальные преобразования первого, второго, третьего и 
более высоких порядков; 

• трансформация «резинового листа» ((англ.) гиббегѕ ћееї) [3]. 


1.2. Анализ требований 
к источникам исходных данных 


1.2.1. Анализ требований 
к данным дистанционного зондирования Земли, 
получаемым от сенсоров космических аппаратов 


В качестве источников данных ДЗЗ можно рассматривать КА, 
соответствующие следующим критериям: 

• достаточная информативность; 

• доступность спутниковых снимков российскому потребите- 
лю (так, например, для некоторых КА действуют ограничения на 
точность данных, распространяемых среди потребителей-нерези- 
дентов страны, владеющих данными КА); 

° запуск КА в промежутке с 2012-го года по 2014-й год и не 
выработка КА расчетного срока службы. 
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Следующие факторы также являются существенными для 
выбора источника данных ДЗЗ: 

• СТОИМОСТЬ; 

• частота обновления; 

• возможность съемки на заказ интересующей области земной 
поверхности; 

• способ распространения данных. 

Под достаточной информативностью спутниковых снимков 
понимается соответствие снимков следующим критериям: 

• разрешение должно быть достаточным для определения 
отдельных малых населенных пунктов (деревень и хуторов). Как 
правило, таким разрешением является 30...60 метров на пиксель; 

• количество и тип спектральных каналов должны быть доста- 
точными для выделения на снимке лесонасаждений, полей, пусты- 
рей, дорог, населенных пунктов и водных объектов. Как правило, 
для выделения на снимке перечисленных объектов достаточно 
наличие синего (В-(англ.) Ыое, длина волны около 450...515 нм), 
зеленого (С-(англ.) отееп, длина волны около 525...605 нм), 
красного (К-(англ.) гей, длина волны около 630...690 нм) и 
нескольких ближних инфракрасных каналов (МК—англ.) пеаг- 
Іпітагеа, длин волны около 750...1400 нм). 

Таким образом, спутниковый снимок должен состоять из теле- 
визионного (ТВ; спектральные каналы В, С, В) и тепловизионного 
(ТПВ; спектральные каналы К и МК) изображений подстилающей 
поверхности с разрешением не более 30...60 метров на пиксель. 

Дистрибьютор данных ДЗЗ должен выполнять следующую пред- 
варительную обработку распространяемых спутниковых снимков: 

• радиометрическая и геометрическая коррекция спутниковых 
снимков; 

• привязка спутниковых снимков к системе координат \О5-84 
и проецирование спутниковых снимков в проекцию ОТМ. Пос- 
кольку данные система координат и проекция являются общеупот- 
ребительными, их поддержка имеется в большинстве современных 
ГИС (АгеС15, Егааѕ Ітаріпе, Оџапіит 01$, СКАЅЅ и другие), боль- 
шинство дистрибьюторов данных ДЗЗ распространяют спутни- 
ковые снимки. 

В настоящее время российские дистрибьюторы данных ДЗЗ 
(компания «Совзонд», ИТЦ «СканЭкс», НЦОМЗ и другие) рас- 
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пространяют продукцию ряда КА, из которых следующие могут 
рассматриваться как перспективные в плане получения данных ДЗЗ 
для решения поставленной задачи: 

• МогіаМіеу-1; 

• МУогіаміем-2; 

• Котрѕаѓ-2; 

• СеоЕуе-1; 

• АГО$ (камера АУМІК-2); 

• (Сагіоѕаѓ-2; 

• Кар:Еуе (группировка КА); 

• Ресурс-П; 

• Гапзае 7; 

• Гапза+ 8. 

Информация о КА «МогіаУіеу-1» и «МогіаУіем-2» сведена 
в табл. 1.1. Информация о КА «Котрѕаі-2» и «СеоЕуе-1» 
сведена в табл. 1.2. Информация о КА «АГО$» и «Ресурс-П» 
сведена в табл. 1.3. Информация о КА «Саоѕаі-2» и 
«КаріаЕуе» сведена в табл.1.4. Информация о КА «Г апӣѕаѓі7» и 
«Гап4за{8» приведена в табл. 1.5. 


Таблица 1.1. «КА М№огіаміем-1» и «“Уогаміеу-2» 


Уона\!ем -1 Уона\У!ем -2 


о 8 октября 2009 года 


Расчетный 7 лет 7 лет 
срок 

пребывания на 

орбите 


Режим съемки, Панхрома- Панхроматический: 500...900 нм 
спектральные тический: Мультиспектральный: Соаѕќа] 
диапазоны 500...900 нм («прибрежный»), В, С, ҮеПом 
(желтый), В, КеаЕаре (красный 
граничный), МІК (3МІКІ и МІК2) 


Пространст- Панхроматический режим: 0,5 м 
венное Мультиспектральный 
разрешение режим: 1,8 м 
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Окончание табл. 1.1. 


Компания 
Владелец РіркаіС1іоБе Компания Р1еиаіСіоБе (США) 
(США) 


Дистрибьютор Компания Компания «Совзонд» ИТЦ 


«Совзонд» ИТЦ «СканЭкс» НЦОМЗ 
Предваритель- 


«СканЭкс» 
НЦОМЗ 

ная обработка 

данных 


Радиомет- 
рическое 


Радиометрическая, 
сенсорная и 
геометрическая 
коррекция 
Система \О5-84/ОТМ и \!О5-84/ОТМ и другие 
координат, другие системы системы координат и проекции 


проекция координат и 
СеоТІЕЕ, МТЕ 


проекции 
СеотТІЕЕ, МТЕ 
1-4 дней (в зависимости от 
широты области съемки) 


Радиометрическая, сенсорная и 
геометрическая коррекция 


Формат 
предостав- 

ления данных 
Периодич- 
НОСТЬ 


1,7-5,9 дней (в 
зависимости от 
выполнения широты области 
съемки съемки) 


Срок предос- | день 
тавления дан- 
ных из архива 


Цена, Архив: 25 км Архив: 25 км“ (около 600 $) 
минимальный (около 600 $) Съемка на заказ: 25 км? 
объем заказа Съемка на заказ: (600-900 $) 

25 км? (600-900 $) 
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Таблица 1.2. «КА Котрѕаѓќ-2» и «СеоЕуе-1» 


28 июля 2006 года 6 сентября 2008 года 


Расчетный 5 лет 7 лет 
срок 

пребывания на 

орбите 


ГРРЕРРИНЕА Панхроматический: Панхроматический: 
500...900 нм 450...900 нм 
спектральные у 5 
Мультиспектральный: Мультиспектральный: 


диапазоны 


В, Сб, К, МК 
Панхроматический 


В, С, В, МВ 
Панхроматический режим: 


Пространст- 


венное режим: Ім 0,41 м Мультиспектральный 
разрешение Мультиспектральный режим: 1,65 м 
режим: 4м 


Радиометричес 10 бит на пиксель 11 бит на пиксель 
кое разрешение 
Агентство КАКІ 


Дистрибьютор | Компания «Ѕройтаре» Компания «Соорв1е» 
(Франция) Компания Компания «Ракурс» 
«Ракурс» Компания Компания «Совзонд» ИТЦ 
«Совзонд» ИТЦ «СканЭкс» 
«СканЭке» НЦОМЗ 


Предваритель- Радиометрическая, Радиометрическая, 
ная обработка сенсорная и сенсорная и геометрическая 
данных геометрическая коррекция 


коррекция (ІА, 2А, 26) 
\О5$-84/ОТМ и другие 
системы координат и 
проекции 


Сео ТТЕЕ 


Система 
координат, 
проекция 


УС95-84/0ТМ и другие 
системы координат и 
проекции 


Сео ТТЕЕ 


Формат 
предостав- 
ления данных 
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Цена, Архив: 50 км” (375- Архив: 49 км” (около 680$) 
минимальный 800$). Съемка на заказ: Съемка на заказ: 100 км? 
объем заказа 225 км? (3000-5700 $) (4000-5700$) 


Окончание табл. 1.2. 


Периодичность 
выполнения 
съемки 


1-3 дня (в зависимости от 
широты области съемки) 


Срок 
предоставле- 
ния данных из 


Таблица 1.3. «КА А08» и «Ресурс-ДК» 


Дата запуска 24 января 2006 года 25 июня 2013 года 


Расчетный 5 лет 5 лет 
срок пребыва- 
ния на орбите 


Режим съемки, Панхроматический: Панхроматический: 

спектральные 520...770 нм 580...800 нм Мультиспект- 

диапазоны Мультиспектральный: ральный: В, С, К, МК 
В, С, К, МВ (720...800 нм, 800...900 нм) 


Гиперспектральный: 96-255 
каналов (40...1100 нм) 


Пространст- Панхроматический Панхроматический режим: 
венное режим: 2,5 м 0,7 м Мультиспектральный 
разрешение Мультиспектральный режим: 3-4 м Гипер- 
режим:10 м спектральный режим: 25 м 

Радиометри- 8 бит на пиксель 10-14 бит на пиксель 
ческое 
разрешение 

Агентство ЈАХА Государственный научно- 


(Япония) производственный ракетно- 
космический центр «ЦСКБ- 


Прогресс» 
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Окончание табл. 1.3. 


Дистрибьютор Компания «Совзонд» Компания «Совзонд» 
ИТЦ «СканЭкс» НЦОМЗ 


Предваритель- Радиометрическая и В зависимости от 
ная обработка геометрическая дистрибьютора 
данных коррекция (1В2+КРС) 


Система \О5$-84/ОТМ и другие \05-84/ОТМ и другие 
координат, системы координат и системы координат и 


предоставле- 

ния данных 

Периодич- 46 суток 3 дня 
НОСТЬ 

съемки 

Срок предо- 3 дня В зависимости от 
ставления дан- дистрибьютора 
ных из архива 


Цена, Архив:1225 км? (около 
минимальный 1700 $). Съемка на 
объем заказа заказ: невозможна 


Таблица 1.4. «КА Сагіоѕаѓ-»2 и «КарійЕуе» 


КА Сагќоѕаќ -2 КарійЕуе 
(группировка КА) 


10 января 2007 года 29 августа 2008 года 


"кб ы 


пребывания на 
орбите 


Режим съемки, Панхроматический Мультиспектральный: 
спектральные : 500...850 нм В, С, К, Кеа Едзе (красный 
диапазоны раничный), МК 


Пространственное 1 м 5м 
разрешение 

Радиометрическое | 10 бит на пиксель 12-16 бит на пиксель 
разрешение 
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Окончание табл. 1.4. 


Владелец Агентство І5$КО Компания 
(Индия) «КаріаЕуеАС»(Германия) 


Дистрибьютор Компания «Ракурс» Компания «Совзонд» 
Компания 
«Совзонд» ИТЦ 
«СканЭкс» 
Радиометрическая и 
геометрическая 
коррекция. Приве- 
дение к картографи- 
ческой проекции 
Система ҮС8-84/0ТМ и \О5-84/ОТМ и другие 
координат, другие системы системы координат и 
проекция координат и проекции 

проекции 


СеоТІЕЕ Сео ЕЕ, МІТЕ 


Архив: 2500 км“ (цена дого- 
ворная). Съемка на заказ: 
5000 км? (цена договорная) 


Предварительная 
обработка данных 


Радиометрическая, 
сенсорная и геометрическая 
коррекция. Приведение к 
картографической проекции 


Формат 
предоставления 
данных 
Периодичность 
выполнения 
съемки 
Срок 
предоставления 

данных из архива 
Цена, 
минимальный 
объем заказа 


(590 $). Съемка на 
заказ: 25 км? (780 $) 


Таблица 1.5. «КА Гапф$а* 7» и «Гапф$а 8» 


15 апреля 1999 года 11 февраля 2013 года 


т В . 


срок 
пребывания 
на орбите 
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Окончание табл. [.5. 


Режим съемки, 


Панхроматический: Панхроматический: 


спектральные 520...900 нм 500...680нм Мультиспект- 
диапазоны Мультиспектральный: ральный: Соаѕќа! 
В, С, В, МІК, 5МІК («прибрежный и аэрозоли»), 
(760...900 нм В, С, В, МІК, 5%1К(1360- 


и 2080...2350 нм) 1390 нм, 1560...1660 нм, 
2100...2300 нм) 
Панхроматический 

режим:15 м 
Мультиспектральный 
режим: 30 м 
8(12 изначально) бит на 


пиксель 


МАЗА, 0505 МАЅА, 0805 
Дистрибьютор 0865 08565 


Предваритель- Радиометрическая, Радиометрическая, 
ная обработка сенсорная и геометри- | сенсорная и геометрическая 
данных ческая коррекция, при- коррекция, приведение к 
ведение к картографи- | картографической проекции 
ческой 
\О5-84/ОТМ и другие 
системы координат и 
проекции проекции 


Сео ТЕЕ СеоТТЕЕ 


Архив: бесплатно Архив: бесплатно 
Съемка на заказ: Съемка на заказ: 
невозможна невозможна 


Панхроматический 
режим: 15 м 
Мультиспектральный 
режим: 30 м 
8 бит на пиксель 


Пространст- 
венное 
разрешение 


Радиомет- 
рическое 


Система 
координат, 
проекция 


МС5-84/0ТМ и другие 
системы координат и 


Формат 
предоставле- 
ния данных 
Периодич- 
НОСТЬ 
выполнения 
съемки 
Срок 
предоставлени 
я данных из 


Цена, 
минимальный 
объем заказа 
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Из рассмотренных выше КА данные следующих КА доступны 
бесплатно: 

• данные, полученные от КА «СеоЕуе-1», доступны с помощью 
Интернет-портала «СоорІе Марѕ» [4]; 

• данные, полученные от КА «Іапаѕаї 7» и «Гапаѕаї 8», доступ- 
ны с помощью Интернет-портала «Еапйћ Ѕ$сіепсе аа Іпќегѓасе (Е$р1) 
а ће С1ІобаІ Гапа Соуег Еасіїїу» [5]. 

Таким образом, для решения исследовательских задач, не требую- 
щих регулярной съемки интересующей области земной поверхнос- 
ти и съемок на заказ, наилучшим выбором представляется исполь- 
зование данных ДЗЗ, полученных от КА «Іапдѕаї 8», по следую- 
щим причинам: 

• достаточное количество спектральных каналов приемлемого 
разрешения; 

• бесплатность; 

• наличие радиометрической и геометрической коррекции; 

• распространение данных с привязкой к системе координат 
\О5-84 в проекции ОТМ; 

• запуск КА в 2013-м году и невыработка им расчетного срока 
службы. 


1.3. Методы и алгоритмы совмещения 
цифровых карт местности с данными 
дистанционного зондирования Земли 


1.3.1. Краткая характеристика 
существующих методов и алгоритмов совмещения 


В настоящее время существует несколько подходов к решению 
задачи совмещения данных ДЗЗ и ЦКМ, среди которых можно 
выделить следующие. 


Методы совмещения, основанные наручном выделении 
на спутниковом снимке и ЦКМ опорных точек 


Одним из наиболее распространенных подходов к совмещению 
данных ДЗЗ и ЦКМ является подход, основанный на выделении на 
спутниковом снимке и ЦКМ множества пар опорных точек, по 
которым в дальнейшем выполняется расчет результирующего пре- 
образования совмещения. 


22 Разработка и анализ математической модели 
процесса совмещения разнородных изображений 


Парой опорных точек называются такие точки, для которых 
соответствующие позиции были идентифицированы как на спут- 
никовом снимке, так и на ЦКМ. Поскольку опорные точки, прос- 
тавленные вручную, должны точно отображаться пространствен- 
ным преобразованием, они являются ограничением, которому дан- 
ное пространственное преобразование должно удовлетворять. 

В качестве пространственного преобразования, определяющего 
отображение спутникового снимка на ЦКМ (и соответственно 
опорных точек, выделенных на спутниковом снимке, в соответст- 
вующие точки, выделенные на ЦКМ), могут выступать преобразо- 
вание из группы движения, преобразование подобия, проективное 
преобразование, аффинное преобразование, трансформация «рези- 
нового листа». 

Одним из методов, реализующих данный подход, является ме- 
тод, описанный в [6] и применяемый на первом этапе техноло- 
гической схемы пространственной привязки спутниковых снимков 
по ЦКМ, описанной в [6], для выполнения первоначального (гру- 
бого) совмещения. Данный метод состоит в выделении пар опор- 
ных точек вручную оператором ГИС, после чего система выпол- 
няет расчет результирующего преобразования. 

Недостатком данного метода является его неавтономность из- 
за необходимости присутствия оператора, выполняющего проста- 
нов опорных точек на спутниковом снимке и ЦКМ. В отличие от 
данного метода совмещения, нейросетевые методы являются пол- 
ностью автоматическими. 


Методы совмещения, основанные 
на полуавтоматическом и автоматическом поиске 
на спутниковом снимке и ЦКМ особых точек 


Развитием методов, основанных на выделении множества 
опорных точек, являются методы, решающие задачу выделения 
опорных точек полуавтоматически или полностью автоматически. 

Полуавтоматические методы совмещения выполняют выделе- 
ние множества опорных точек, верификацию которого производит 
оператор вручную. Данные методы реализованы в большинстве 
современных ГИС (например, в Егдаз Ітаріпе) и являются доста- 
точно популярными, однако не позволяют выполнять автоматичес- 
кое совмещение данных ДЗЗ и ЦКМ. 
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Методом совмещения, предполагающим полуавтоматическое 
выделение на спутниковом снимке множества опорных точек, 
является метод, описанный в [7] и основанный на построении 
множества опорных точек по наземным ориентирам (в [7] в 
качестве такового ориентира выступает вагон-лаборатория). Дан- 
ный метод предполагает поиск на спутниковом снимке объектов, 
точные координаты которых на ЦКМ известны заранее, форми- 
рование из точек, соответствующих данным объектам, множества 
опорных точек и расчет по сформированному множеству результи- 
рующего преобразования. 

Для достижения полной автоматизации процесса совмещения 
данных ДЗЗ и ЦКМ применяется несколько подходов -— в 
частности, часто реализуется подход, основанный на применении 
детекторов особых точек для выделения оных на спутниковых 
снимках и данных ДЗЗ. 

В качестве детекторов особых точек применяются детекторы 
ЪГЕТ, ОКЕ, Бедета, Фестнера, Харриса [8, 9]. Перечисленные де- 
текторы реализуют эвристические алгоритмы выделения особых 
точек, основанные на поиске изгибов контурных линий, резких 
перепадов яркости и прочих особенностей изображений. 

Недостатком методов совмещения, основанных на применении 
перечисленных детекторов, является их низкая точность при реше- 
нии задач совмещения разнородных изображений, коими являются 
спутниковые снимки и данные ДЗЗ. Это объясняется, прежде всего, 
тем, что перечисленные детекторы ориентированы на использова- 
ние при решении задач совмещения кадров видео-последователь- 
ностей или изображений, одно из которых является изображением 
части сцены, присутствующей на другом. 

В отличие от рассмотренных методов совмещения, нейросете- 
вые методы реализуют полностью автоматический подход к совме- 
щению, позволяющий выполнять совмещение разнородных изобра- 
жений, одно из которых (изображение карты) является искусствен- 
ным и может быть выполнено в искусственных цветах. 
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Методы совмещения, основанные 
на автоматическом поиске одноименных объектов 
на снимке и векторной карте 


Логическим развитием методов совмещения данных ДЗЗ и 
ЦКМ по множеству опорных точек является метод, основанный на 
автоматической или полуавтоматической идентификации объектов 
карты [10]. 

Данный метод предполагает поиск совпадений между образами 
объектов, выделенными на спутниковом снимке, и объектами, при- 
сутствующими на ЦКМ. Образ объекта и объект, совпадение меж- 
ду которыми идентифицировано с высокой вероятностью, позволя- 
ют включить в множество опорных точек соответствующие точки, 
лежащие в границах образа и объекта, например, точки изгибов их 
контурных линий. 

В качестве искомых объектов могут выступать объекты 
различных классов. Так, в [10] рекомендуется искать совпадения 
между образами и объектами гидрографии (речная сеть, береговые 
линии морей и озер). 

Данный подход к совмещению реализует нейросетевой метод 
совмещения, основанный на классификации образов. В отличие от 
метода, рассмотренного в [10], данный метод использует для поис- 
ка совпадений между образами объектов и объектами, присутст- 
вующими на ЦКМ, искусственные нейронные сети и другие 
концепции теории искусственного интеллекта. 


Корреляционно-экстремальные методы совмещения 


Корреляционно-экстремальные методы совмещения [11] основа- 
ны на решении задачи оптимизации целевой функции-меры сходства 
между спутниковым снимком и ЦКМ. 

Для предотвращения ложных срабатываний методов, возни- 
кающих из-за изменения значений пикселей в поисковой области, 
мера сходства нормируется. Однако, как показано в [12], данные 
методы обладают слабой устойчивостью к геометрическим искаже- 
ниям (за исключением сдвига). Для корректной работы методов 
часто необходимы дополнительные данные, так, например, для 
повышения точности совмещения данных ДЗЗ и ЦКМ предвари- 
тельно необходимо выполнить высокоточное совмещение с ЦКМ 
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базового спутникового снимка, по совмещению с которыми будет 
выполняться дальнейшее совмещение других, перекрывающих его, 
спутниковых снимков. 


Контурно-корреляционные методы совмещения 


Частным случаем корреляционно-экстремального подхода к 
совмещению является подход, реализующий контурное корреля- 
ционное совмещение. 

Данные методы совмещения основаны на выделении контуров 
образов объектов на спутниковом снимке В , контуров объектов на 
ЦКМ М и дальнейшем сопоставлении (корреляции) контуров друг 
с другом. 

Метод совмещения, описанный в [1], предполагает выделение 
на изображении, совмещаемом с ЦКМ, контуров с применением 
оператора Собеля, после чего коррелятор, основанный на поиске 
характерных точек изгибов контуров, выполняет сопоставление 
отдельных участков контуров образов объектов, выделенных на 
изображении, контурами объектов на ЦКМ. 

Как было отмечено в [1], данный метод совмещения может 
быть использован при решении достаточно узкого круга задач, та- 
ких, как совмещение данных, поступающих от камер, установлен- 
ных на летательном аппарате (ЛА), при заходе ЛА на посадку. 
Применение данного метода при решении задач совмещения 
данных ДЗЗ и ЦКМ ограничено и требует существенной ручной 
настройки параметров. 

В [13] предложен алгоритм совмещения, основанный на поиске 
на изображении В окон, содержащих достаточно точное (в смысле 
некоторой меры расстояния) контурное представление окна, выде- 
ленного на ЦКМ. Для означенного поиска в [13] предлагается ис- 
пользовать генетический алгоритм, применение которого повыша- 
ет автоматизм процесса совмещения, однако существенно увеличи- 
вает временные затраты на решения данной задачи. 

В [6] описывается метод совмещения, основанный на поэтап- 
ном сканировании спутникового снимка В с некоторым окном, 
выделении контуров с помощью простейшего дифференциального 
оператора и сопоставлении очередного окне окном, выделенным на 
ЦКМ М. В качестве меры сходства в [6] предлагается исполь- 
зовать меру, учитывающую протяженность контура, его толщину и 
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направление. Данный подход хорошо поддается распараллели- 
ванию, поэтому его реальные временные характеристики позво- 
ляют использовать его для оперативного решения задач совмеще- 
ния, однако точность данного подхода непосредственно зависит от 
качества выделения контуров и алгоритма фильтрации помех, 
вызванных облаками, дымкой и другими факторами. 

В [14] рассмотрен подход к совмещению спутниковых снимков 
В, полученных от геостационарных КА, основанный на коорди- 
натной привязке снимков В по ЦКМ М и контурным точкам дис- 
ка Земли. Подход основан на поиске соответствий между контур- 
ными линиями водных объектов на ЦКМ и контурами образов 
объектов на спутниковом снимке в некоторой зоне (области) 
поиска, задаваемой автоматически в ходе выполнения совмещения, 
конфигурация которой зависит от направления контура, его прог- 
нозируемой толщины, его конфигурации. Для минимизации влия- 
ния облачности на качество совмещения в [14] предлагается ис- 
пользовать подходы, основанные на идентификации облаков на 
снимке, оценке степени перекрытия ими контуров водных объектов 
и принятии решения об использовании или исключении из рас- 
смотрения отдельных участков контуров, перекрытых облаками. 

Данный подход обеспечивает высокую точность координатной 
привязки и, следовательно, высокую точность совмещения данных 
ДЗЗ и ЦКМ, однако он не может быть использован при обработке дан- 
ных ДЗЗ, полученных от КА, расположенных на негеостационарных 
орбитах, что существенно ограничивает применимость подхода. 

Перечисленные методы совмещения данных ДЗЗ и ЦКМ име- 
ют следующие отличия от нейросетевых подходов: 

• от нейросетевого подхода, основанного на поиске образов 
объектов, присутствующих на ЦКМ, на спутниковом снимке, рас- 
смотренные методы отличаются принципом работы – нейросете- 
вой подход выполняет поиск соответствий между образами объек- 
тов на спутниковом снимке В и объектами, присутствующими на 
ЦКМ М, не рассматривая структуру образов и объектов, тогда как 
анализ структуры границ образов и объектов является основной 
особенностью перечисленных методов; 

• от нейросетевого подхода, основанного на поиске корреля- 
ции между контурами образов, выделенных на спутниковом сним- 
ке, и границами объектов, присутствующих на ЦКМ, рассмотрен- 
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ные методы отличаются тем, что для выполнения контурной корре- 
ляции нейросетевой подход использует нейросетевой коррелятор, 
который позволяет достичь различных эффектов, не достижимых 
рассмотренными методами, в частности, нейросетевой коррелятор 
способен автоматически восстанавливать мелкие разрывы конту- 
ров, что положительно сказывается на качестве совмещения. 

В целом, применение методов совмещения, основанных на ис- 
пользовании систем искусственного интеллекта, позволяет осу- 
ществлять более точное совмещение из-за возможности учета раз- 
личного рода помех, таких как зашумленность спутникового сним- 
ка, влияние облаков и дымок и т.д., в автоматическом режиме, 
тогда как рассмотренные методы контурно-корреляционного сов- 
мещения требуют для достижения похожих результатов исполь- 
зовать дополнительные эвристические алгоритмы. Кроме того, кон- 
турно-корреляционные методы совмещения имеют больше наст- 
роек, чем предлагаемые нейросетевые подходы, что понижает уро- 
вень автоматизма процесса решения задачи совмещения. 

Однако необходимо отметить, что временные характеристики 
(без учета времени поиска параметров запуска алгоритмов) рас- 
смотренных контурно-корреляционных подходов к совмещению 
чуть лучше, чем временные характеристики нейросетевых подходов. 


Метод стабилизации изображений, основанный 
на преобразовании Фурье и вейвлет-преобразованиях 


Методы стабилизации изображений выполняют поиск парамет- 
ров и вида совмещения путем сравнения результатов различных 
интегральных преобразований между совмещаемыми спутниковым 
снимком и ЦКМ [15]. 

Основным недостатком данных методов является их низкие кор- 
ректирующие способности — данные методы позволяют, как прави- 
ло, компенсировать только сдвиг спутникового снимка относительно 
ЦКМ, при чем компенсация сдвига выполняется в достаточно узком 
диапазоне. Кроме того, данные методы, как правило, содержат в себе 
алгоритмы повышения разрешения, что обусловливает высокие 
требования по отношению сигнал/шум на спутниковом снимке. 

В отличие от методов стабилизации, нейросетевые методы 
позволяют компенсировать не только сдвиг спутникового снимка 
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относительно ЦКМ, но и другие, более сложные, преобразования – 
такие как, например, поворот и масштабирование. 


1.3.2. Методы и алгоритмы совмещения, основанные 
на использовании искусственных нейронных сетей 
и нейроподобных концепций 


Рассмотренные подходы к решению задачи совмещения 
спутниковых снимков и ЦКМ имеют ряд недостатков. В настоящей 
работе рассматриваются системы искусственного интеллекта, 
алгоритмы моделирования систем и численные методы обучения 
систем, реализующие методы совмещения данных ДЗЗ и ЦКМ, 
свободные от перечисленных недостатков. 

Оптимальным методом совмещения разнородных изображе- 
ний с точки зрения минимизации средней ошибки на большой 
выборке тестовых данных является метод, реализуемый зритель- 
ными системами высокоорганизованных животных и человека [16]. 
На основании данного утверждения оказывается, что наилучшими 
подходами к автоматическому совмещению данных ДЗЗ и ЦКМ с 
точки зрения качества совмещения являются методы, максимально 
правдоподобно имитирующие деятельность зрительных систем вы- 
сокоорганизованных животных и человека. К таким методам отно- 
сятся подходы, основанные на моделировании систем искусст- 
венного интеллекта, включающих в себя искусственные нейронные 
сети (ИНС) различного типа и другие нейроподобные концепции 
теории искусственного интеллекта. 

ИНС и нейроподобные конструкции обладают следующими 
достоинствами: 

• способность к обобщению информации. 

ИНС, выполняющая классификацию образов и обученная на 
репрезентативной выборке тестовых образов, вследствие своей 
нелинейной природы оказывается способна [17-19] к обобщению 
информации, что выражается в выделении общего набора характер- 
ных скрытых признаков для каждого класса образов. Данная спо- 
собность ИНС и нейроподобных конструкций оказывается полезна 
при выполнении совмещения путем поиска соответствий между 
образами на спутниковом снимке и объектами на ЦКМ; 

способность к запоминанию и восстановлению информации. 
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ИНС и нейроподобные конструкции, реализующие ассоциатив- 
ную память, способны [17-19] к восстановлению образов, ранее за- 
помненных ими в процессе обучения. 

• возможность распараллеливания вычислений, выполняемых в 
процессе обучения и использования ИНС и нейроподобных конст- 
рукций. Модульная структура ИНС и нейроподобных конструк- 
ций позволяет реализовывать [20] параллельные алгоритмы их обу- 
чения и использования для практически любых типов параллель- 
ных систем, к которым относятся параллельные системы, взаимо- 
действующие через разделяемую память (многопроцессорные парал- 
лельные системы, системы, основанные на использовании техно- 
логии общих вычислений на видеокартах ((англ.) гепега|-ригрозе 
етарһісѕ ргосеѕѕіпе ипіїѕ, СРСРО, технология выполнения вычис- 
лений на графическом процессоре, не предназначенном для реше- 
ния задач вывода компьютерной графики)), параллельные системы, 
взаимодействующие путем передачи сообщений (вычислительные 
кластеры, в том числе, вычислительные кластеры, организованные 
по технологии облачных вычислений). 

Таким образом, интеллектуальные системы, построенные с ис- 
пользованием ИНС и нейроподобных структур, оказываются спо- 
собны к горизонтальному масштабированию, ограниченному толь- 
ко контекстом вычислений (решаемой задачей). 

Необходимо отметить возможность распараллеливания подхо- 
дов к выбору структур ИНС и нейроподобных структур. Одним из 
таковых подходов является подход, основанный на использовании 
генетического алгоритма, который сам по себе достаточно хорошо 
поддается распараллеливанию [19-22]. 

К таковому аппаратному обеспечению относятся специализи- 
рованные нейропроцессоры (подклассцифровых сигнальных процес- 
соров (англ.) Ча] ѕірпа! ргосеѕѕог, ЮЅР), к которым относится 
семейство процессоров МеигоМайлх [23]. Использование нейропро- 
цессоров позволяет реализовывать ИНС и нейроподобные конструк- 
ции на аппаратном уровне, что приводит к существенному увели- 
чению производительности системы искусственного интеллекта. 

Существуют следующие подходы к решению задачи совме- 
щения данных ДЗЗ с ЦКМ, основанные на использовании ИНС и 
нейроподобных концепций теории машинного обучения. 
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Подход, основанный на поиске образов объектов, 
присутствующих на ЦКМ, на спутниковом снимке 


Данный подход предполагает выполнение следующих дейст- 
ВИЙ: 

1) выделение на спутниковом снимке характерных образов 
(1.14); 

2) классификация характерных образов (1.15); 

3) поиск соответствий между образами, выделенными на спут- 
никовом снимке, и объектами, присутствующими на ЦКМ: 

Р = па (В, М,сооға) = 


сопјоғт(В, М ‚соога ,сіаѕѕійсайоп(ѕерағайоп(В)), 4е5с); 


Н = ѕерағайоп(В); (1.14) 

С = сІаѕѕійсайоп(Н ‚4езс); (1.15) 

Р а, = сопрюгт(В, М ,сооға,С); (1.16) 
Н={й}; РСВ; 
Сее сен; 


С теб с еи 
у = деѕс(й); у= {у,У,,...,у)}5 Уу ЄК; Өє\; (1.17) 
п єс; ЙЕН. 
Функция 5ерагайоп выполняет выделение образов объектов на 
спутниковом снимке В. Образы л формируют множество выде- 
ленных образов Н. 


Множество образов Н является аргументом функции класси- 
фикации с/аѕѕійсаііоп . Другим аргументом функции с/аѕ51йсаііоп 


является функция 4е5с (1.17), выполняющая расчет векторов приз- 
наков у образов й, выделенных на спутниковом снимке. Резуль- 
татом функции с/аѕѕ Йсайоп является множество классов С = {с}. 


Классы не пересекаются друг с другом, то есть каждый образ й 
должен быть отнесен только к одному классу. 

Множество классифицированных образов (множество классов) 
С является аргументом функции соијоғт . Функция сопюгт 


выполняет поиск соответствий между образами из классов С и 
объектам и на ЦКМ М с использованием информации о классах 
образов и их положении на спутниковом снимке. Результатом 
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функции сопрюгт является множество пар совпадающих точек 
Р, 


гам ейи ? 
пар точек Р, (1.11) и функции сотр" (1.12). 
Функция 5ераганоп может быть основана на использовании 


следующих методов: 

• сегментация на основе выделения границ [10]; 

• алгоритм порогового преобразования [9, 10]; 

° алгоритмы наращивания областей (алгоритм центроидного 
связывания, алгоритм слияния-расщепления) [10]; 

• алгоритм морфологического водораздела [9]; 

• метод релаксационной разметки [10]; 

• кластеризация изображения с использованием многомерной 
нейронной карты Кохонена [20, 24, 25]; 

• кластеризация изображения с использованием алгоритма 
Кк-теапѕ [26]; 

• некоторые другие методы. 

Сегментация на основе выделения границ, алгоритм порогово- 
го преобразования, алгоритмы наращивания областей, алгоритм 
морфологического водораздела и метод релаксационной разметки 
существенным образом зависят от освещения сцены. Например, 
некорректно выбранные параметры алгоритма выделения границ 
для сегментации на основе выделения границ, зависящие от осве- 
щенности сцены, ведут к выделению отдельных элементов образов 
объектов сцены, либо к объединению отдельных образов объектов 
сцены вследствие нечеткости их границы. Для учета условий осве- 
щенности сцены данные алгоритмы используют дополнительные 
эвристические подходы к расчету своих параметров, что приводит 
к усложнению процесса совмещения и увеличению вычислитель- 
ной (а следовательно, и временной) сложности выделения образов. 

Кластеризация изображений с использованием многомерной 
нейронной карты Кохонена и алгоритма А-теапѕ в меньшей сте- 
пени зависят от освещенности сцены. Алгоритм А-теап5, по срав- 
нению с многомерной нейронной картой Кохонена, предоставляет 
меньше возможностей по своей настройке (алгоритм А-теапѕ ис- 
пользует только Евклидово расстояние, тогда как многомерная 
нейронная карта Кохонена может работать с любой вещественной 
метрикой). 


по которым в дальнейшем выполняется расчет множества 
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Таким образом, наилучшим выбором для выделения образов на 
сцене оказывается многомерная нейронная карта Кохонена. 

В качестве функции с/а55саНоп могут быть использованы 
ИНС, нейроподобные конструкции и прочие концепции теории 
машинного обучения, обучаемые с учителем и решающие задачу 
построения регрессии: 

а = Г(й,4езс); (1.18) 
ає\, 
где 4 – номер класса, к которому относится образ й. 

Классификация выполняется по вектору признаков у = деѕс(й) 
образа й, что позволяет говорить о процессе классификации как о 
процессе разбиения на непересекающиеся подмножества (классы) 
векторного пространства Г = {у}, содержащего векторы признаков 
образов Й. 

Различные классификаторы, обучаемые с учителем, выполня- 
ют разбиение векторного пространства Г на классы различными 
способами. Параметры разбиения рассчитываются по обучающему 
множеству образов Н, — множеству заранее классифицирован- 
ных образов й. Расчет параметров называется обучением класси- 
фикатора. Качество выполненного обучения проверяется на 
тестовом множестве образов Н — множестве заранее классифи- 


1е51 
иН 


зентативными подмножествам и множества Н 


цированных образов й. Множества Н являются репре- 


теас 1651 


заранее 


Іеасћ 1еѕі 


классифицированных образов й, при чем выполняются условие: 
Ни» УН, = Ә . Множества Н, и Н, формируются путем 


Іеасћ 1е51 


случайной выборки из Н 


теасй 1е51 * 


Для классификации могут быть использованы следующие 
простые (не составные) классификаторы: 

• дерево принятия решений [9]. 

Дерево принятия решений, используемое для решения слож- 
ных с точки зрения количества классов и размерности пространст- 
ва признаков классификационных задач, состоит из достаточно 
большого количества узлов, требующих тонкой первоначальной 
настройки и дальнейшей подстройки в зависимости от особеннос- 
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тей задачи, что ведет к необходимости эвристической настройки 
классификатора практически для каждой сцены; 

• наивный Байесовский классификатор [8, 9]. 

Наивный Байесовский классификатор позволяет получить 
результат, оптимальный со статистической точки зрения. Данный 
классификатор позволяет выполнять линейное разделение классов, 
чего в ряде задач недостаточно из-за многомерности пространства 
признаков и сложности линейного разделения искомых классов; 

• скрытая Марковская модель [26]. 

Основными проблемами при использовании скрытых Марковс- 
ких моделей являются их колеблющаяся производительность, 
неприемлемая для решения поставленной задачи совмещения, и 
сильные структурные ограничения модели — используемые допу- 
щения обусловной независимости могут затруднить введение неко- 
торых относительно простых условий; 

• нейронная сеть прямого распространения без обратных связей 
[21, 22] (частный случай — многослойный перцептрон [17, 18, 26]). 

Существенным недостатком нейронной сети прямого рас- 
пространения и ее частного случая — многослойного перцептрона — 
является отсутствие алгоритма обучения, позволяющего достичь 
глобального минимума функции ошибки выполняемой классифика- 
ции. Действительно, алгоритмы обучения нейронной сети прямого 
распространения — алгоритм градиентного спуска, квазиньютоновс- 
кие методы обучения (например, алгоритм Левенберга-Марквард- 
та), алгоритм сопряженных градиентов и прочие алгоритмы обуче- 
ния позволяют достичь лишь локального минимума и требуют до- 
полнительных эвристик для повышения вероятности схождения в 
глобальный минимум функции ошибки [27]. Таким образом, ней- 
ронные сети прямого распространения плохо подходят для решения 
сложных задач классификации – задач, линейное разделение классов 
в которых является достаточно сложным либо невозможным вовсе. 

Влияние данного недостатка на точность классификации мини- 
мизируется с помощью специальных эвристических алгоритмов 
инициализации весов синапсов ИНС (алгоритм Нгуэна-Видроу, 
алгоритм Ле Куна), а также путем выбора структуры ИНС с помо- 
щью генетического алгоритма [21, 22]. 
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Достоинством нейронной сети прямого распространения явля- 
ется низкая временная сложность алгоритмов ее обучения и рабо- 
ты, а также возможность распараллеливания их выполнения; 

• сверточные ИНС-когнитрон, неокогнитрон, сверточная ней- 
ронная сеть Ле Куна, обученная с использованием методики «глу- 
бокого» обучения ((англ.) деер Іеагпіпе; (англ.) еер БейеЁ пебмогК) 
[28-32]. 

Обучение сверточных сетей является достаточно ресурсоемкой 
задачей и не гарантирует схождение функции ошибки к глобаль- 
ному минимуму [29]: 

. машина опорных векторов ((англ.) ѕиррогі уесіог тасћһіпе, 
ЅУМ) [17, 26]; 

• прочие классификаторы. 

Таким образом, наилучшим выбором для решения задачи 
классификации образов оказываются нейронные сети прямого 
распространения без обратных связей, машины опорных векторов 
и сверточная нейронная сеть Ле Куна. 

Нейронные сети прямого распространения, обладая низкой 
временной сложностью алгоритмов обучения и работы, могут быть 
эффективно использованы в составе классификатора, решающего 
задачу классификации образов в контексте поставленной задачи. 
Увеличение вероятности схождения сети к глобальному минимуму 
функции ошибки возможно с помощью эвристических алгоритмов 
инициализации весов синапсов нейронов и путем подбора струк- 
туры нейронной сети с использованием, например, генетического 
алгоритма. 

Машина опорных векторов позволяет выполнять оптимальное 
с точки зрения достижения глобального минимума функции ошиб- 
ки нелинейное разбиение классифицируемого векторного прост- 
ранства на два класса. 

Простые классификаторы могут быть объединены в составные 
классификаторы, содержащие функциональные блоки комбиниро- 
вания выходов отдельных классификаторов (например, несколько 
классификаторов, значения на выходах которых есть вероятности 
принадлежности образа к каждому из соответствующих классов, 
могут быть объединены в составной классификатор, выходом 
которого будет номер класса, вероятность принадлежности образа 
к которому максимальна). Как правило, алгоритм классификации 


Глава 1 35 


предполагает использование нескольких простых классификаторов, 
объединяемых в составной (сложный) классификатор [17]. Данное 
объединение может быть осуществлено следующими способами: 

объединение классификаторов в двоичное дерево [33, 34]. Объе- 
динение простых классификаторов в двоичное дерево позволяет 
минимизировать сложность задачи классификации, решаемой клас- 
сификатором, помещенным в узел дерева. Действительно, каждый 
узел классификатора должен решить простейшую из возможных 
задач — задачу бинарной классификации; 

• объединение классификаторов с использованием Воо$(-алго- 
ритмов [17]. В случае объединения классификаторов в составной с 
использованием Вооѕі-алгоритмов результаты работы классифи- 
каторов уравновешивают друг друга, предоставляя на выходе Во- 
0${-алгоритма наиболее вероятный результат классификации, что 
также существенным образом повышает качество классификации. 

Недостатком Вооѕі-алгоритмов является необходимость допол- 
нительных временных затрат на поиск оптимальной модели объе- 
динения классификаторов; 

• организация экспертной системы (например, использование 
системы голосования) [17]. 

Экспертная система является наиболее эффективным средством 
объединения классификаторов с точки зрения качества работы сос- 
тавного классификатора, однако в процессе ее использования тре- 
буются существенные временные затраты на ее создание и обучение, 
что в контексте задачи оперативного совмещения данных ДЗЗ и 
ЦКМ делает экспертную систему очевидно неприменимой; 

• другие способы комплексирования классификаторов. 

Наилучшим средством объединения простых классификаторов 
в составной с точки зрения качества классификации и времени 
обучения классификатора оказывается объединение простых клас- 
сификаторов в дерево — таким образом, каждый узловой классифи- 
катор будет решать простую задачу разделения векторного под- 
пространства с выхода предыдущего узла дерева на два класса. 

Расчет векторов признаков классифицируемых образов (функ- 
ция деѕс ) может основываться на следующих алгоритмах: 

расчет минимального (усредненного) максимального вектора 
спектральных яркостей пикселей образа [35]. 
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Недостатком данного подхода является невозможность описа- 
ния им пространственной структуры текстур образов. Данный под- 
ход, будучи примитивным, но в то же время вычислительно эф- 
фективным, может быть использован для классификации образов в 
простейших случаях – например, в случае от деления класса обра- 
зов, соответствующих водным объектам, до класса образов, соот- 
ветствующих лесу, то есть в случаях, когда основная информация 
для работы классификатора содержится в спектральных яркостях 
пикселей образов и, следовательно, структура текстур образов 
может быть проигнорирована; 

• расчет текстурных признаков пикселей: 

— описание текстур с помощью расчета статистических харак- 
теристик выходов блоков пространственных фильтров [26]. 

К недостаткам данного подхода к описанию текстур относят- 
ся интуитивность составления блока пространственных фильт- 
ров, сложность формализации подходов по автоматическому 
составлению блоков и, вообще говоря, отсутствие каких бы то 
ни было критериев оптимальности составления блоков, кроме 
тестирования блоков на реальной или модельной задаче. 

Отдельно необходимо упомянуть проблемы, возникающие в 
случае, если блок фильтров слишком мал или, наоборот, слишком 
велик. В первом случае возможно возникновение ситуации, при 
которой блок не даст достаточно информации для качественного 
решения поставленной задачи классификации образов, выделенных 
на изображении. Во второй ситуации, наоборот, возможен избыток 
информации, при котором происходит переобучение классифика- 
тора, выраженное в том, что количество векторов описаний, нахо- 
дящихся на границе разделения различных классов, становится 
пренебрежимо малым по сравнению с векторами описаний образов, 
типичных для каждого из классов, итак, граница разделения классов 
в существенной степени размывается, что негативно влияет на 
качество работы классификатора. Кроме того, в случае большого 
блока фильтров существенно увеличиваются временные затраты на 
получение описания текстур, что обусловливается необходимостью 
свертки исходного изображения В с каждым из фильтров; 

— использование энергетических характеристик Лавса [8, 36, 37]. 
Энергетические характеристики Лавса обладают следующи- 
ми достоинствами: 


Глава 1 


37 


возможность компактного описания большинства воз- 
можных характерных паттернов структуры текстуры 
таких, как периодическое повторение элементов текс- 
туры, множество не связных мелких элементов, отдель- 
ные относительно большие пятна и прочие; 
возможность компактного описания уровня яркости 
текстуры, масштабированного к одинаковому уровню 
освещения по всему спутниковому снимку; 
Недостатком данного подхода является необходимость 
свертки исходного изображения с набором пространст- 
венных фильтров, ведущая к существенным времен- 
ным затратам на расчет описания текстуры; 
— использование текстурных признаков Харалика [37, 38]. 
Достоинством текстурных признаков Харалика является их 
большая информативность и большая гибкость характера 
описания по сравнению с рассмотренными подходами. 
Недостатком текстурных признаков Харалика является их 
меньшая компактность, по сравнению с энергетическими 
характеристиками Лавса, что объясняется необходимостью 
расчета одних и тех же характеристик для нескольких 
матриц вхождений; 
— прочие способы описания образов, выделенных на изоб- 
ражении. 
Таким образом, в контексте поставленной задачи наилуч- 
шим выбором является использование следующих подхо- 
дов к описанию образов, выделенных на изображении: 
— использование энергетических характеристик Лавса. 
Данный подход к описанию образов используется на 
верхних уровнях дерева классификаторов, то есть на тех 
уровнях, на которых происходит разделение существенно 
отличных друг от друга классов, например, класса водных 
объектов от класса лесных насаждений; 
— использование текстурных признаков Харалика. 
Данный подход к описанию образов используется на ниж- 
них уровнях дерева классификаторов, поскольку он хорошо 
зарекомендовал себя для описания текстур в задачах 
разделения классов, если граница между классами является 
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существенно нелинейной (например в случае разделения 
класса хвойного леса и класса лиственного леса). 

Необходимо также отметить тот факт, что, как правило, © 
(размерность вектора признаков) является достаточно большой, что 
ведет к увеличению временной сложности процесса классифи- 
кации, поэтому целесообразно использование того или иного мето- 
да понижения размерности пространства признаков. К таковым 
методам относятся: 

• попарное одношаговое или многошаговое усреднение векто- 
ров признаков. Попарное одношаговое или многошаговое усредне- 
ние векторов признаков может быть использовано в контексте 
решения достаточно простых задач-задач, в которых возможно 
линейное разделение классов, поскольку усреднение не позволяет 
учитывать и регулировать потери информации, происходящие при 
понижении размерности пространства признаков; 

• применение преобразования Карунена-Лоева, матрица кото- 
рого может быть рассчитана с использованием нейронной сети, 
реализующей анализ главных компонентов ((англ.) ргіпсіра! сотро- 
пепї апа[уѕ1ѕ, РСА) [17, 18, 39, 40]. 

Для понижения размерности пространства признаков можно 
использовать преобразование Карунена-Лоева, поскольку оно поз- 
воляет оценивать и регулировать потери информации при пониже- 
нии размерности пространства признаков. 

Функция соијоғт описывает многоэтапный алгоритм, состоя- 


щий из поискового алгоритма и алгоритма вероятностной оценки 
соответствия образа й, выделенного на спутниковом снимке В, 
объектам ЦКМ М , расположенным вблизи точек { сооға((/, /)) }; 


УВ, ЕЙ [41]. 


Рассмотренный подход к совмещению обладает следующими 
достоинствами: 

° низкие временные затраты на выполнение совмещения, что 
связано с широкими возможностями по распараллеливанию отдель- 
ных шагов алгоритма вследствие их нейросетевой природы, а так- 
же с тем, что для совмещения достаточно, как правило, нахож- 
дения не более нескольких десятков пар образов и объектов; 
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• высокая точность совмещения, обусловленная гибкой струк- 
турой классификатора, реализующего статистически оптимальную 
классификацию образов, выделенных на спутниковом снимке. 


Подход, основанный на поиске корреляции между 
контурами образов, выделенных на спутниковом снимке, 
и границами объектов, присутствующих на ЦКМ 


Данный подход предполагает выполнение следующих 
действий: 

1) предварительная обработка спутникового снимка с целью 
устранения несущественных (слабых или прерывистых) контуров 
(1.19); 

2) выделение контуров на спутниковом снимке (1.20); 

3) описание контуров (1.21); 

4) поиск совпадений контуров образов на спутниковом снимке 
с границами объектов на ЦКМ (1.22). 


Ром = Лт (В, М,сооға) = 


= соггёанНои( М ,саеѕс(сопіоиғ( ргераге(В)))); 


В „„ = р'ерағе(В), (1.19) 

В =сотоиг(В„„); (1.20) 

С = саезс(В „); (1.21) 
Ри = сОТТ@аноп(М,С); (1.22) 


В ЕЕ ИЕ ЕО: ОЕ 


Исходный спутниковый снимок В является аргументом функ- 
ции предварительной обработки ргераге ‚ результатом которой 


является изображение В на котором, в отличие от спутникового 


ргр ? 
снимка В, отсутствуют несущественные (слабые или прерывис- 
тые) контуры. 


Изображение В, является аргументом функции сотоиг , ко- 


торая выполняет выделение контуров на изображении. Результатом 
функции сопюиг является контурная матрица В. , каждый пик- 


спі ? 


сель В которой содержит меру силы контура, проходящего 


сп и 


через пиксель В исходного изображения В„,. Под мерой силы 


р'р іј 
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контура в (1.20) понимается численная мера, принимающая зна- 
чения из диапазона [0;1], где 0 соответствует отсутствию контура, а 


| – максимально сильному контуру. В простейшем случае В 


ст 


является бинарной матрицей. 

Контурная матрица В, является аргументом функции опи- 
сания контуров саеѕс ‚ результатом которой является множество С 
описаний = контуров, выделенных на изображении В„,. 


Множество С вместе с ЦКМ М являются аргументами функ- 
ции соггеаНоп , выполняющей корреляционное совмещение ЦКМ 
и спутникового снимка. В ходе корреляционного совмещения 
контурам, описание которых присутствует в С, ищутся соот- 
ветствия на ЦКМ (Границы объектов), по результатам чего 
формируется множество пар совпадающих точек Р по кото- 


га ейи ? 
рым в дальнейшем выполняется расчет множества пар точек Р, 
(1.11) и функции сотр" (1.12). 

Алгоритм выделения контуров, описываемый функцией 
сопіоиғ, должен удовлетворять следующим требованиям: 

низкие временные затраты; 

• точность, достаточная для дальнейшего автоматизированного 
совмещения спутникового снимка с ЦКМ. 

Последнее требование предполагает выделение не всех, хоро- 
шо различимых для человеческого глаза, контуров, но выделение 
наиболее протяженных контуров и их дополнительное (в зави- 
симости от характера контуров) сглаживание. Так, например, 
сглаживание может потребоваться при выделении контура водного 
объекта, имеющего растительность, произрастающую на его берегу 
и отсутствующую на карте. Для обеспечения такового выделения 
необходимо выполнение предобработки спутникового снимка, для 
чего используется функция ргераге. 

Существуют следующие алгоритмы выделения контуров на 
изображении: 

• алгоритмы, основанные на дифференциальных операторах: 

— алгоритм, основанный на использовании перекрестного 
градиентного оператора Робертса [9]; 


Глава 1 41 


— алгоритм, основанный на использовании оператора 
Собеля (или частной версией оператора Собеля — оператора 
Щура) [9]; 

— алгоритм, основанный на использовании оператора 
Превитт [9]; 

— детектор краев Кэнни [8, 9]. 

Данные алгоритмы выделяют большое количество слабых 
контуров (в том числе, контуров, наличие которых на изоб- 
ражении обусловлено разного рода помехами), непригод- 
ных для решения задачи автоматизированного совмещения 
спутникового снимка с ЦКМ, что ведет к увеличению 
вычислительной и временной сложностей алгоритма совме- 
щения и снижению качества совмещения; 

• метод Хафа. 

Данный метод предполагает существенные вычислительные 
затраты, связанные с большим количеством операций с плавающей 
ТОЧКОЙ. 

Рассмотренные алгоритмы выделения контуров не удовлетво- 
ряют в полной мере перечисленным выше условиям. 

Методами выделения контуров, свободными от недостатков 
перечисленных алгоритмов, являются: нейросетевой подход к вы- 
делению контуров, основанный на использовании вейвлета Габора 
[42], и нейросетевой алгоритм выделения контуров, основанный на 
использовании нелинейных ячеистых искусственных нейронных 
сетей ((англ.) поп Ппеаг сеПшаг пеига! пебмогк, М.СММ) [43- 46]. 

Для описания саеѕс контуров, выделенных на изображении, 
могут использоваться следующие подходы: 

• формирование вектора пар координат пикселей спутникового 
снимка В , занятых контуром [9]. 

Недостатком подхода, описывающего контур путем формиро- 
вания вектора пар координат пикселей, занятых им, является его 
чрезмерная информационная перегруженность, что особенно за- 
метно на изгибах контуров, следствием которой является снижение 
временной эффективности работы алгоритма совмещения и высо- 
кая чувствительность к шумам различной природы; 

расчет цепочечного кода контура [8, 9]. К недостаткам цепо- 
чечного кода контура относится его существенная абстрагирован- 
ность от конкретной задачи — действительно, в контексте корреля- 
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ционного совмещения с использованием ассоциативной памяти ис- 
пользование абстрактного кода приводит к необходимости его 
адекватного масштабирования в диапазон вещественных значений, с 
которыми работает ассоциативная память, что не всегда возможно; 
составление множества подматриц матрицы В, фиксирован- 


ного размера (множества окон), содержащих контур [46, 47]. 

Данный способ описания контуров не учитывает минимальные 
изменения на его протяжении, но оказывается достаточно инфор- 
мативным, простым (с временной точки зрения) для расчета и 
хорошо подходящим для совмещения с использованием ассоциа- 
тивной памяти; 

• прочие способы описания контуров. 

Для описания контуров на изображении можно применять 
подход, основанный на составлении множества подматриц матри- 
цы В. фиксированного размера (множества окон), содержащих 


контур. 

Множество С описаний контуров, выделенных на спутнико- 
вом снимке В, является одним из аргументов функции соггеайоп. 
Другим аргументом функции соггёаНоп является ЦКМ М . Функ- 
ция соғғеіаїоп выполняет корреляционное совмещение спутнико- 
вого снимка В с ЦКМ М, которое заключается в поиске соот- 
ветствий между контурами, описания которых наличествуют в С, 
и границами объектов, присутствующих на ЦКМ М. 

Существует ряд подходов к корреляционному совмещению: 

• алгоритм предварительного совмещения изображений на ос- 
нове выделения преобладающих углов характерных линий, описан- 
ный в [1]. 

Методы корреляционного совмещения, один из которых был 
писан в [1], обладают рядом недостатков, наиболее существен- 
ными из которых являются необходимость дополнительной наст- 
ройки методов для каждого спутникового снимка вследствие су- 
щественного влияния погодных условий, условий освещенности, 
разнообразных прочих помех на результат совмещения и низкое 
качество результата совмещения, что позволяет использовать дан- 
ные алгоритмы лишь при решении узкоспециализированных задач, 
которыми являются, например, задачи определения характеристик 
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(класс, аттрибутивная информация) объекта на спутниковом сним- 
ке по ЦКМ; 

• алгоритм совмещения, основанный на использовании поэтап- 
ного сканирования спутникового снимка [48]. 

Поэтапное сканирование предполагает фактически простой 
перебор возможных множеств Р с расчетом качества совме- 


гам’ ейи 


щения для каждого. Простой перебор в данном случае предпола- 
гает существенные временные затраты на выполнение совмещения 
спутникового снимка и ЦКМ. 

Существуют несколько эвристических подходов к оптимиза- 
ции алгоритма поэтапного сканирования, однако ни один из них 
существенным образом не улучшает временные характеристики 
алгоритма; 

• алгоритм совмещения, основанный на использовании генети- 
ческого алгоритма [13]. 

Одним из подходов к улучшению временных характеристик 
алгоритма поэтапного сканирования является подход, основанный 
на использовании генетического алгоритма для совмещения дан- 
ных ДЗЗ и ЦКМ. Генетический алгоритм, обладая временными ха- 
рактеристиками, лучшими по сравнению с простым перебором, все 
равно не применим в задачах, требующих оперативного совмеще- 
ния. Кроме того, генетический алгоритм предъявляет существен- 
ные требования к объему оперативной памяти вычислительной 
системы и, таким образом, не применим во встраиваемых системах; 

• подходы, основанные на использовании адаптивной памяти: 

— подходы, основанные на использовании адаптивно-резо- 
нансной теории [29], подходы, основанные на использова- 
нии ИНС Хопфилда и ИНС Хемминга [17, 18]. 

Концепции адаптивно-резонансной теории, ИНС Хопфилда 
и ИНС Хемминга не могут рассматриваться в качестве серь- 
езного инструмента для построения высокоточных систем 
совмещения данных ДЗЗ и ЦКМ, поскольку ИНС работают с 
двоичными образами (и, следовательно, требуют специаль- 
ных преобразований между вещественным описанием кон- 
тура и его двоичным описанием) и способны запоминать 
лишь относительно не большое количество восстанавливае- 
мых образов. Оба недостатка обусловливают существенные 
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затруднения в процессе использования данных нейросете- 
вых концепций и провоцируют горизонтальный рост сис- 
темы искусственного интеллекта, реализующей совмещение; 
— подходы, основанные на использовании нейронной сети 
Коско (двунаправленная ассоциативная память) [29, 49]. 
Нейронная сеть Коско представляет интерес, в большей 
степени, как модель человеческого мозга; 

— подходы, основанные на использовании многослойного 
перцептрона с обратными связями ((англ.) гесштепі ти @- 
[ауег регсерігоп, КМГР) [18]. 

Основным применением многослойного перцептрона с 
обратными связями является удаление шумов из сигналов, 
хотя теоретически КМГР может быть использован для 
реализации адаптивной памяти; 

— подходы, основанные на использовании ограниченной 
стохастической машины Больцмана ((англ.) геѕігісіеа Вой2- 
тапп тасһіпе, ВВМ) [17, 46, 47]. 

Данная концепция обеспечивает высокую точность восста- 
новления, способна запоминать существенно больше обра- 
зов по сравнению с концепциями адаптивно-резонансной 
теории, ИНС Хопфилда и ИНС Хемминга и может быть 
эффективно реализована в системе, поддерживающей па- 
раллельные вычисления. 


Таким образом, единственной концепцией, которая может быть 
использована для решения задачи совмещения спутниковых снимков 
и ЦКМ, является ограниченная стохастическая машина Больцмана. 

Рассмотренный подход, так же как и подход, основанный на 
поиске объектов, присутствующих на ЦКМ, на спутниковом сним- 
ке, обладает следующими достоинствами: 

° низкие временные затраты на выполнение совмещения, что 
связано с хорошими возможностями по распараллеливанию от- 
дельных этапов совмещения вследствие их нейросетевой природы; 

• высокая точность совмещения, что обеспечивается нелиней- 
ностью процесса выполнения совмещения. 
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МЕТОДЫ СОВМЕЩЕНИЯ РАЗНОРОДНЫХ 
ИЗОБРАЖЕНИЙ ПУТЕМ ПОИСКА 
И СОПОСТАВЛЕНИЯ ОБРАЗОВ ОБЪЕКТОВ 


2.1. Алгоритм моделирования системы 
искусственного интеллекта 


На вход системы искусственного интеллекта, выполняющей 
совмещение данных ДЗЗ и ЦКМ путем поиска образов объектов, 
присутствующих на ЦКМ, на спутниковом снимке, поступают 
следующие данные: 

В – исходный спутниковый снимок; 

М - исходная ЦКМ; 

сооға — функция преобразования координат (1, /) пикселя 65,, 


спутникового снимка В в координаты (х,, у, ) точки на ЦКМ М 


в соответствии с системой координат и проекцией спутникового 
снимка В. 
На выход системы искусственного интеллекта поступает мно- 


жество пар совпадающих точек Р (1.10), по которым в даль- 


гам’ еди 
нейшем выполняется расчет множества пар точек Р., (1.11) и 
функции сотр’ (1.12), описывающей результирующее совмещение 


спутникового снимка В и ЦКМ М. 

Система искусственного интеллекта состоит из следующих мо- 
дулей: 

1) модуль выделения на спутниковом снимке характерных об- 
разов (2.1); 

2) модуль классификации характерных образов (2.2); 

3) модуль поиска соответствий между образами, выделенными 
на спутниковом снимке, и объектами, присутствующими на 
ЦКМ (2.3). 


2а ейи — Лпа а, (В, М,сооға) = 


сопјоғт(В, М, соога ,сіаѕѕійсаіоп(ѕерағайоп(В),аеѕс)); 


46 Методы совмещения разнородных изображений путем поиска 
и сопоставления образов объектов 


Н = 5ераганои( В); (2.1) 
С = сІаѕѕійсайоп(Н ,аеѕс); (2.2) 
Рл = сопютт(В,М ‚соога,С); (2.3) 
Н ={№}, йСВ; 
Са сен; 
ЕС сес, еее; 
у = 2еѕс(й); у= {у,У,,...,у,}; Му ЄК; Өє\, (2.4) 
й ЕС; ЙЕН. 


Модуль 5ерагаНоп выполняет выделение образов объектов на 


спутниковом снимке В. Образы л формируют множество выделен- 
ных образов Н. 

Множество образов Н подается на вход модуля классификации 
сІаѕѕійсаіоп . Другим входным параметром модуля классификации 


сІаѕѕійсаіоп является модуль деѕс (2.4), выполняющий расчет век- 


торов признаков у образов й, выделенных на спутниковом снимке. 
На выход модуля с/аѕѕійсаіоп поступает множество классов 


С = {с}. Классы не пересекаются друг с другом, то есть каждый 


образ й должен быть отнесен только к одному классу. 
Множество классифицированных образов (множество классов) 
С подается на вход модуля сопюгт. Модуль сопогт выполняет 


поиск соответствий между образами из классов С и объектами на 
ЦКМ М с использованием информации о классах образов и их 
положении на спутниковом снимке. На выход модуля соп/огт 


поступает множество пар совпадающих точек Г. 


га’ ейи ? 


по которым в 


дальнейшем выполняется расчет множества пар точек Р, (1.11) и 


функции сотр' (1.12). 

В $ 1.3.2. приведен детальный анализ возможных подходов к 
реализации каждого из модулей системы. Окончательно были 
выбраны следующие подходы: 

• модуль 5ераганоп — многомерная нейронная карта Кохонена; 


• МОДУЛЬ аеѕс — алгоритмы описания образов, основанные на 
использовании энергетических характеристик Лавса и текстурных 
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признаков Харалика, алгоритм понижения размерности прост- 
ранства признаков, основанный на использовании РСА-сети; 
• модуль с/аѕѕійсайоп — древовидный классификатор, ИНС 


прямого распространения без обратных связей и машина опорных 


векторов; 


Получение 
В, М 


Получение 
настроек 
системы 


Понижение размер- 
ности простран- 
ства признаков 


Классификация 
образов 


Поиск соответ- 
ствий между 
образами на Ви 
объектами на М 


Выделение на 
В характер- 
ных образов Н 


Возврат 


Р.а еди 


Рис. 2.1. Алгоритм моделирования 
системы искусственного интеллекта 


• модуль сопјоғт — многоэтапный алгоритм, состоящий из по- 
искового алгоритма и алгоритма вероятностной оценки соот- 
ветствия образа й, выделенного на спутниковом снимке В, объек- 
там ЦКМ М, расположенным вблизи точек {соога ((ї, /))}; УБ, ЕЙ. 
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Схема алгоритма моделирования системы искусственного 
интеллекта приведена на рис. 2.1. Алгоритм моделирования сос- 
тоит из следующих этапов: 

1) выделение на спутниковом снимке характерных образов 
(2.1) с использованием многомерной нейронной карты Кохонена; 

2) описание выделенных образов (2.4) с использованием энер- 
гетических характеристик Лавса и текстурных признаков Харалика; 

3) понижение размерности пространства признаков с использо- 
ванием РСА-сети; 

4) классификация характерных образов (2.2) с помощью древо- 
видного классификатора; 

5) поиск соответствий между образами, выделенными на спут- 
никовом снимке, и объектами, присутствующими на ЦКМ (2.3). 


2.2. Выделение образов на спутниковом снимке 
2.2.1. Структура, принципы обучения и функционирования 


Существует несколько подходов к выделению образов на спут- 
никовом снимке (функция ѕерағаѓіоп в (2.1)). Как было показано в 


$ 1.3.2, наиболее эффективным с точки зрения времени выполне- 
ния и требуемых вычислительных ресурсов и позволяющим полу- 
чить при этом результат качества, достаточного для дальнейшего 
совмещения, является подход, основанный на использовании мно- 
гомерной нейронной карты Кохонена [24]. 

Пусть Г = {у} – векторное пространство, определенное над 
[-1;] и содержащее векторы спектральных яркостей (| у |= К). Для 
спутникового снимка В справедливо, таким образом, Вс И. 

Алгоритм выделения образов состоит из следующих этапов: 

1) классификация векторов спектральных яркостей пикселей 


спутникового снимка В с помощью многомерной нейронной кар- 
ты Кохонена; 


2) кластеризация спутникового снимка В с получением мно- 
жества (= {4} кластеров. 


Каждый кластер 4 является замкнутым множеством пикселей 
из В, отнесенных к одному и тому же классу. 
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Алгоритм кластеризации состоит из построчного сканирования 
В, классификации каждого пикселя и постепенного наращивания 
формируемых кластеров; 

3) фильтрация выделенных кластеров с формированием на 
выходе множества Н выделенных образов. Для каждого кластера 
4 рассчитывается количественная оценка диаШу полезности 


кластера для дальнейшего совмещения, с помощью которой фор- 


мируется множество Н образов (2.5), выделенных на спутниковом 
снимке: 


Н = {4 | диаііу(а) > > Т ТЕ р (2.6) 
Н ону Є [0; Ц, 
где ГРЕТА пороговое значение, определяющее минимальную 


оценку полезности кластера 4 , при которой кластер будет включен 


в множество Н выделенных образов. Значение Г является 


диайу 
одним из настраиваемых параметров алгоритма совмещения. 
К диаШу предъявляются следующие требования: 


• оценка должна быть определена всюду на Г и непрерывна на 
всей области определения; 

областью значений оценки является диапазон [0; 1]; 

• значение оценки должно быть тем больше, чем выше полез- 
ность кластера 9 для дальнейшего совмещения. 


В качестве дама» могут быть использованы различные оцен- 


ки. Экспериментальные исследования показали, что наилучшей с 
точки зрения временных затрат и качества результата является 
оценка, основанная на вычислении площади кластера д: 
диау(а) = 221097, ол) 
04 
а= 175 

где © — коэффициент нормировки, равный количеству пикселей на 
спутниковом снимке В (т.е. равный площади спутникового снимка 
В ‚ выраженной в пикселях). 

Ключевым элементом алгоритма выделения образов является 
многомерная нейронная карта Кохонена, выполняющая классифи- 
кацию пикселей спутникового снимка В. 
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Многомерная нейронная карта Кохонена выполняет разбиение 
исходного векторного пространства Г на набор классов { У, }; 


Г. < И, каждый из которых состоит из векторов, ближе всего рас- 


положенных к центральному вектору у, данного класса в смысле 


используемой меры расстояния 4. 

Нейронная карта представляет собой многомерную сетку 
нейронов, каждый из которых рассчитывает меру расстояния меж- 
ду входным вектором и центральным вектором соответствующего 
класса. 

Пусть 5 = {5\,...,5,›...5}; 5, ЄМ – вектор размерностей карты. 
Каждому нейрону, лежащему в узловой точке карты с коорди- 
натами "={\....,\,,....\.} ‚ можно поставить в соответствие 


вектор у, , центральный для Г, -го выделяемого класса, где 
ү 
№ = + У (и, – 1) – номер класса. Важно отметить, что жє М. 


п=2 
Вектор 5 является одним из параметров алгоритма совмещения и 
выбирается оператором. 

Нейронная карта относится к классу ИНС, обучаемых без 
учителя. Обучение нейронной карты осуществляется с помощью 
алгоритма самоорганизации нейронной карты Кохонена. 

Алгоритм классификации пикселей спутникового снимка В с 
помощью многомерной нейронной карты Кохонена состоит из 
следующих этапов: 

1) инициализация нейронной карты — выбор параметров, ини- 


циализация векторов у, весов нейронов; 

2) самоорганизация карты на множестве ВсИ; 

3) прогон карты для каждого пикселя ре В, отнесение пикселя 
к одному из допустимых классов /,. 


Принцип работы многомерной нейронной карты Кохонена 
проиллюстрирован рис. 2.2. На рис. 2.2 трехмерная карта описы- 
вает двухмерное векторное пространство в соответствии с его ста- 
тистическими характеристиками. Несоответствие размерности (кар- 
та большей размерности описывает пространство меньшей раз- 
мерности) позволяет достичь нелинейности описания, что поло- 
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жительно влияет на качество выполняемой классификации И 
соответственно на качество сегментации пространства. 


Нелинейное 
описание картой 


двухмерного 
пространства 
Х Трехмерная карта Самоорганизация 
Кохонена нейронной карты 


Рис. 2.2. Сегментация двухмерного пространства 
с помощью трехмерной карты Кохонена 


В качестве меры расстояния 4 на векторном пространстве Г 
могут быть использованы различные метрики. Экспериментально 
оказалось, что использование Евклидова расстояния (2.8) в 
качестве таковой метрики является лучшим выбором с точки 
зрения временных характеристики с точки зрения оценки качества 
совмещения. 


т 
аб) = У ә). (2.8) 
к=] 

В случае использования Евклидова расстояния в качестве 
метрики 4 задача обучения нейронной карты сводится к построе- 
нию разбиения Воронова множества В – т.е. к разбиению, мини- 
мизирующему дисперсию в результирующих классах относительно 
соответствующих центральных векторов. 


2.2.2. Инициализация весов нейронов карты 
с помощью датчика псевдослучайных чисел. 
Выбор датчика псевдослучайных чисел 


Существуют следующие подходы к инициализации весов ней- 
ронов нейронной карты Кохонена: 

• инициализация весов нейронов нулевыми векторами; 

• инициализация весов нейронов случайными значениями из 
диапазона [1',, Г... |; 
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• инициализация весов нейронов случайными значениями из 
диапазона [М,-Т., М +Г. |. 


тіл ° тах 
По результатам экспериментальных исследований для инициа- 
лизации весов нейронов нейронной карты Кохонена в контексте 
решения задачи совмещения данных ДЗЗ и ЦКМ была выбрана 
инициализация весов случайными значениями из диапазона [7.,, 
Т |. Принимая во внимание ограничение у = {у, }; у, Е[-1;1] обус- 


тах 
лавливаемое характером предметной области, в качестве ТГ. и 
Г.„. были выбраны следующие значения: 

1 0 

1: =0,[: 

В качестве датчика псевдослучайных чисел предлагается исполь- 
зовать генератор псевдослучайных чисел (ГПСЧ), получивший наз- 
вание «вихрь Мерсенна» ((англ.) Мегѕеппеѓуіѕїег) [50]. 

Вихрь Мерсенна обладает следующими достоинствами: 

• большой период, равный 2'?”' - 1; 

• низкие временные затраты на генерацию псевдослучайных 
чисел и обновление внутреннего состояния ГПСЧ; 

• хорошие статистические свойства ГПСЧ, о чем свидетельст- 
вует прохождение генератором тестов ОІЕНАКР” [51]. 


2.2.3. Обучение 


Алгоритм самоорганизации нейронной карты 
по алгоритму Кохонена 


Алгоритм самоорганизации нейронной карты был предложен 
финским ученым Тойво Кохоненом в 1982-м году [52]. Алгоритм 
обучения нейронной карты Кохонена состоит из следующих 
этапов: 

1. Для каждого вектора Ве В выполняется: 

1) вектор Б подается на вход нейронной карты; 


2) нейронная карта определяет нейрон-победитель у 


еитег 


нейрон, вектор весов которого находится ближе всего ко 
входному вектору; 
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3) значения элементов векторов у 


о ? 


соответствующих ней- 


ронам карты, корректируются (2.9) в сторону значений 
компонент входного вектора В: 


Ув У, 6 +ПО(В, У, ); (2.9) 
К=ЬК; ПЕ(0;1] 
с = паюйроиг(\, у, с Е[0;1], (2.10) 
где: 7 — коэффициент, определяющий скорость обучения. На пер- 


утпег ); 


вых итерациях процесса обучения у принимается достаточно ма- 


лым и возрастает к последним итерациям обучения до |; с - 
коэффициент, определяющий степень соседства между нейроном- 
победителем и нейроном, вектор весов которого подвергается 
коррекции. 

Коэффициенты у и с позволяют уменьшить вероятность 


появления «пустых» нейронов в обучаемой нейронной карте, т.е. 
нейронов, слабо участвовавших в процессе обучения, в результате 
чего классы, соответствующие данным нейронам, оказываются 
разреженными. Для пространств, распределение векторов в кото- 
рых имеет ярко выраженные «точки притяжения» (Т.е. для статис- 
тически неоднородных пространств), появление таковых пустых 
векторов представляется существенной проблемой, вызывающей 
снижение эффективности применения к данным пространствам 
нейронных карт Кохонена. 
К функции пеейбоиг предъявляются следующие требования: 


• функция должна быть непрерывна и определена всюду над М, 
так как Ум; ис №; 

• чем ближе расположен нейрон-победитель к нейрону, вектор 
весов которого подвергается коррекции, тем ближе значение о 
должно быть к | (и достигать | в случае, если выполняется кор- 
рекция значений весов нейрона-победителя, чего требует правило 
Уэбба); 

2. К счетчику итераций цикла обучения прибавляется единица 
и, если счетчик не достиг максимального значения Г, 


тегайоп Койоптеп 1еасй ? 
выполняется переход на шаг 1). 

Рассмотренный алгоритм обладает следующими особеннос- 
ТЯМИ: 
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• необходимость выбора между нормализацией векторов или ее 
отсутствием. Проблема нормализации векторов состоит в выборе 
между выполнением и невыполнением нормализации векторов 
пространства Г перед обучением и функционированием нейрон- 
ной карты Кохонена. Нормализация векторов приводит к гаранти- 
рованному схождению алгоритма (т.е. рано или поздно на очеред- 
ной итерации цикла обучения ни один из векторов исходного 
пространства не сменит свой класс), но при этом коллинеарные 
векторы будут рассматриваться как эквивалентные без учета воз- 
можного существенного различия их модулей. В случае отсутствия 
нормализации векторов при вычислении расстояний между 
векторами будут учитываться также их модули, однако в данном 
случае сходимость алгоритма Кохонена не гарантирована. 

Выбор между выполнением и невыполнением нормализации 
векторов исходного пространства перед обучением и функциониро- 
ванием нейронной карты Кохонена производится, исходя из сущест- 
вующих практических соображений. В описываемом алгоритме сов- 
мещения не выполняется нормализация входных векторов нейронной 
карты Кохонена, а алгоритм ее обучения останавливается при дости- 
жении максимально допустимого количества итераций; 

фактор появления на карте «пустых» нейронов. 

Проблема «пустых» нейронов состоит в появлении на карте 
нейронов, не участвовавших в процессе обучения, из-за чего век- 
торы, соответствующие данным нейронам, оказываются в разре- 
женных областях исходного пространства векторов. 

Для минимизации негативного влияния данного фактора по 
рекомендации [18] в (2.9) используется коэффициент с. 


Расчет коэффициентов соседства нейронов 


Важным элементом алгоритма самоорганизации нейронной 
карты Кохонена является функция пеейбоиг (2.10). От выбора 
данной функции зависят скорость сходимости алгоритма и его 
вычислительная сложность. Существуют несколько подходов к 
определению функции иеірйроиғ ‚ среди них: 

Дельта-функция. В статье [24] в качестве функции пеейроиг 


была использована дельта-функция (2.11) в соответствии с принци- 
пом «победитель получает все» ((англ.) улппег (аКез а, МТА): 
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|; Уу = ТУЕ ) 
= (2.11) 
0; ии 


мчтпег* 


с = пеірћроиғ(у, у, „„„.) 

В данном случае корректируется только вектор весов нейрона- 
победителя, что соответствует принципу Уэбба обучения нейрон- 
ных сетей. 

Проблемой данного подхода является нерастяжимость карты. 
На очередной итерации обучения коррекции подвергается только 
центральный вектор класса, ассоциированный с нейроном-победи- 
телем, что нарушает связи данного нейрона с соседними нейрона- 
ми и снижает четкость границ между классами. Логично было бы 
предположить, что для повышения качества обучения карты необ- 
ходимо, чтобы карта была «растяжима», т.е. на очередной итерации 
алгоритма обучения коррекции подвергались бы все нейроны сети, 
а степень коррекции определялась мерой близости нейрона к ней- 
рону-победителю. 

Достоинством данного подхода является его малая временная 
СЛОЖНОСТЬ. 

Радиальная функция. Альтернативой использованию дельта- 
функции в качестве функции пеейБоиг является использование 


радиальной функции (2.12): 


уз 1 Ш АТТА 


СЫА ИШИНЕ. 1НРННЕГ Роа а 2 


О = пеівһроиғ(у, М ние) = 


© 


В данном случае достигается плавное уменьшение коэффи- 
циента соседства по мере удаления от нейрона-победителя и, сле- 
довательно, карта оказывается «растяжимой». 

Проблемой данного подхода является его существенная времен- 
ная сложность по сравнению с использованием дельта-функции. 


Выбор функции расчета 
коэффициентов соседства нейронов 


В ходе проведенных экспериментов использование радиальной 
функции в качестве функции пеейроиг при обучении нейронной 


карты Кохонена позволило получить большую точность совмеще- 
ния, чем использование дельта-функции. Существенная временная 
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сложность радиальной функции в значительной степени нивелиро- 
валась возможностью распараллеливания ее вычисления. 


2.3. Классификатор 


2.3.1. Структура, принципы обучения 
и функционирования классификатора 


Классификатор – это набор методов, алгоритмов и заранее оп- 
ределенных параметров алгоритмов, решающий задачу разбиения 
некоторого множества образов на набор непересекающихся под- 
множеств-классов. 

Классификатор описывается функцией с/аѕѕ/йсаіоп 


С = (аз; Тсапой(Н ,аеѕс); (2.13): 
Ней С =} се: ЧЕЙ 
Усс ЕС; с, пс, =, 
где: Н – множество образов Й таких, что йс В; С - результиру- 
ющее множество классов с,, где 4 — индекс класса; деѕс — функ- 
ция, описывающая алгоритм построения векторного пространства 


признаков образов: 
у = аезс(Й), (2.14) 


у = {у И ны Е у, Е К, 


где у — вещественный вектор признаков образа Й. 

Фактически, задача классификации сводится к задаче построе- 
ния регрессии С по Н на основе набора обучающих данных – 
обучающем множестве образов Н, и тестовом множестве об- 


Іеас 


разов Н, для которых результат классификации (множества 
а. иа...) уже известен. 


іеасһ 1651 


Функция с/аѕѕіЙйсайоп имеет вид: 
а = сІаѕѕійсаііоп(ћ, аеѕс). (2.15) 

Таким образом, задача классификации предполагает маркиров- 
ку входных образов целочисленными маркерами-индексами клас- 
сов, к которым относится каждый из образов. 

В качестве алгоритма построения векторного пространства 
признаков образов можно использовать алгоритм, основанный на 
использовании энергетических текстурных характеристик Лавса, 
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текстурных признаков Харалика и ИНС, реализующей анализ глав- 
ных компонентов. Данный алгоритм описан в $ 2.3.2. 

В качестве классификатора предлагается использовать слож- 
ный классификатор, состоящий из набора простых (бинарных) клас- 
сификаторов. Структура классификатора приведена на рис. 2.3. 

Принцип работы простого классификатора представляется как: 


р, ма, = СЇаѕѕійеғ(Н,, ); (2.16) 


пап 


Н„енН; Н, = Н, оН, ; т,пє\; 


тп тп + тп 
Оу Е (а: а мау Е 

где: т — порядковый номер уровня дерева, на котором расположен 

классификатор; и — порядковый номер классификатора на уровне 

т; Н, — подмножество образов классифицируемого множества 

образов Н , отнесенное классификатором на предыдущем уровне к 

одному из двух классов; Н, – подмножество образов, отнесен- 


тп 


тп + 
ных классификатором к классу, условно маркируемому как 1 (т.е. 
пеН, _ <> сіаѕѕійе(һ)=1); Н, _ – подмножество образов, отне- 


тп - тп - 


сенных классификатором к классу, условно маркируемому как —1 
(т.е. РЕН, . <> с1а55Тег(р) =-1); Р, – набор бинарных мар- 


тп + 
керов, определяющих класс образов лє Н 


тагу 


в смысле данного 


тп 


бинарного классификатора. 
Множества Н и Н передаются на следующий уровень 


тп + тп - 


классификатора по соответствующим ребрам дерева. 
Если узел (т, п) дерева является стоком, множество образов 


Н„„ формирует соответствующий класс с,. 
Обучающее множество образов Н 


іеасһ ? 


тестовое множество 


образов Н, и соответствующие им множества эталонных откли- 


1е5ї 


ков Чиа 


1651 


являются параметрами алгоритма совмещения. 


Существуют две различные версии алгоритма совмещения. 
Одна из версий использует в качестве простых классификаторов 
ИНС прямого распространения без обратных связей, другая версия 
использует машины опорных векторов. Глава 4 содержит рекомен- 
дации по использованию той или иной версии предлагаемого 
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алгоритма совмещения, выработанные по результатам эксперимен- 
тальных исследований. 


Н 


Рис. 2.3. Структура классификатора 


2.3.2. Алгоритм построения векторного пространства 
признаков образов 


Общее описание алгоритма построения 
векторного пространства признаков образов 


Алгоритм построения векторного пространства Г признаков 
образов лє Н описывается функцией деѕс (2.14). 


Алгоритм построения векторного пространства признаков 
образов состоит из двух этапов: 
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1) расчет вектора признаков у'є /' для каждого образа лє Н 
с помощью того или иного алгоритма описания образов, в качестве 
которых можно использовать: 

– алгоритм, основанный на использовании энергетических 
характеристик Лавса (2.17); 
– алгоритм, основанный на использовании текстурных 
признаков Харалика (2.18). 
у' = Іамѕ(А) (2.17) 
у' = йагайск(И) (2.18) 

Энергетические характеристики Лавса используются на верх- 
них уровнях дерева классификаторов, то есть на тех уровнях, на 
которых происходит разделение существенно отличных друг от 
друга классов, например, класса водных объектов от класса лесных 
насаждений. Гекстурные признаки Харалика используются на ниж- 
них уровнях дерева классификаторов, поскольку они хорошо заре- 
комендовали себя для описания текстур в задачах разделения клас- 
сов, если граница между классами является существенно нелиней- 
ной (например, в случае разделения класса хвойного леса и класса 
лиственного леса); 

2) понижение размерности векторного пространства У' 
признаков с использованием ИНС, реализующей анализ главных 
компонентов (РСА-сеть) (2.19). 

у = рѕа(у') (2.19) 

На данном этапе на базе векторного пространства Г“! приз- 
наков образов й с помощью преобразования Карунена-Лоева фор- 
мируется векторное пространство У. Расчет матрицы преобразова- 
ния Карунена-Лоева выполняет РСА-сеть. Понижение размернос- 
ти векторного пространства /' необходимо для улучшения вре- 
менных характеристик и повышения показателей качества работы 
классификатора за счет перераспределения информации в вектор- 
ном пространстве Г. 


Алгоритмы описания образов 


Алгоритм описания образов, основанный 
на использовании энергетических характеристик Лавса 


К.И. Лавс в [36] предложил использовать для описания текстур 
специальные энергетические характеристики. Алгоритм описания 
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образа й с помощью энергетических характеристик Лавса состоит 
из следующих этапов: 
1) формирование набора пространственных фильтров Ё, ; 


Е, = 25: 
Е = {Е Е хє {1,Е,5,0, К}; ує{1,Е,5,0, К); 


еп еп х у 
аре ое, аен 
Л г = 34;6;4;}; Ј/, к =1-1-50;231; 
Ли; ={—150;2;0;-1}; Л, и ={-12;0;-2;1); 
в =;-4;6;—4;}, 
где векторы Л ,, /,, ғ, Лис» Лии» Л в Предполагают выделение 
отдельных характеристик текстуры: 
• уровень спектральной яркости (вектор /., , ); 
• границы элементов текстуры (вектор /,, ;); 
• отдельные пятна на текстуре (вектор /., ; ); 
• периодически повторяющиеся перепады яркости («волны»; 
вектор Хи }; 
• текстурную «рябь» (вектор #,, ,; «рябь» – ситуация, когда на 


текстуре повторяются с некоторой периодичностью отдельные, со- 
ставляющие ее, элементы). 


Комбинирование векторов /, ,, ие, Лис» /.„ ж» Ли в ПОЗВО- 
ляет получать отдельные характеристики текстуры, например, 
фильтр Е, ,„ позволяет получить информацию о насыщенности 


вертикальных краев отдельных элементов текстуры. Каждый фильтр 
Е„ х > Кроме Ё, ‚, , имеет нулевую сумму составляющих его весов; 


еп ху 


2) применение к образу й набора фильтров Е 


еп > 


по результатам 
чего составляется набор матриц /,, : 
и = 40, 5 (2.20) 
хєЄ{І,Е,5,0,К}; ує{1,Е, <, ,К};; 
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3) составление по набору матриц й, набора энергетических 


карт О. и набор матриц РО, , содержащих характеристики текс- 
тур, покрывающих образ й: 


(+, /+и’) р, (і, 7) 
п ху | 


Би (1) = (2.21) 


(1.7) (1-а, и): ей й И. (1, Ј) 
О. =20+; О, >=5; (2.22) 
В = р. = {2 В. 


Іам 


р ху — Р.» ху + Дея ж› Хх у 


еп 
—^ 
р, Ва я 22 5 


і=1, 7 ј=1,7; 2Е{ЁЕ‚5,Й’, В}, 
где 62 — размер окна, в котором рассчитываются энергетические 
признаки (выбирается оператором). Фильтр Е, ,, является норми- 


еп [ 


рующим фильтром – с его помощью из описания текстур удаляется 
влияние не равномерности освещения на спутниковом снимке; 


4) ассоциирование с каждым пикселем 6,, єй вектора энерге- 


тических характеристик 4, ,: 


а, , = {Р,, (Ё 1} М2 р. 2 (1, 1}; (2.23) 
5) расчет вектора энергетических характеристик у' образа Й: 
У а 


Ге) 
Ь ЕЙ 
=. (2.24) 
А| 


Алгоритм описания образов, основанный 
на использовании текстурных признаков Харалика 


Р.М. Харалик в [38] предложил описывать текстуру характе- 
ристиками специальных матриц вхождений Н , рассчитываемых 
по образу А. Каждая матрица вхождений описывает структуру 
текстуры в соответствующем направлении от одного края обраба- 
тываемого образа к другому. По результатам расчета характе- 
ристик матриц Н составляется вектор текстурных признаков у' 
образа И ‚ которые могут рассматриваться как характеристики 
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текстуры, покрывающей каждый пиксель образа. Алгоритм расчета 
текстурных признаков Харалика у, для К -го спектрального канала 
образа И состоит из следующих этапов. 

1. Выбор нескольких направлений о = {а}; х= 1,Х; а, є[0,2л), 
где ©, суть есть угол, описывающий направление соответствую- 
щей матрицы вхождения Н. В [37, 38] выполняется расчет матриц 


а А лл Зл 
вхождений для четырех направлений: о = УТРЕ 
2. Выбор размера О, области соседства каждого пикселя обра- 
за. В [38] было выбрано ©=3. 


3. Заполнение матриц вхождений Н = {Н}; х=1,Х размером 


А на А элементов нулями, где А есть количество уровней кван- 
тования диапазона спектральных яркостей. Так как первоначально 
спектральные каналы спутникового снимка В являются, чаще 
всего, 8-ми битными без знаковыми целыми, то А = 256. 

4. выполнение для каждой пары «пиксель», «матрица 


вхождений»: (5,,, Н,); 6,, ЕЙ следующих операций: 
4.1. «Наложение» на пиксель 65,, окна размером ©, на ©), 
пикселей. 
4.2. Расчет координат (1,,/,) крайнего пикселя 6,, по 
направлению а, : 


* для с, =0; 0, =3: 1 =; Ј, = ј, +1; 


для а =—; О =3: 1 = 1; Ј = ј, +1; 


для 2, =—; 6, =3:1.=1- | ДЕЛ; 


• для СЕ О Зет уу 


4.3. Если координаты (і, , /„. ) крайнего пикселя по направ- 
лению а, выходят за пределы изображения В, то выпол- 


няется переход к следующей паре «пиксель», «матрица 
вхождений». 
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4.4. Расчет: 
Н (6, 


ЛК? ЯР 


і) = Н, (0, „55, а) +1. (2.25) 


5. Масштабирование матриц вхождений Н, к сумме своих 
элементов, равной единице, для чего выполняется преобразование: 


Н. (б, М бк) 


рхи 2 оН, (т, п) 
ИЖ 


6. Расчет (2.27) для каждой матрицы Н, характеристик Н х(р), 


(2.26) 


«(В рк»б Т = 


И И Ш, ы — средние по строкам и столбцам соответственно, 


С, „И С, „м — дисперсии по строкам и столбцам соответственно: 


Н.(р) = У У и РЕМ (2.27) 


т=1{ п=| 
А 


П =>", (т,п); 


т=] п=] 
И, соі - Уну нти), 
п=| т=] 


ДА 


О. „ом = У (т- И, в ин, (т, п); 


т=| 


А 
О реа = У (п и- 0, „УН (тп), 


п=] т=1 


где р — коэффициент, зависящий от гистограммы образа. Для 
б А 
упрощения А вычислений р Е =127,5. 
7. Расчет для каждой матрицы Н, следующих характеристик: 


Уту, (т,п))>; (2.28) 


А-1 А А 


ћ, = У; Н,((т,п)); (2.29) 


| п=] 
| | 


т-п\= р 
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А А 
У У отин, (тп) -— „сы. 
р) 


= (2.30) 
ОО д 
А А 
= У тУ (т-и, „Н, (т,п), (2.31) 
т=1 п=] 
У 2.32 
А. = Н (т,п); | 
аз У Эчис „От, п) (2.32) 
Р, - п.т) (2.33) 
2л д ЕЕ, 
1=2, тУ (т-ћ,,) Н. (т); (2.34) 
т=2 п=] 
А А 
а= У У Н, (т,п)1ов(Н,(т,п)); (2.35) 
т=| п=] 
2А о Ре 
һо = У Н.(т)іор(Н.(т)+е); е<10“, (2.36) 
т=2 


где: й, — однородность; й,, — контраст; ћ,, – корреляция; №,, 
дисперсия; /,; — гомогенность; й, — среднее по суммам; й; – 
дисперсия по суммам; й, — энтропия; й. — энтропия по суммам 
(коэффициент е устраняет неопределенность логарифма при 
Н.(т)=0; 

8. Составление для каждого пикселя 0, ;; ,; ей вектора приз- 


наков 4, {И}; х=1,Х; у =1,9. 


9. Расчет вектора признаков у, для К -го спектрального канала 


образа И: 
Уа, 
' Ь ЕЙ 
у ==. 2.87 
ЛИТІ (2.37) 


Наконец, вектор признаков у' образа Й рассчитывается как 


К ' 
У" =], . 
к=] 


Глава 2 65 


Алгоритм понижения размерности векторного 
пространства признаков, основанный 

на использовании ИНС, 

реализующей анализ главных компонент 


Структура, принципы обучения и функционирования 


Анализ главных компонент ((англ.) ргіпсіра! сотропепі апа]у- 
515, РСА) – является способом понижения размерности векторного 
пространства с регулируемой потерей информации. 

Пусть имеется векторное пространство Г'={у}; |Г'Е=О’'. 


Необходимо рассчитать однозначное преобразование И' — 7, 
где Г - векторное пространство такое, что Г = {у}; 7|ЕО<О'. 
Преобразование рѕа (2.19) должно выполняться с учетом возмож- 


ных потерь информации, поэтому необходимо наличие некоторого 
порогового значения Г,„, Е [0,1), определяющего минимальный 


объем информации, который должен «содержаться» в Г изГИ'. 
Величину Т = 1007, можно рассматривать как процентную 


пјо рғос 
долю информации из /И', сохраняемой в И. 
В качестве функции рѕа анализ главных компонент исполь- 


зует преобразование Карунена-Лоева [39, 40], матрицу которого 
можно рассчитать с помощью специальной ИНС-РСА-сети. Таким 
образом, функцию рѕа можно записать в виде: 


у = рѕа(у') =Й"'у'; (2.38) 
У =[0. ]; т=1,0', п=1,0', (2.39) 


тп 
где И’ – матрица преобразования Карунена-Лоева, являющаяся 
также матрицей весов нейронов РСА-сети. 


Функция рѕа описывает процесс работы РСА-сети матричное 


умножение И’'у’ соответствует расчету взвешенных сумм входов 
нейронов сети (векторы весов которых являются строками матрицы 
У ), а функцией активации в данном случае является прямая у= х. 
Алгоритм обучения РСА-сети состоит из следующих этапов: 
1. Центрирование векторного пространства /' (2.40) с получе- 
нием на выходе векторного пространства И: 


сешег * 
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Я Е Ота } | Усешеғ = ОРЕ Е, У сешег 42° У ешег о) > (2 .40) 
В ' 
У ть, = У— МУ. 


2. Расчет вектора весов у первого нейрона РСА-сети с помо- 


Таким образом, Му 


сещег 


щью правила Ойя [53], являющегося модифицированной версией 
правила обучения Уэбба: 


2.1. Установ компонентов вектора у в нулевое значение. 


2.2. Пересчет вектора у, (2.41) для каждого у, ЕТ: 


сешег сешег 


(2.41) 
где у коэффициент скорости обучения, Е 7(0,1]. 
2.3. Повторение пункта 2.2 алгоритма для различных зна- 
чений у. 
По результатам экспериментальных исследований в качестве 
стратегии изменения у было выбрано постепенное увеличение 7 с 


шагом 0,05 от 7, = 0,05 до 7,„, =1. 


3. Расчет вектора весов у, т-го нейрона РСА-сети с помо- 
щью правила Сенгера [54]: 
3.1. Установ компонентов вектора у, в нулевое значение. 


3.2. Пересчет вектора у, (2.42) для каждого у, єҮ,,„,„: 


сешег сешег ° 


Ут = "я ИУ т он ); (2.42) 
и и (2.43) 


т сепіеғ ? 


т-1 


У сетер т Е У етеу т уз У й У, ? 


ћ=] 
где 7 — коэффициент скорости обучения такой, что 7 є (0,1]. 
3.3. Повторение пункта 3.2 алгоритма для различных значений 
7 ‚ способ выбора которых приведен в пункте 2.3 алгоритма. 


4. Повторение пункта 3 для т=2,О". 
5. Удаление из нейронной сети нескольких нейронов (и соот- 


ветственно, удаление из матрицы И’ нескольких строк). Векторы 
и, являются собственными векторами для матрицы автокорре- 


ляции векторного пространства Г 


сешег ? 


причем | у, |= 1. Для каждо- 
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го вектора у, рассчитывается соответствующее собственное число 


п 


еівеп, , после чего вектор собственных чисел еіреи = {еіреп, } нор- 
мируется. 

Строки у, в матрице И’ изначально упорядочены по убы- 
ванию их собственных чисел. Нормированные собственные числа 
еіреп можно рассматривать как относительный объем информации 
из векторного пространства /', «сохраняемый» собственными век- 
торами у, в процессе преобразования Карунена-Лоева. Таким об- 


разом, в И должны быть оставлены только первые М строк, где 
М такое, что выполняется условие: 


т 


) еівеп, > Ти (2.44) 


т=] 
Алгоритм обучения РСА-сети сохраняет в пространстве / как 
минимум ТГ % информации, имеющейся в /', в смысле 


тр рғос 
собственных векторов и собственных чисел матрицы автокорре- 
ляции векторного пространства Г'. 


2.3.3. Машина опорных векторов 
Структура, принципы обучения и функционирования 


Машина опорных векторов ((англ.) зиррой уесюг тасһіпе, ЗУМ) 
выполняет линейное разделение пространства образов на два класса. 
Пусть дано пространство Г (2.46), состоящее из векторов у: 


у= {у,...,У,, „у%0}5 9=1,0; 
у =є К 


В контексте решаемой задачи пространство Г образуется из 
векторов признаков классифицируемых образов. 

Ключевая особенность машины опорных векторов состоит в ее 
способности выполнять нелинейное разбиение классифицируемого 
векторного пространства Г путем отображения Г в векторное 
пространство /' такое, что |Г У '|, и линейное разделение кото- 


рого возможно или производится с меньшей ошибкой, чем линей- 
ное разделение исходного пространства. 
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Векторное пространство Г" (2.46) формируется следующим 

образом: 
И '= {у}; (2.46) 
УЕ,» 0} 9 =1,0"; 
! 
у' ЕК. 
ф:И ЭИ" 
еа 
у'= ф(у). 

Разделяющая поверхность / рассматривается в пространстве 
Г'. Разделение, таким образом, является линейным в отношении 
Г'. Требуемая нелинейность разделения пространства Г обеспе- 
чивается отображением ф . Размерности пространств Г иТ' не- 
обязательно совпадают и в общем случае имеет место быть сле- 
дующее соотношение: О < 0'. 


Принцип работы машины опорных векторов проиллюстриро- 
ван рис. 2.4. 


Нелинейное 
отображение 
В. В Ез 


Обученис 


машины 
опорных 
Обучающая выборка векторов 


6 Везделяющая 
х поверхность 


Рис. 2.4. Принцип работы машины опорных векторов 
Машина опорных векторов в процессе своей работы не ис- 
пользует отображение ф — ЗУМ оперирует ядрами скалярных про- 


изведений векторов пространства /'. Уравнение разделяющей по- 
верхности Ё, описываемой машиной опорных векторов, имеет вид: 


2. 
1(8,у) = У иКегпе[(з,,у)+Т; (2.47) 


1=1 


Кегпе($,„у) = ф(5)ф(у); 
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5 = {5,}; 1=1,0°%; 5сЙ;; 
Й ={и}; жм є; 
ТєК, 
где: Г — разделяющая поверхность; 5 — множество опорных векторов; 
И - вектор коэффициентов при опорных векторах; кеғие/(у,,у,) — 
ядро скалярного произведения — скалярное произведение образов 
векторов 5 и у в пространстве /'; Т – пороговое значение. 


В качестве ядра скалярного произведения могут быть исполь- 
зованы функции, удовлетворяющие теореме Мерсера. 


Пусть Кеғие/(у,,у,) — непрерывное симметричное ядро, опреде- 
ленное на интервале у,,у, є[а, В]. 


Ядро Кегие[ описывается рядом (2.48): 
КегпеКу,,у,)= х; (у, )6, (у;)6, (У, ); (2.48) 
к=] 


А, > 0. 


Для обеспечения корректности этого разложения и его абсо- 
лютной и равномерной сходимости необходимо и достаточно 
выполнение условия (2.49): 


| Г етецу,у,)е (у) (уа 210; (2.49) 
МС [000 <. 


Доказательство теоремы приведено в [55]. 

Бинарная кластеризация, выполняемая машиной опорных век- 
торов, состоит в том, чтобы рассчитать значение /(5,у) для неко- 
торого входного вектора у. Если [(5,у) > 0, то вектор у относится к 
классу, условно маркируемому +1, в противном случае [(5, у) <0 – 
вектор у относится к классу, условно маркируемому —1. Ситуация 
І(5,у) = 0 рассматривается как некорректная, вектор у маркируется 
как «неклассифицированный», а образ, соответствующий вектору, в 
дальнейшем не рассматривается. 

Таким образом, цель обучения машины векторов состоит в 
определении следующих ее параметров: 

• множества 5 опорных векторов; 
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• вектора И’ весов при опорных векторах; 
порогового значения Г. 


Ядро скалярного произведения векторов 


Существуют следующие ядра Кегие[ скалярного произведения 
векторов: 


• полиномы натуральной степени Кегие/„„, (2.50); 
• радиальная базисная функция Кегие! „, (2.51); 


• гиперболический тангенс Кегпе[, (2.52). 
КеғпеІ (У.У) = (уу, +12; рє№; (2.50) 


Кегпе/[ „(у ,у,)= |515 м1} СЕВ*; (2.51) 
с” 


Кегпе[, (у, у.) = (ругу, +); р,Кє К, (2.52) 
где р,с, р, к – параметры соответствующих функций, выбирае- 
мые оператором. 

Экспериментальные исследования показали, что выбор ра- 
диальной базисной функции Кегпе[ „„ в качестве ядра скалярного 
произведения приводит в конечном счете к более высокой точ- 
ности совмещения по сравнению с полиномами Кегие/ „, и гипер- 


болическим тангенсом Кегпе/, 
Обучение 
Цели и задачи обучения 


Как было показано в [17], задача обучения машины опорных 
векторов есть задача условной квадратичной оптимизации. Необ- 
ходимо максимизировать функцию Лагранжа О при выполнении 


условий Каруша- о е ыы 54) и (2.55) [56]. 


О(а) = Уа -15 аа а,4 кетеКу, ,у’); (2.53) 


1=1 =] 
а, =0; (2.54) 


0<а <С,,; (2.55) 
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У ЕЙ; ИСИ; 
а Е {-1;+1}; 
С є, 


5ит 


где: /’- обучающее множество векторов; О” — мощность обучаю- 
о * 
щего множества; 2 = {4}; і= 1,0 - эталонные отклики для векто- 


ров из обучающего множества; 4А = {а}; і=1,0° – множество мно- 
жителей Лагранжа, по которым выполняется оптимизация функции 
О; С. – некоторая положительная константа — параметр алго- 
ритма, 

Вектор коэффициентов Й’ при опорных векторах рассчитыва- 
ется по множителям Лагранжа А и эталонным откликам Р; 

и = а,1,. (2.56) 

В процессе обучения большинство множителей Лагранжа при- 
нимает значение 0, и таким образом, во множество 5 опорных 
векторов входит лишь часть векторов, составляющих обучающее 
множество: 5 с У”. Множество 5 опорных векторов представляет 
собой, в конечном счете, набор векторов из обучающего множест- 
ва, ближе всего лежащих к поверхности Ё, разделяющей прост- 
ранство на два класса. 

Пороговое значение Т рассчитывается (2.57) следующим 
способом: 

3. У и’ Кегие[(5, ,у')– а 
Т Ету (2.57) 

Для решения задачи условной оптимизации функции Лагранжа 

О используются различные алгоритмы, сравнение которых было 


выполнено в [41]. В контексте задачи совмещения лучшими 
показателями по времении качеству результата обладает алгоритм 
последовательной минимальной оптимизации ((англ.) ѕедиепіа! 
шипа! орітіғхайоп, ЗМО) [57]. 


Алгоритм последовательной минимальной оптимизации 


Алгоритм последовательной минимальной оптимизации осно- 
ван на идее разбиения задачи оптимизации целевой функции О по 
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набору переменных А на несколько задач оптимизации по под- 
наборам переменных. При этом имеет место сходимость алгоритма 


к глобальному максимуму целевой функции О’... 
Алгоритм $МО состоит из следующих этапов: 


1) формирование обучающего множества И” векторов и мно- 
жества 2) соответствующих им эталонных откликов; 

2) установ значений элементов множества множителей Лагран- 
жа А в0; 

3) инициализация дополнительных переменных, используе- 
мых, в том числе, для улучшения временных характеристик алго- 
ритма за счет введения дополнительных кэшей в оперативной 
памяти вычислительного устройства; 

4) выбор пары множителей Лагранжа; 

5) максимизация функции О по выбранным множителям 


Лагранжа; 
6) расчет новых значений элементов кэшей; 
7) проверка на выполнение критерия останова алгоритма: 
— если критерий останова выполнен, выполняется формиро- 
вание множества 5 опорных векторов (куда попадают толь- 


ко те векторы из обучающего множества у; , для которых 
выполняется условие а, * 0 ), множества коофициентов И 


при опорных векторах, а также вычисляется значение порога 
Т, после чего происходит останов алгоритма обучения; 
– если критерий останова не выполнен, то выполняется пе- 
реход к этапу 4 данного алгоритма. 

Алгоритм 5МО останавливается в следующих случаях: 

° превышено максимально допустимое количество итераций 
алгоритма; 

• выполнены условия останова, специфичные для используе- 
мого алгоритма выбора пар множителей Лагранжа. 

Для выбора пар множителей Лагранжа для очередной итерации 
алгоритма последовательной минимальной оптимизации предлага- 
ется использовать алгоритм Фана-Чена-Линя [58], предложенный 
группой тайваньских ученых: Коп?-Еп Еап, Раі-Нѕиеп Сһеп, Сһ1һ- 
Јеп [лп. 
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Алгоритм Фана-Чена-Линя базируется на использовании ин- 
формации о второй производной функции ошибки Е 


5ит ° 
Е (У) = МУ",У)- а. (2.58) 
Первый множитель Лагранжа а, в паре множителей, кандида- 
тов на ен выбирается последующему принципу: 


Тр =4# (а <С„ Уа, =По(а, > 00а, = –1}; (2.59) 
іє тах = {— АУЕ(у, )|1ЕТ,,}. 
1=1,0` 


Второй множитель Лагранжа а, в паре множителей – кандида- 


тов на чеч выбирается по следующему принципу: 
= (| (а <С.„ Уа = -1) (а > 00а = 1}; (2.60) 


ре 


ра Мп = 4 |1 ЕТ, 
1.0° (94 


с = Кегие(у, ,у; ) + КегпеКу' у’) 2КегпеКу, ,у’)а,а 


‚Аг 


С; с> 0; 
7 1072; с<0; 
у= а УЕ, „(у )+4 УЕ, (у); 
в В=у; у>0; 
В= 0; у<0. 


Благодаря С достигается максимальный «разброс» между па- 
рой множителей по второй производной, что приводит к максими- 
зации эффекта от шага оптимизации. 

Алгоритм $МО, использующий для выбора пары множителей 
Лагранжа на очередной итерации алгоритм Фана-Чена-Линя, оста- 
навливается также в случае, если выполняется условие: 


АУЕ(у,) +А.УЕ(у,) <е; (2.61) 


еЕ В *, 


где е — параметр алгоритма, выбираемый оператором. 
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2.3.4. Искусственная нейронная сеть прямого 
распространения без обратных связей 


Структура, принципы обучения и функционирования 


В качестве бинарного классификатора с/аз5 Тег (2.16) наравне 


с машиной опорных векторов была использована ИНС прямого 
распространения без обратных связей [22]. 

Функционирование ИНС прямого распространения без обрат- 
ных связей можно представить как: 


И. а<0; 
Ыпагу — 1; а >= 0; 
а= } (у) = Ј) (аеѕс(ћ)), 


где функция {ў описывает ИНС. Каждый нейрон сети описывается 


(2.62) 


в виде функции иеиғо : 
у, = пеиго, (у) = Г. (8,)), (2.63) 


Р =1,Н,; М. м | № = {у}; 
Ут > Мт э У, Є К, 
где: г — индекс (порядковый номер) нейрона; У, — количество 


внутренних нейронов в сети (настраиваемый параметр алгоритма); 
у, – входной вектор нейрона; и — вектор весов нейрона размером 


М, элементов; И’ — множество векторов весов нейронов ИНС; у, — 


выходное значение нейрона (уровень возбуждения нейрона); ў, — 


функция активации нейронов сети. 
Структура связей внутренних нейронов ИНС прямого распрост- 
ранения без обратных связей представляется в виде ориентиро- 


ванного графа без циклов С ‚ где ребро графа описывает синапс 
ИНС, вес ребра графа определяет вес синапса, узел графа описывает 
нейрон ИНС и предполагает применение к взвешенной сумме 
возуждений синапсов, входящих в узел, функции активации. Пример 
представления ИНС прямого распространения без обратных связей в 
виде ориентированного графа без циклов приведен на рис. 2.5. 
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Частным случаем рассматриваемой ИНС является многослой- 
ный перцептрон, однако, в отличие от последнего, ИНС прямого 
распространения без обратных связей не обязательно является 
полно связной (и, следовательно, четко структурируемой на нес- 
колько слоев), что улучшает качество обучения сети в виду допол- 
нительной нелинейности, вызванной отсутствием жестких требова- 
ний к структуре сети. 


Внутренние 
нейроны 


Выходной 
нейрон 


Рис. 2.5. Ориентированный граф внутренних нейронов 
ИНС прямого распространения без обратных связей 


Граф С описывается с помощью матрицы смежности С: 


811 8, 51н, 
С == 8,1 5, 8,н, (2.64) 
1 105 нур тт А 
8640,1}, (2.65) 


где элемент =, кодирует наличие (1) или отсутствие (0) ребра, 


выходящего из вершины г и заходящего в вершину р, т.е. синап- 
са, исходящего из нейрона г и заходящего в нейрон р. 
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Особым подвидом внутренних нейронов сети являются нейро- 
ны, соответствующие истокам графа С . На вход данных нейронов 
поступают не выходы других нейронов, а входной вектор у сети. 
Для данных нейронов М =О=у|. 

Кроме внутренних нейронов ИНС располагает также единст- 
венным выходным нейроном. На вход данного нейрона поступают 
выходы нейронов, соответствующих стокам графа С . Выходной 
нейрон рассчитывает среднее арифметическое своих входов, кото- 
рое и используется в качестве выходного значения сети. 

Алгоритм работы с ИНС прямого распространения без обрат- 
ных связей состоит из следующих этапов: 

1) выбор структуры ИНС (матрицы С ) с помощью грамматик 
графовой генерации Китано и генетического алгоритма ($ 2.3.4); 

2) инициализация векторов весов у нейронов ИНС с помо- 


щью алгоритма Нгуэна-Видроу, описанного в [59]; 
3) обучение внутренних нейронов ИНС в пакетном режиме; 
4) применение ИНС для решения практической задачи. 


Функция активации нейронов 


В качестве функции активации ў 


асі 
следующие функции [17]: 
логистическая функция: 


могут быть использованы 


] 


Ја 6) = Труа (2.66) 
те 


–тіор ё » 

тіовє В”. (2.67) 
Логистическая функция имеет вид, приведенный на рис. 2.6 
(верхняя кривая; т1ое = 10). К недостаткам логистической функции 


относятся ее асимметричность (устраняется вводом дополнитель- 
ного вычитания 0,5 в 2.66), сложность поиска подходящего значе- 
ния т1ор и большой диапазон возможных значений тор; 


• гиперболический тангенс: 
ГОТ (5) = (ат (т, 26) (2.68) 
Е = (2.69) 


Гиперболический тангенс имеет вид, приведенный на рис. 2.6 
(нижняя кривая; т, =1, т, , = 2). К достоинствам гиперболичес- 
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кого тангенса относятся симметричность и простота расчета произ- 
водной. К недостаткам – несколько настраиваемых параметров. 

В качестве функции активации можно выбрать гиперболи- 
ческий тангенс (2.68) за счет его симметричности и простоты 
расчета производной. 


—— Логистическая функция 
— Гиперболический тангенс 


Рис. 2.6. Функции активации 


Обучение 


ИНС прямого распространения без обратных связей относится 
к классу ИНС, обучаемых с учителем. 
Обучение ИНС осуществляется по множеству Н, образов, 


для которых известны эталонные отклики а 


геасй * 
Существуют несколько подходов к обучению ИНС: 
• индивидуальный подход, в котором коррекция векторов весов 


у нейронов сети выполняется для каждого образа й, є Н, ; 


• пакетный подход, в котором коррекция векторов весов у 
нейронов сети выполняется единожды для каждой итерации прого- 
на сети над всем множеством Н, образов после завершения 


теас 
прогона. 

В предлагаемом алгоритме совмещения для обучения исполь- 
зуется пакетная стратегия и алгоритм обучения состоит, таким 
образом, из следующих этапов: 
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1. Расчет у, 
21-Е 


геасћ = 0 А 


3. рта, = 0, где ртаа, – градиент Е, „(Ў ). 


= деѕс(й 


(еасй 


) для каждого Й ‚ЕН 


еас геасй * 


4. Повторение для каждого Й „ЕН 


теас 


4.1. Расчет а = У (У). 
4. 2. Е геасћ — Е аср га (бо а) ? . 
4.3. ртай, = етай, + етаа 


ситеп Е? 


‚ следующих действий: 


геас 


где егаа 


ситет Е 


– градиент 
Е „(Ў ) на текущем прогоне ИНС. 
и Н 


ситешт Е 


Приращения ұғаа рассчитываются с 


ситет Е 
помощью алгоритма обратного распространения ошибки. 
Е 


_ Чеасй __ 


е2 


геасћ — Зас 
2 Н іеасћ | 


б. рай, ЕЕН 


| Е 
7. Коррекция весов нейронов ИНС с использованием рға, с 
помощью того или иного алгоритма коррекции весов. 
По завершении обучения сети выполнятся тестирование ка- 
чества обучения, для чего рассчитывается Е, по тестовому мно- 


тез; 


жеству Е, образов, для которых известны идеальные отклики 4 


1е51 


>> Е„„|. Тестирование качества обуче- 


ния применяется в процессе выбора структуры ИНС и выбора па- 
раметров обучения ИНС. 


Основным требованием к Е, и Е, является необходимость 
репрезентативности данных множеств. 


ИНС. Как правило, | Е 


геасй 


Алгоритм обратного распространения ошибки 


Рассмотрим нейрон г , которому соответствует вектор весов у. 
Для данного нейрона необходимо рассчитать те компоненты 
вектора 2таа и матрицы Н которые соответствуют 


ситет Е сштеш Е » 


компонентам вектора у, в Ў. 
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Пусть вектор А, ={",} содержит индексы нейронов ИНС, вы- 
ходы которых поступают на вход нейрона (| К, |= М). Пусть век- 
тор К, = {/,„} содержит индексы нейронов ИНС, на вход которых 


поступает выход рассматриваемого нейрона (| А, |= М). При 


ош 


этом предполагается, что компоненты вектора ғр” аа и матри- 


ситеш Е 
цы Н. ғ» Соответствующие нейронам из А, уже рассчитаны 
(обратное распространение ошибки расчет градиента и гессиана от 


стоков графа С к его истокам). 
Для нейрона г необходимо рассчитать: 


у 27аЯ „мен Е т? т = | М, 
ду д 
тайы в а 202-02 0з. 0.70) 
ду, бу г с, дм, 
• атаа,; т=1,М, 
огай, __5Е _ дЕ ду, д8, (2.71) 
А От бу, Ё, 9и, 
В (2.70) и (2.71): ==. 
5 - ОЕ 
Для нейронов, соответствующих стокам графа С а = 
= а га = У іеасћ 


‚ где Н | - количество нейронов, соответствующих 


ош 


Н 


ош 


стокам графа С , то есть фактически частной производной функции 


—^ 


ошибки по выходу стока графа С, является отклонение выхода 
сети от идеального выхода, нормированное по количеству стоков 


графа С. 
о. ағай, 


ош є Кош 


я прочих нейронов: —— = ж, 
Для пр р бу, м 


оиі 


ОЕ р 
Е, (СА) 
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СЧА 
би ” 
‚6 _ 
ге тт 


Алгоритмы коррекции весов ИНС 


Для коррекции весов ИНС могут быть использованы следую- 
щие алгоритмы: 

• алгоритм градиентного спуска. 

Использует для коррекции весов только рғаа; . 


Пусть Ж, — оптимальное множество весов сети (т.е. множество, 


дающее глобальный минимум Ё,(Й)). Тогда имеет место быть 
следующее разложение в ряд Тейлора (2.72) в предположении, что 
Е (Ў ) имеет сколь угодно много производных (среднеквадра- 
тичное отклонение, используемое в качестве Е, „„(Й), удовлетво- 
ряет этому условию): 

ай Ио) = Е 


па) + етаа. ди + К, (2.72) 
где К, — остаточный член ряда Тейлора. Процесс обучения пред- 
ставляет собой итерационный процесс постепенного приближения 
к Й,, поэтому (2.72) сводится к: 

Е „КЎ ) = Е) + втаа ОЙ? + К.. (2.73) 


пех 
Для сходимости алгоритма к Ў, необходимо, чтобы на каждой 


геас 


итерации процесса обучения выполнялось условие Е, „(И )> 
>Е „(И ), для чего достаточно выбрать 60 =-ета4,; , тогда 


—вта4; <0 всегда, и Е’ „(Й) будет гарантировано убывать. Таким 
образом, алгоритм градиентного спуска на каждой итерации своей 
работы выполняет коррекцию Ж, = И -п,етай,, где п, є(0,1] – 


пех! 


коэффициент скорости обучения, зависящий от номера итерации #; 

• метод Гаусса-Ньютона, квазиньютоновские методы. 

Данное семейство алгоритмов использует для коррекции весов 
нейронов не только градиент ета, , но также и гессиан Н; (т.е. 


информацию о второй производной функции ошибки). Очевидно, 
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что для среднеквадратичного отклонения, используемого в качест- 
ве Ё Н, будет всегда положительной квадратной формой. 


ѓеасћ ? 


Тогда: 


К, = УР" Н,бТ7 + Р; 


Е И-И) = Е и (И) + втаа 0’ ИОД + №,. (2.74) 
Очевидно, что при И – Ў, : етаа; = и Н, =0, тогда БИ можно 


| рассчитать, исходя из утверждения руаа.дЎ + 26 И' Н.Б =0, 
что 2 можно привести к виду: 
Ртаа, + б'н, = 0), (2.75) 


предполагая, что на очередной итерации алгоритма ОИ + 0. Тогда 
5 = -Н;!рғаа, . 

Таким образом, метод Гаусса-Ньютона на каждой итерации 
выполняет коррекцию Ў, = И – 1, Н; ғай, . 


пем 
Расчет гессиана Н, представляется вычислительно сложной 
задачей, поэтому существует ряд методов [27], апроксимирующих 
гессиан Н, некоторой формой и получивших название квазинью- 
тоновских методов. 
Одним из таких методов является метод Левенберга-Марк- 
вардта, который использует для расчета приближения гессиана Н, 


якобиан „Л; и вектор коэффициентов, рассчитываемых по собст- 
венным числам „Л; ; 


• эвристические методы. 

Существует ряд эвристических подходов к коррекции весов 
ИНС прямого распространения без обратных связей. К таким 
методам относятся: 

• метод сопряженных градиентов; 

• алгоритм КРКОР; 

• алгоритм Ои1сКРКОР; 

другие алгоритмы. 
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В контексте поставленной задачи совмещения данных ДЗЗ и 
ЦКМ для обучения ИНС прямого распространения без обратных 
связей был использован алгоритм градиентного спуска, поскольку 
данный алгоритм обеспечивает высокую производительность и 
приемлемое качество [особенно в сочетании с эвристическими 
алгоритмами инициализации весов нейронов ИНС, что су- 
щественно снижает вероятность схождения процесса обучения к 
локальному минимуму, позволяя достичь при этом глобального 


минимума функции ошибки Е „(Ж)]. 


Кроме того, градиентный спуск, в отличие от методов, осно- 
ванных на использовании информации о второй производной функ- 


ции ошибки ЕЁ, ,(Й’) (метод Гаусса-Ньютона, квазиньютоновские 


еасй 
методы), менее требователен к точности представления чисел, тог- 
да как означенные методы при расчете гессиана или приближения 
к нему требуют высокоточного представления вещественных чисел 
(лучше, чем двойная точность, то есть, чем максимальная точность, 
реализуемая аппаратно на современных вычислительных устройст- 
вах), что приводит к дополнительным временным задержкам и рас- 
ходам оперативной памяти. 


Алгоритм выбора структуры ИНС 
с помощью генетического алгоритма 
и грамматик графовой генерации Китано 


В общем случае подбор структуры ИНС прямого распростра- 
нения представляется сложной задачей. Существует ряд подходов к 
ее решению, из которых наиболее развитым является подход с 
использованием генетического алгоритма [19, 21] для поиска наи- 
лучшей структуры нейронной сети с точки зрения среднеквадра- 
тичной функции ошибки Ё, : 


А 1651 


У (апп(у,, Мека „) 


Е „апп, Г, ) = К (2.76) 


1е5ї 


№ 1е51 = 7. 


е5ї |, 


где: У. – репрезентативное множество тестовых векторов; а == 


1е51 1е5ї 


множество эталонных откликов для векторов из у 


1е51 * 
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Кодирование и декодирование структуры ИНС 


Важной проблемой, решаемой в процессе выбора структуры 
ИНС с помощью генетического алгоритма, является проблема ко- 
дирования структуры ИНС в виде непрерывной последовательнос- 
ти символов — последовательности хромосом, составляющих единст- 
венный ген генома каждой особи. 

Н. Киапо в [60] предложил способ кодирования структуры 
ИНС, называемый грамматиками графовой генерации Китано. 
Имеет смысл использовать одно из развитий грамматик графовой 
генерации Китано [21], которое заключается во введении допол- 
нительного параметра Ф, определяющего глубину разбора началь- 
ного правила грамматики и позволяющего задавать с помощью 
грамматик, содержащих одно и тоже количество правил, сети 
произвольного размера. 

Структура ИНС прямого распространения описывается ориен- 
тированным графом без циклов. Каждая дуга в данном графе соот- 
ветствует синапсу и имеет собственный вес, устанавливаемый в 
процессе обучения сети. Каждая вершина графа соответствует од- 
ному нейрону; с каждой вершиной связано пороговое значение 
соответствующего нейрона, настраиваемое в процессе обучения сети. 

Имеющиеся в графе дуги могут быть заданы с помощью 
матрицы смежности С (2.64) ИНС, в которой единица кодирует 
наличие ребра, направленного от вершины с номером, равным 
номеру строки, к вершине с номером, равным номеру столбца, а 
ноль — соответственно, отсутствие такового ребра. Матрицу С 
можно закодировать с помощью набора правил порождающей 
грамматики А: 


А={$,5,,5.,Е.,Е,,Е.,Е}; (2.77) 
$ = {0,11}; 5, =1{9,}; 1=1,16; 
$.={5.}; т=1, Н; Веб; 

1 0 


| 
Е =41 ‚ 0 ‚ № 0; 
1 1 0 0 
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а В 
Е. ={Е,„}; п} А; а, В,у,ӧ є(Е; Рр), 


где: 5 – набор терминальных символов; 5, – набор предтерми- 
нальных символов; 5. — набор нетерминальных символов; Ё; – 
набор правил вывода для терминальных символов; Ё, — набор пра- 


вил вывода для предтерминальных символов; Ё. — набор правил 


вывода для нетерминальных символов; А — стартовый нетерми- 
нальный символ. 


Грамматика А называется грамматикой графовой генерации 
Китано. Одна грамматика А соответствует одной особи популяции 
генетического алгоритма. 


Элементы 5, 5,, 5, Ё, Ё. фиксированы для любой А, и, 
следовательно, ген особи кодирует только набор правил Ё`.. Набор 
правил Ё. кодируется по схеме, приведенной на рис. 2.7. Каждое 
правило К), отображается в четыре хромосомы @„, В„, 7», ди, 
следующих друг за другом и за хромосомой д, –1 в единственном 
гене особи. 


нее 


Рис. 2.7. Схема кодировки набора правил 
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Алгоритм восстановления матрицы С из грамматики А состо- 
ит из следующих этапов: 
1) создается начальная матрица С = К; 


2) выполняется у =1 шаг применения правил Ё, ЕЁ и Ё: 


ар В 
С, =5—0(, =“ |; 
Ур Обь 
3) шаг 2 повторяется для каждого из правил @,, Вь, У», 02: 


< 3—1 
С, ә Ц. ; размерность С, оказывается, таким образом, равной 2909 


на 20) = 4 на 4 правила; 


4) шаг 3 повторяется (© - 2) раз – до получения матрицы С, ; 
5) в матрице С, любой символ 2#є 5} 5; рассматривается 


как символ у и преобразуется по соответствующему правилу в 0; 
6) матрица С. преобразуется в матрицу С по формуле (2.78) — 


таким образом, в матрице С гарантированно будут отсутствовать 
циклы: 


0, “> р, 
С=(8,,); С, =(2,,); 8, = (2.78) 
Е ТЕД 

Параметр 4 алгоритма восстановления матрицы смежности (© 
определяет максимально возможное количество нейронов, могущих 
быть в результирующей ИНС-2* '. нейронов. Учитывая тот факт, что 
число входов и выходов ИНС фиксировано, ток ИНС, описываемой 
соответствующим графом, присоединяются несколько нейронов: 

• входные нейроны в достаточном количестве связываются си- 
напсами с истоками графа; данные синапсы направлены от вход- 
ных нейронов к истокам графа; 

• выходные нейроны в достаточном количестве связываются 
синапсами со стоками графа; данные синапсы направлены от сто- 
ков графа к выходным нейронам. 

Мощность Н. множества < ,. оказывает, равно как и параметр 
Ф, существенное влияние на качество работы генетического алго- 
ритма, поскольку именно Н. определяет длину гена (как 4Н 
хромосом) и, соответственно, степень «разнообразия» правил, ко- 
дируемых геном. Выбор Н.. осуществляется эвристическим спо- 
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собом – небольшое значение Н. приведет к вырождению популя- 
ции, и, следовательно, схождению алгоритма к локальному опти- 
муму, тогда как большее значение Н ., наоборот, увеличит разно- 


образие особей, но также увеличит время схождения генетического 
алгоритма возможно, к глобальному оптимуму. 

По сравнению с традиционной версией грамматик графовой ге- 
нерации Китано, предложенная версия позволяет одинаково ком- 
пактно описывать нейронные сети с любым максимальным коли- 
чеством нейронов, что достигается за счет параметра Ч. 


Генетический алгоритм 


Генетический алгоритм, используемый для выбора структуры 
ИНС и использующий для кодирования структуры ИНС грамма- 
тики графовой генерации Китано, состоит из следующих шагов: 

1) случайным образом создается начальная популяция особей 


О = {8,}; 2=1Но; 


2) из множества У, случайным образом выделяются подмно- 


апп 


жества „и Г.,; для множеств Г, ИГ. из множества 4 


е5ї апп 
выделяются соответствующие им подмножества д, и а; 
3) для каждой особи #, восстанавливается соответствующая 
грамматика А, , производящая нейронную сеть аин, ; 
4) для каждой особи я. рассчитывается функция приспособ- 


ленности 5и (р, ): 


ѕий(е,) = Р ее В*; е=10"; (2.79) 


(атп,,/„„)+е 


апп 


5) особи е, сортируются в порядке убывания функции 
приспособленности 5ий(е.); 


6) особи с индексами 2є[1; Но,„| образуют множество особей 
Новы > ©соби с индексами 2е(Но»„,Новь»„ + Но» | образуют 


при этом должно выполняться 


множество особей Н ‚ особи с индексами (Нь. + Н, 


О сго55 ? О сго$5 ? 


образуют множество особей Н,, 


сго55 ? 


условие Нор, 22; 
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7) особи аЕО 


СГО55 


заменяются на особи 4», „ = С"055(9,,9,) ; 
4,9, ЕО,.„› то есть на особи – результаты скрещиваний случай- 


ных особей из множества (О, лучших особей популяции; 


8) особи Ч9ЕО„„ заменяются на особи О. = 


= тщашоп(а,); 9, ЕО,.., то есть на особи — результаты мутаций 

случайных особей из множества О,,, лучших особей популяции; 
9) если лучшая особь 4, популяции была лучшей и соответст- 

вовала условию останова Е 21, на менее чем Г 


апп 4 апп Ета апп веп тег 
предыдущих итераций алгоритма подряд, то выполняется переход 
на шаг 2, иначе алгоритм останавливается. Параметр Т есть 


апп Ета, 

уровень, определяющий достаточное качество обучения ИНС. 
По останову выполнения генетического алгоритма в качестве 
лучшей ИНС принимается ИНС, генерируемая грамматикой А 


(особь 4, окончательной популяции). Генетический алгоритм обла- 
дает следующими параметрами: 

У |Г..| = определяют размеры обучающей и тестовой 
выборок векторов; значения данных параметров зависят от 
предметной области и сложности задачи классификации, решаемой 
искомой ИНС; 

"Но, Но 
няемых при переходе к следующей итерации генетического алго- 
ритма, заменяемых на результаты скрещивания и мутации; 

• сто55, тшашоп — функции скрещивания и мутации соот- 


еасй е5! | 


Н 


О сго55 


— определяют количество особей, сохра- 


Без ? 


ветственно; 

°Т„к. — Определяет верхний уровень качества обучения луч- 
шей ИНС, достаточный для останова алгоритма; 

°1Т тень — Определяет достаточное для останова алгоритма 


количество итераций, на которых подряд должны повторяться 
условия останова для лучшей особи; очевидно, что чем больше 
‚ тем больше вероятность схождения алгоритма к 


апп реп тег 
устойчивому оптимуму – к ИНС, дающей одинаковый по качеству 
обучения результат на произвольных У, и, 


1е51 ° 


88 Методы совмещения разнородных изображений путем поиска 
и сопоставления образов объектов 


2.4. Поиск соответствий между образами, 
выделенными на разнородных изображениях 


Алгоритм формирования множества Р 


’'еи Пар совпадающих 
точек по множеству классифицированных образов С , выделенных 
на спутниковом снимке ДВ, состоит из следующих этапов: 
1. Для каждого образа ѓєс; сє С выполняются следующие 
действия. 
1.1. С помощью функции соогА рассчитывается множество 


В, є М точек, таких что рє Р, <> р=соога(,); 6, ЕЁ. 


1.2. Из множества Р, точек выделяется подмножество 


Р, Іеғреѕї с В ‚ такое что Уре В, Іегре5і относятся к одному И 
тому же объекту на карте и выполняется условие: 
УР, 5иф5е! = В, 
Р, Іеғееѕї < В, 5иб5е! а о; | В, Іеғреѕт р Р, ѕиЬѕе! | у (2.80) 


1.3. Если объект ЦКМ, описываемый множеством Р, 


Іеғреѕі ? 
и образ Е относятся к одному и тому же классу, то гео- 
метрический центр (х, у) объекта, содержащего Р, И 


1егвез! ? 


центроид (1, /) образа # включаются в множество Риме = 
= (00, /), (ә, у). 
2. Для множества Р,„‹г’ рассчитывается совместная гисто- 
грамма отклонений дх, иду: 
(х', у’) = сооға((і, Ј)); 
0х=|х'-х|; (2.81) 
бу=|у-У|. (2.82) 
3. Так как искомое геометрическое преобразование единствен- 


но, то гистограмма будет иметь одну ярко выраженную моду (а 
совместное распределение дх, и бу будет нормальным) – пары 


точек ((/, /),(х, у), формирующие моду, включаются в множество 
Р 


ғау еди 


пар совпадающих точек с использованием правила 3-х сигм. 
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МЕТОДЫ КОНТУРНОГО КОРРЕЛЯЦИОННОГО 
СОВМЕЩЕНИЯ РАЗНОРОДНЫХ ИЗОБРАЖЕНИЙ 


3.1. Алгоритм моделирования 
системы искусственного интеллекта 


На вход системы искусственного интеллекта, выполняющей 
совмещение данных ДЗЗ и ЦКМ с помощью контурного корреля- 
ционного совмещения, основанного на использовании нейросете- 
вой ассоциативной памяти, поступают следующие данные: 

В – исходный спутниковый снимок; 

М - исходная ЦКМ; 

соот – функция преобразования координат (х, у ) пикселя В, 


спутникового снимка В в координаты (х,, у, ) точки на ЦКМ М в 


соответствии с системой координат и проекцией спутникового 
снимка В. 

На выход системы искусственного интеллекта поступает 
множество пар совпадающих точек Р (1.10), по которым в 


гам’ ейи 
дальнейшем выполняется расчет множества пар точек Р, (1.11) и 
функции сотр" (1.12), описывающей результирующее совмещение 


спутникового снимка В и ЦКМ М. 

Система искусственного интеллекта состоит из следующих 
модулей: 

1) модуль предварительной обработки спутникового снимка с 
целью устранения несущественных (слабых или прерывистых) 
контуров (3.1); 

2) модуль выделения контуров на спутниковом снимке (3.2); 

3) модуль описания контуров (3.3); 

4) модуль поиска совпадений контуров образов на спутнико- 
вом снимке с границами объектов на ЦКМ (3.4). 

Р = лпа.„(В,М,соота)= 


гам’ еди 
соггёаной(М ‚саезс(сопюиг( рғерағе(В)))); 


В, = ргераге(В); (3.1) 
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Ви = сотоиг(В„); (3.2) 
С = саезс(В„); (3.3) 
Р „„ =соттеаной(М,С); (3.4) 


Виа}; 1=ЬГ = Б, [0]; С= {2}. 


Исходный спутниковый снимок В представляет собой вход- 
ной параметр модуля предварительной обработки ргераге , резуль- 


татом работы которого является изображение В’, на котором, в 


р'р ? 
отличие от спутникового снимка В, отсутствуют несущественные 
(слабые или прерывистые) контуры. 


Изображение В„, подается на вход модуля сотоиг , который 


выполняет выделение контуров на изображении. На выход модуля 
сопюиг поступает контурная матрица В, каждый пиксель В 


спі ? ст іј 


которой содержит меру силы контура, проходящего через пиксель 
№. у исходного изображения В„,. Под мерой силы контура в (3.2) 


ст іў 
понимается численная мера, принимающая значения из диапазона 
[0;1], где 0 соответствует отсутствию контура, а | – максимально 
сильному контуру. В простейшем случае В, является бинарной 
матрицей. 

Контурная матрица В, представляет собой входной параметр 
модуля описания контуров сес , результатом работы которого явля- 
ется множество С описаний © контуров, выделенных на изобра- 


жении В... 


Множество С вместе с ЦКМ М поступает на вход модуля 
соттёаноп ‚ выполняющего корреляционное совмещение ЦКМ и 
спутникового снимка — контурам, описание которых присутствует 
в С , ищутся соответствия на ЦКМ (границы объектов), по 
результатам чего формируется множество пар совпадающих точек 
Р по которым в дальнейшем выполняется расчет множества 


гам’ ейи ? 
пар точек Р, (1.11) и функции сотр' (1.12). 


В $ 1.3.2 приведен детальный анализ возможных подходов к 
реализации каждого из модулей системы. Окончательно были 
выбраны следующие подходы: 
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Модуль рғерағе – ИНС, выполняющая анализ независимых 


компонентов ((англ.) шдереп4епе сотропепі апа]уѕіѕ, [СА), и алго- 
ритм, основанный на фильтре низких частот и фильтрации по 
Тихонову; 

Модуль сопюиг – нелинейная ячеистая нейронная сеть, ней- 
ронная сеть, основанная на применении вейвлет-преобразования 
Габора; 

Модуль саезс — метод, основанный на составлении множества 
подматриц матрицы саезс фиксированного размера (множества 
окон), содержащих контур; 

Модуль соггёаНои – ограниченная стохастическая машина 
Больцмана. 

Схема алгоритма моделирования системы искусственного ин- 
теллекта приведена на рис. 3.1. Алгоритм моделирования состоит 
из следующих этапов: 


Описание контуров 


Получение 


настроек 
системы 


Шоиск совпадающих 


Предварительная 
обработка спутни- 
кового снимка В 


Р, гош еди 


Выделение кон- 


туров на спут- 
никовом снимке 


Рис. 3.1. Алгоритм моделирования системы 
искусственного интеллекта 
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1) предварительная обработка спутникового снимка с целью уст- 
ранения несущественных (слабых или прерывистых) контуров (3.1); 

2) выделение контуров на спутниковом снимке (3.2); 

3) описание контуров (3.3); 

4) поиск совпадений контуров образов на спутниковом снимке 
с границами объектов на ЦКМ (3.4). 


3.2. Предварительная обработка спутникового снимка 
3.2.1. Описание алгоритма 


Исходный спутниковый снимок В вследствие различных при- 
чин, к которым относятся разнообразные шумы, недостаточное 
разрешение аппаратуры и прочие, может содержать большое коли- 
чество мелких и прерывистых контуров, непригодных для исполь- 
зования в процессе дальнейшего совмещения. 

Наличие таковых контуров требует выполнения дополнитель- 
ной предобработки спутникового снимка и постобработки (фильт- 
рации) выделенных контуров с целью устранения контуров, непри- 
годных для дальнейшего совмещения. 

Модуль предварительной обработки спутникового снимка 
ргераге минимизирует влияние на результат совмещения несу- 
щественных (слабых или прерывистых) контуров путем выполне- 
ния следующих действий: 

1) фильтрация іса аддитивных помех, имеющихся на спутни- 
ковом снимке В (3.5), для чего используется ИНС, выполняющая 
анализ независимых компонентов ((англ.) шдереп4епе сотропепї 
апа1уз1$, ІСА); 

2) фильтрация ИЙег несущественных (слабых или прерывис- 


тых) контуров (3.6), для чего используется алгоритм, основанный 
на применении фильтра низких частот и фильтрации по Тихонову. 


В, = ргераге( В) = ПЙег(са(В)); 
В’. = іса(В); (3.5) 
В„„ = ЛИег(В,,). (3.6) 
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3.2.2. Алгоритм устранения слабых 
и прерывистых контуров 


Описание алгоритма 


Алгоритм фильтрации на спутниковом снимке слабых и пре- 
рывистых контуров [61] имеет вид: 


В „„ = втеу( ЛИег(аеМиттив ([оухВ,,))); (3.7) 
В, „ = Іои(В, ); (3.8) 

В ы = ЯеБигғіпе (В, ); (3.9) 

В„„ = втеу(Виы ) (3.10) 


и состоит, таким образом, из следующих этапов: 

1) Јоу (3.8) – поканальное применение к изображению В,, 
фильтра низких частот Јоу, следствием чего является сильное 
размытие В,., – изображение В, 


ИГ 

После применения к изображению фильтра низких частот 
(ФНЧ) все линии сглаживаются, но основная информация об 
объектах на изображении не теряется. Вся информация перераспре- 
деляется по некоторому закону и может быть восстановлена, за 
исключением мелких объектов и тонких линий, которые несущест- 
венны для дальнейшего совмещения; 

2) аеМиттия (3.9) – поканальное восстановление четких гра- 
ниц размытых объектов ((англ.) де иггип?). 

После того, как с помощью ФНЧ тонкие линии и мелкие объек- 
ты были практически устранены из изображения В, , выполняется 


ом ? 
восстановление четких границ крупных объектов с получением на 
выходе изображения В, ; 


3) етеу (3.10) – преобразование изображения В, в однока- 
нальное изображение В, для упрощения дальнейшего выделения 
контуров, значимых для решения задачи совмещения изображений. 


Фильтр низких частот 


В [61] было выполнено сравнение ФНЧ, которые могут быть 
использованы в составе описываемого алгоритма, по результатам 
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которого наилучшим выбором для реализации функции /оу был 
признан медианный фильтр. 

Медианный фильтр является пространственным фильтром, не 
описываемым ядром свертки. Значения отсчетов внутри окна 
фильтра сортируются в порядке возрастания (убывания), и значе- 
ние, находящееся в середине упорядоченного списка, поступает на 
выход фильтра. В случае четного числа отсчетов в окне выходное 
значение фильтра равно среднему значению двух отсчетов в сере- 
дине упорядоченного списка. Окно перемещается вдоль фильтруе- 
мого сигнала, и вычисления повторяются. 

Центральный элемент заменяется медианой всех элементов 
изображения в окне. Медианой дискретной последовательности 


№М+1 
{и,...и„..и. } для нечетного № является элемент ДЛЯ 


которого существуют элементов и меньших или равных ему 


по значению, и 


больших или равных ему по значению. 


Восстановление четких границ размытых объектов 


По результатам выполненных в [61] экспериментов в качестве 
алгоритма восстановления границ был выбран алгоритм, основан- 
ный на сглаживающей фильтрации методом наименьших квадратов 
со связью (фильтрация по Тихонову) [9]. 

Алгоритм восстановления границ выполняет поканальное 
восстановление искаженного изображения В, с помощью модели 


ом 
процесса искажения дејоғт . 
В качестве модели процесса искажения аеютт в [61] был 


выбран ФНЧ Гаусса: 


а? +? 
аејоғт(а,р) = Е и 4 ) (3.11) 
2ло 
где параметр с влияет на остроту пика производящей функции 
Гаусса. 


По модели процесса искажения дејоғт создается искажаю- 


щий пространственный фильтр К. „„ · Сам процесс искажения 


Глава З 95 


имеет вид (3.12) в пространственной области или как (3.13) в 
частотной: 


В. = Кае рт" В, (3.12) 
Ра 5 г. Јогт ; В;; (3 : 13) 


В; | ыы ай (В, Е ); 
Ку Јоғт С ай (К де јот ); 
В; = а1(В,), 
где Ай – прямое дискретное преобразование Фурье. 


Для восстановления изображения необходимо решить задачу 
нахождения экстремума (минимума) некоторого сглаживающего 
функционала. В методе, основанном на фильтрации по Тихонову, в 
качестве такого функционала Р можно использовать квадрат нор- 
мы лапласиана: 


Р = у Ув, ў (3.14) 


т=0 п=0 
с дополнительным ограничением (связью) вида: 
– К 


а ОДЕ 
АА =р?. (3.15) 
пєК. 
Так как алгоритм фильтрации контуров предполагает приме- 


нение ФНЧ к исходному спутниковому снимку (т.е. известна точ- 
ная модель процесса искажения), коэффициент 7 в (3.15) прини- 


мает значение 7 =0 и выражение (3.15) сводится к: 
АА = 0. (3.16) 
Решение оптимизационной задачи (3.14) при условии (3.16) в 
частотной области имеет вид: 


* 


Ка 
* = Јогт * ; 3 17) 
ай к * 2 * |2 [оу К > ( ; 
| К 2 Јотт | ия | Г 
7ЕК; 


Г = 41), 
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оо 
І=|-1 4 1, (3.18) 
оо 


где: коэффициент у должен быть выбран таким образом, чтобы 
выполнялось условие (3.16); матрица Г – оператор Лапласа (3.18). 
< * ұғ 

После того, как найдено решение В, , оптимизационной 


задачи (3.17), рассчитывается В, ,: 


Вы, = (Ву 1), (3.19) 
где функция і27 есть обратное преобразование Фурье. 


Преобразование многоканального изображения 
к одноканальному 


Изображение В, необходимо преобразовать к одноканально- 
му изображению В,, с целью улучшения временных характерис- 
тик процесса выделения контуров и повышения качества его 


работы. 
Преобразование многоканального изображения В, к однока- 


нальному В, имеет вид: 


В„„ = втау( Вы }; (3.20) 


К 
уз" ИК 
= — к=] 
В рӯ Б 2ғау(В у) Га 
3.3. Выделение контуров на спутниковом снимке 


3.3.1. Описание алгоритма 


На вход модуля выделения контуров соиюиг поступает изоб- 


ражение В,,, , полученное с выхода модуля предобработки спутни- 


кового снимка В. Результатом работы модуля выделения контуров 
является бинарная одноканальная контурная матрица В, мар- 


спі ? 


кирующая наличие или отсутствие в каждом пикселе спутникового 
снимка контура, пригодного для использования в ходе дальней- 
шего совмещения В и ЦКМ М. 
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Под «пригодностью» в данном случае понимаются следующие 
требования к выделенным контурам: 

• достаточная протяженность контура — чем длиннее контур, 
тем больше информации он предоставляет для дальнейшего совме- 
щения спутникового снимка и ЦКМ; 

° толщина контура должна быть равномерной на всем его 
протяжении. 

Для выделения контуров на спутниковом снимке В, удовлет- 
воряющих перечисленным требованиям, применяются следующие 
методы: 

1) выделение контуров с помощью нейронной сети, исполь- 
зующей вейвлет-преобразование Габора; 

2) выделение контуров с использованием нелинейной ячеистой 
ИНС. 


3.3.2. Алгоритмы выделения контуров 
на спутниковом снимке 


Алгоритм выделения контуров 

на спутниковом снимке 

с использованием нелинейной 

ячеистой искусственной нейронной сети 


Нелинейная ячеистая ИНС [(англ.) поп Ппеаг сешаг пеџга] 
пебмогк, МЕСММ] [43—46] является двухмерной искусственной 
нейронной сетью размером / на .Л нейронов (т.е. размер МЕСММ 


совпадает с размером изображения В,,, ). 


Каждый (1, /)-й нейрон сети рассчитывает свой выход У, В 


момент времени / по набору входных данных с учетом обратной 
связи, имеющейся внутри нейрона. Существуют различные подхо- 
ды к организации обратной связи для выделения контуров. В нас- 
тоящей работе используется сеть М.СММ, работа которой дискре- 
тизирована по времени — обратная связь оказывает влияние на рабо- 
ту нейрона с задержкой в одну итерацию работы сети. 


Выходы У; нейронов сети образуют матрицу У’ выходов 


размером / строк на Ј столбцов. 


На рис. 3.2 приведена схема структуры нейронной сети 
МЕСММ. 
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7 


тіспт 17 


2. а 


Рис. 3.2. Структура (1 /)-го нейрона нелинейной ячеистой ИНС 


Алгоритм работы нейрона (і, /) на итерации { работы сети 
состоит [43] из следующих этапов. 


< 1 
1. На вход нейрона поступает подматрица В, матрицы В», 


поступают выходы ү нейронов сети на предыдущей итерации (на 


–] 


о, | у] ў ч, 
первой итерации —=0; “й, /:Ү; =0) и вес х; обратной связи, 


при этом В, <В, и т с Ү'', то есть каждый нейрон обладает 


своим уникальным двумерным локальным сенсорным полем т 


разме-ром М’ на № элементов. 
В качестве В, предлагается использовать окно размером ( М", 


№) пикселей на изображении В > с центром в пикселе с коорди- 
натами (1, /). 


2. Для каждого нейрона ({, ј) повторяются следующие 
действия. 


2.1. Рассчитываются Е И В. 


ОЗ С (3.21) 
(ГЛЕР, 

2. А (7), (3.22 
(й, Ј)єР; 


1 а уг 
О = И 


1) ? 
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1 Г > 
где функции Г, и Р; – нелинейные (в общем случае) функции 


< і 1—] 
расчета характера связи нейрона с элементами матриц В, И у, 
2.2. Рассчитывается и 


х= хк, + Рі +Т! 


іј піспп |? 


(3.23) 
где Т' 


піспп іј 


— пороговое значение. 


2.3. Рассчитывается {, : 


б = ит (5), (3.24) 


где /„„ - функция активации нейрона. 


3. Шаг 2 повторяется до тех пор, пока не выполнится одно из 
следующих условий: 

Сеть выполнила достаточно большое количество итераций. 

Выполняется условие: 


РДЭ] = ү" Ё |< Г піспп ? (3.25) 


где пороговое значение Г, >20 и выбирается в зависимости от 


піспп 


характера задачи. 
4. Формируется матрица В,, 


В. = ‘бы е (3.26) 
0, уу =0, 
К І, у; #0. 


Важно отметить, что МЕСММ является самообучающейся ИНС, 
сходящейся к оптимальному с точки зрения поставленной задачи 
результату. 

Кроме того, необходимо отметить также тот факт, что процесс 
работы МЕСММ хорошо поддается распараллеливанию, что обусл- 
овлено ее модульной структурой. 

В [43, 45, 46] обосновывается использование для выделения 
контуров на спутниковом снимке следующих значения параметров 
работы сети: 


М = М =3; 
оер 


піспп іу 
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[585 РОЛ 
Т ОР е. 
86,» (ПЕРО); 
0, (= 7’); 

е В; = 


1—1 +! | 
2у,, (= Ј?). 
В [45, 46] приведены результаты сравнения различных функ- 
ций /,,„„ активации нейронов сети для решения задач выделения 


контуров, исходя из которых в качестве функции /,„ следует 
выбрать гиперболический тангенс (3.27): 


Рот) = Р(х). (3.27) 


Нейросетевой алгоритм выделения контуров 
на спутниковом снимке, основанный 
на использовании вейвлета Габора 


Анализ информации, поступающей в мозг человека от первич- 
ных зрительных органов, происходит с использованием крайне 
сложной системы нейронов. В частности, как показано в [26], для 
получения информации о контурах на изображении мозг человека 
использует сложные преобразования, схожие с вейвлет-преобразо- 
ванием Габора. Поэтому вейвлет Габора используется для достаточ- 
но точного моделирования действий зрительной системы человека 
по выделению и анализу контуров объектов на изображении. 


Двухмерное вейвлет-преобразование Габора имеет вид: 
12 '2 


ваћон(х', у) = ехр| ^^? — |соѕ(Ах'+ оу”; (3.28) 
2ли 


х' = хс0$© + узш®; у’=-хз Ө + усоѕ6. 
На основе вейвлета Габора строится пространственный фильтр 
5 аво (№5. @,©), размер которого 5„»„ и„ адаптивно настраивает- 


а 
ся в зависимости от масштаба изображения. 


Материнская функция вейвлета Габора (3.28) обладает следую- 
щими параметрами: 


° и – стандартное отклонение распределения Гаусса, входя- 
щего в состав вейвлета Габора; 
• А4, о – коэффициенты при координатах хи у; 
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• © – угол относительно оси ОХ, относительно которого оце- 
нивается сила контуров. 


Для выделения контуров на изображении необходимо сгене- 
рировать набор фильтров Ё: 


Е ЕТЕ ,.. Јон (А0, (т = ПОВ 7 ариз (3.29) 


М = 5 


• 
вађБог тит ? 


м>0; 89-2. (3.30) 
М 


Применение набора фильтров Е к изображению В, поз- 


спі ғау (3. Э | ), 


каждое из которых содержит оценку силы контуров, идущих под 
углом (т – 1)д0 , в пикселях исходного изображения В,,, 


В ны 


ст газу ={В спі гам’ 1? * спі ғам т? °**%°* сш гам М } у 


воляет получить на выходе набор изображений В 


(3.31) 


Применение набора фильтров Ё к изображению В„, описыва- 


ется как обработка изображения В„ с помощью многослойной 


ИНС Ж. Действительно, каждый из фильтров можно рассматри- 


вать как отдельный нейрон с линейной функцией активации. 
Задача выделения контуров на изображении заключается в вы- 


боре контурных пикселей В„,„, исходного изображения по резуль- 


татам анализа набора изображений В 


спі гаи’ * 


Для принятия решений о нахождении в пикселе 5. контура 


ртр іј 
используется многослойная ИНС Ж, , структурная схема которой 
приведена на рис. 3.3. 


Алгоритм работы ИНС Я), для пикселя Б имеет следующий 


ртр іј 
ВИД. 
1. Для каждого т- го нейрона первого слоя для т=2,М-1| 
выполняются следующие действия. 
1.1. На вход нейрона подается тройка пикселей /, 


Г. = у 12 а Е, (3.32) 
| 


т-1 сш тамт-1/ ? 


Г =6 


т спі га т іу ? 
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т+1 — Эу ғам т+1 1) * 


Второй слой нейронов 


Инн 
о 


о т 


асіімасіоп_ Беа Ь_сп(т, п) 


Нейрон (т + 1) 


Рис. 3.3. Структурная схема многослойной ИНС, 
используемой для принятия решения о нахождении контура в пикселе 


Ь (мт + 1) 


Ь (т + 2) 


| 
1.2. Тройки пикселей суммируются [„ „= х Ре 


п=-] 


1.3. К сумме /,, „ применяется пороговая функция акти- 
вации /„; с порогом Ти: 
=: 0? ) = т РЕ т ЕТ 
т Јас 1 т ѕит2 : 
0; Е (ЕРЕ (3.33) 
Бра Є [0; 1]. 


2. Единственный нейрон второго слоя суммирует свои входы: 
=), у, 
т 


3. К сумме /, применяется пороговая функция активации #.,,, 


с порогами Ти, И Тиз: 

5 І, & [Т афог 25 ПАРЕ |. 
Ву = Лиз = М (3.34) 
0; 2 «Та ^1 >Т 


вађБоғ3* 


+ 
Во 2» 1 Борз Є Х Е 
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Первый слой ИНС Ж, оценивает силу контуров на изоб- 


ражении В„, в соответствующих диапазонах направлений. Второй 


слой ИНС И, оценивает наличие контуров в каждом из пикселей 


изображения — небольшое значение суммы [., предполагает нали- 
чие в пикселе слабого контура, относящегося, скорее всего, к 
шумовым контурам; слишком большое значение суммы [, ука- 


зывает на нахождение пикселя на пересечении нескольких кон- 
туров — пиксель может находится на пересечении двух или трех 
контуров, но, когда количество таковых контуров больше трех, 
данный пиксель можно удалить скорее всего, каждый из контуров, 
проходящих через пиксель, является достаточно широким и анали- 
зируемый пиксель попадает в центр линии контура. 

Окончательно алгоритм выделения контуров на изображении 
состоит из следующих основных этапов: 


1) прогон изображения В„, через ИНС Ж, с получением на 


выходе набора изображений Вии: Ви =Й' (Вы, 1, А, 0,5 ово тит? 
О тег) , 
2) прогон набора изображений В.„„„ через ИНС Ж, с получе- 


нием на выходе бинарного контурного изображения В, : В, = 
=, (ВТ 


спі гам’ ? > вафог | > Т вафог 2? Т ғађоғ 3 ). 
Предложенный алгоритм выделения контуров обладает сле- 
дующими параметрами: 


® ИНС И И, А, @2 5 льо утэ 
ИНОЙ РЕ, 


вађог1? > рађБог 2?“ рађоғ3 * 


К) 


Лек ? 


Параметры используемых ИНС подбираются эксперименталь- 
но с целью ускорить процесс подбора. 

Описанный алгоритм выделения контуров хорошо поддается рас- 
параллеливанию. Каждый нейрон ИНС И/ и И, является обособлен- 


ным модулем, работающим независимо от прочих нейронов. Парал- 
лельная работа нейронов используемых ИНС была орга-низована на 
базе видеокарты, поддерживающей технологию СРСРО, поскольку 
видеокарты данного класса располагают большим количеством неза- 
висимых вычислительных ядер в отличие от много-ядерных процес- 
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соров общего назначения, располагающих в лучшем случае 16-ю и 
32-мя вычислительными ядрами. 


Выбор алгоритма выделения контуров 


Каждый из рассмотренных алгоритмов выделения контуров ис- 
пользуется в составе описываемого алгоритма совмещения спутни- 
ковых снимков и ЦКМ. 

Нейросетевой алгоритм выделения контуров на спутниковом 
снимке, основанный на использовании вейвлета Габора, за счет 
более низких временных показателей используется для выполнения 
оперативного совмещения спутниковых снимков и ЦКМ. В тех же 
случаях, когда временные характеристики процесса совмещения 
оказываются менее важны, чем качественные, для выделения кон- 
туров на спутниковом снимке должен использоваться алгоритм, 
основанный на применении нелинейной ячеистой искусственной 
нейронной сети для выделения контуров. 


3.4. Описание контуров 


По результатам анализа способов описания контуров, выпол- 
ненного в $ 1.3.2, в качестве средства описания контуров был выб- 
ран способ, основанный на составлении множества подматриц мат- 
рицы В., фиксированного размера (множества окон), содержащих 
контур. 

Алгоритм формирования множества С, ={#„} описаний №, 


контуров, присутствующих на спутниковом снимке В, состоит из 
следующих этапов (здесь и далее символом г обозначается кон- 


тур, окно, в котором находится контур (контурное окно), и вектор 
описания контура, получаемый разверткой контурного окна пост- 
рочно): 


1) матрица В,, разделяется на множество Св = {2„ } непере- 


ст 


секающихся квадратных окон размером 5 на 5 


зла ѕіге па ѕіге 


элементов; 
2) окно 8, Е Сь <> Тс ЗУ вы 575 тах» ГД сы И сыа 


минимальное и максимальное количество элементов окна @ ви › 
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занятых контуром, достаточные для использования контурного ок- 
на в процессе дальнейшего совмещения. 

Пороговое значение Г.„„ Позволяет исключать из рассмот- 
рения контурные окна, содержащие небольшие (несущественные с 
точки зрения дальнейшего совмещения) контуры и фрагменты 
контуров. 


Пороговое значение у: позволяет исключать из рассмот- 


тах 


рения контурные окна, содержащие избыточную информацию о 
контурах, приводящую к усложнению (и зачастую к невозмож- 
ности) сопоставления контуров, выделенных на спутниковом 
снимке, с контурами объектов на ЦКМ. 

Протяженные контуры, выделенные на спутниковом снимке 
В, для которых сторона минимального прямоугольника, содержа- 


щего контур, превосходит значение 9 пикселей, оказываются 


упа ѕіге 
искусственно разделенными на несколько контурных окон, что 
положительно сказывается на качестве выполняемого совмещения, 
так как позволяет получить больше достоверного контурного 
материала для выполнения совмещения, поскольку чем протяжен- 
ней контур, тем больше вероятность того, что соответствующий 
ему объект (или группа объектов) присутствует на ЦКМ и форми- 
рует соответствующий контур на ЦКМ. 

Множество С, = {е,„} описаний контуров объектов, при- 


сутствующих на ЦКМ М , формируется аналогичным образом. 
3.5. Поиск совпадений контуров образов 


3.5.1. Структура коррелятора, 
принципы его обучения и функционирования 


Существует несколько подходов к решению задачи контурного 
корреляционного совмещения (задачи поиска соответствий между 
контурами, выделенными на спутниковом снимке В, и контурами 
объектов, присутствующих на ЦКМ М ). Как было показано в 
$ 1.3.2, наиболее эффективным с точки зрения времени выполне- 
ния и требуемых вычислительных ресурсов и позволяющим полу- 
чить при этом результат качества, достаточного для дальнейшего 
совмещения, является подход, базирующийся на использовании 
нейросетевой ассоциативной памяти, основанной на ограниченной 
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стохастической машине Больцмана ((англ.) геѕігісіеа Воіїгтапп 
тасһіпе, КВМ) [46, 47]. 

Пусть С, ={#„} – множество описаний контуров &,„„, выде- 
ленных на спутниковом снимке В(| С, |= №,); С, = {2,,„} – мно- 


жество описаний контуров &,,„ объектов, присутствующих на 
ЦКМ (| Си |= А 


Размерности векторов е, и р, очевидно совпадают и равны 
2 


| & вн |= 8 Ми |= ла ѕіе * 
Алгоритм работы коррелятора состоит из следующих этапов: 
1) обучение коррелятора на множестве С, описаний контуров, 


выделенных на спутниковом снимке В, с использованием алго- 
ритма СО-К ((англ.) сопігаѕіуе діуегрепсе), описанного в [62]; 


2) прогон коррелятора для каждого е, є С, и формирование 


пар описаний контуров (2,„,2,,)ЄС,„, таких, что для контура 


Р 


14 


м, Контур #,, является наилучшим соответствием с точки зре- 
ния коррелятора; 
3) формирование множества пар совпадающих точек Р по 


газ еди 


множеству С 


ра * 

Ключевым элементом представленного алгоритма является 
коррелятор, основанный на применении ограниченной стохасти- 
ческой машины Больцмана. 

Ограниченная стохастическая машина Больцмана относится к 
классу ИНС, реализующих ассоциативную память. Обучение сто- 
хастической машины Больцмана сводится к запоминанию входного 
множества векторов, после чего на вход машины подается неко- 
торый входной вектор и машина выполняет ряд итераций до тех пор, 
пока выходы нейронов ее входного слоя не достигнут стабильного 
состояния, которое и будет являться вектором из обучающего мно- 
жества, восстановленного для данного входного вектора. 

Структурная схема ограниченной стохастической машины Больц- 
мана приведена на рис. 3.4. 

ВВМ состоит из двух слоев нейронов: 
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• слой видимых нейронов, векторы весов которых обозначают- 
ся как И где 4 – индекс нейрона. 


уі5іЫе д ? 
| Й чые 4 


т. 
• слой скрытых нейронов, векторы весов которых обозна- 


чаются как Иа. 


Векторы весов нейронов каждого слоя образуют матрицы 


14 и Й агп > Количества строк которых Ж, „„ь И Эль маен 


уе 
являются количествами нейронов В ВИДИМОМ И скрытом СЛОЯХ 
соответственно. 


_ |1 Видимый слой 


Скрытый слой 


Рис. 3.4. Структурная схема 
ограниченной стохастической машины Больцмана 


Так как ограниченная стохастическая машина Больцмана явля- 
ется полносвязаннои ИНС, то ыга Е И т зе И 5 №адеп — 


= иіѕіЫе 9 Е Эн зе - 
Алгоритм работы ИНС по сопоставлению с контуром #3 
контура в, є С, состоит из следующих этапов: 
1) исходный вектор ғ, поступает на вход ВВМ: е„„, = & м; 
2) иные ОИЕ 


3) рассчитывается вектор 8» = {Вы а} › ГДЕ Файт а = 


2 Т : 22 $ 
— И ые 48 %5іЫе 4 1, От уе ? 


108 Методы контурного корреляционного совмещения 
разнородных изображений 


4) рассчитывается вектор &ы, = {8,4} › ГДЕ вые а = 


= Т . и ь 
Бр Й ален 48 таап? Ч = 1, 5 гЬт А4еп ? 


5) рассчитывается ӧу,,,„,: 


т ме 
к )? 
(2 уе 4 8 уіѕіБе ргеу 4 
= 9=1 ; 
Об ые 25 ) (3 35) 
гфт уіѕіЫе 
6) если бе, < Т, зор › ТО выполняется останов алгоритма, 


иначе выполняется переход на шаг 2. 
По завершении выполнения алгоритма очередной вектор & ны. 


принимается в качестве приближения е,’ к вектору #„„. 
+ 
Пороговое значение Т,„„„„, Е К является максимально допус- 


тимым отклонением очередного вектора #„,, от соответст- 


вующего вектора на предыдущей итерации, при котором выпол- 
нение алгоритма останавливается и ИНС считается достигшей 
стабильного состояния. 


—^ 


Поиск вектора #,„’ по его примерному значению р,’ из мно- 


жества С, выполняется по алгоритму, состоящему из следующих 
этапов. 


1) для каждого 2, є С, рассчитывается др, : 


У 5іге __ в 2 
а=1 Каа 5 в„'а ) 
9 


упа уе 


бе» = (3.36) 


2) выбирается вектор ғ, такой, что и'= агртіпдұе,. 


Таким образом, в качестве р, ' выбирается вектор ғ, 'єС,, 


максимально близкий к ©,’ в смысле Евклидовой меры 
расстояния. 
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3.5.2. Формирование множества пар совпадающих точек 


Алгоритм формирования множества Р (1.10) пар совпа- 


гам’ ейи 


дающих точек по множеству С пар совпадающих контуров сос- 


рат 
тоит из следующих этапов: 


1. Для каждой пары контуров (2,„,2„)еЕС, выполняются 


раі" 
следующие действия: 
1.1. Для контурного окна (контура) в,’ , выделенного на 


контурной матрице В соответствующей спутниковому 


сш ? 


снимку В, рассчитываются координаты (1., /,.) его цент- 


рального пикселя 5, , ЁБ. . 


1.2. Для контурного окна (контура) 2, , выделенного на 
ЦКМ М, рассчитываются координаты (х,, у, ) его центра. 
1.3. Рассчитываются отклонения дх и ду: 
(х,у) = соота((ї,,),.)); 
хх. -х, | (3.37) 
бу=у, - У, |. (3.38) 
2. Рассчитывается совместная гистограмма отклонений дхи ду. 


3. Так как искомое геометрическое преобразование единст- 
венно, то гистограмма будет иметь одну ярко выраженную моду (а 
совместное распределение ӧх и ду будет нормальным) — пары 


точек ((1., /.), (х,, у, )), формирующие моду, включаются в 


множество Р. пар совпадающих точек с использованием 


“ ейи 


правила 3-х сигм. 
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ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ МЕТОДОВ 
СОВМЕЩЕНИЯ РАЗНОРОДНЫХ ИЗОБРАЖЕНИЙ 


4.1. Описание программного стенда 


Алгоритмы моделирования систем искусственного интеллекта, 
численные методы настройки их параметров, численные методы 
обучения классификаторов, а также непосредственно системы ис- 
кусственного интеллекта, выполняющие совмещение данных ДЗЗ и 
ЦКМ, были реализованы в составе следующих программных про- 
дуктов: 

• программа совмещения данных ДЗЗ с ЦКМ с использованием 
искусственных нейронных сетей — программная библиотека и поль- 
зовательский интерфейс к ней; 

• программный стенд тестирования алгоритмов предобработки 
и совмещения изображений; 

• программа построения и уточнения карт гарей лесных масси- 
вов на основе данных ДЗЗ с использованием системы искусст- 
венного интеллекта. 

Для каждого из перечисленных программных продуктов име- 
ется свидетельство о его регистрации в ФГБУ «Федеральный инс- 
титут промышленной собственности Федеральной службы по ин- 
теллектуальной собственности, патентам и товарным знакам» (ФГБУ 
«ФИПС»-РОСПАТЕНТ). 

Разработка программных продуктов была выполнена с помо- 
щью следующих инструментов: 

• языки программирования: С++ (стандарт С++11), Руфоп 3, 
Гла, В; 

• программные библиотеки: С16С, ОрепС\У, О 4, Ої 5, $01; 

• спецификации: ОрепОТ,, ОрепМР, ОрепСГ; 

• программные утилиты: Дохуреп, СМаке, Сі, Мт. 
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4.2. Порядок проведения 
экспериментального исследования 


4.2.1. Цели экспериментального исследования 


Целями экспериментального исследования являются: 

• экспериментальное определение наилучших параметров алго- 
ритма моделирования системы искусственного интеллекта, выпол- 
няющей совмещение данных ДЗЗ и ЦКМ путем поиска образов 
объектов, присутствующих на ЦКМ, на спутниковом снимке — 
данная цель достигнута в рамках экспериментов № І, 2; 

• экспериментальное определение наилучших параметров алго- 
ритма моделирования системы искусственного интеллекта, выпол- 
няющей совмещение данных ДЗЗ и ЦКМ путем поиска корреляции 
между контурами образов, выделенных на спутниковом снимке, и 
границами объектов, присутствующих на ЦКМ - данная цель дос- 
тигнута в рамках эксперимента № 3; 

сравнение временных показателей и показателей качества сов- 
мещения, выполненного с помощью систем искусственного интел- 
лекта, и совмещения, выполненного рядом существующих распрост- 
раненных алгоритмов, данная цель достигнута в рамках эксперимен- 
та № 4; 

• оценка пригодности использования разработанных алгорит- 
мов моделирования, численных методов обучения систем искусст- 
венного интеллекта и непосредственно разработанных систем ис- 
кусственного интеллекта для выполнения высокоточного оператив- 
ного совмещения данных ДЗЗ и ЦКМ. 

Результаты совмещения оцениваются с временной точки 
зрения и с точки зрения точности выполненного совмещения. 


4.2.2. Условия проведения 
экспериментального исследования 


Инструменты проведения 
экспериментального исследования 


Экспериментальное исследование было выполнено с помощью 
программного стенда тестирования алгоритмов предварительной 
обработки и совмещения изображений. 
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Для обработки и анализа результатов экспериментов был 
использован язык программирования К. 

Экспериментальное исследование было проведено с использо- 
ванием вычислительного кластера, организованного по технологии 
Меѕѕаве раззше іпіегѓасе (МРІ), реализованной в виде программ- 
ного комплекса Ореп МРТ. Вычислительный кластер был составлен 
из следующих вычислительных систем: 

• управляющая вычислительная система кластера: 

— аппаратная часть: [пќе! Соге 1533370 (частота процессора 
2,7 ГГц; 4 вычислительных ядра; объем кэша 13-3 мегабай- 
та), АМР Кадеоп НР” 8750 М (частота процессора 775 МГц; 
объем ОЗУ - 2 гигабайта (технология ОЮКЗ); количество вы- 
числительных ядер - 192), 4 гигабайта ОЮКЗ ОЗУ; 

— программная часть: ОС СМОЛ іпих (дистрибутив АгсВ 
Глпих; версия ядра Глпих 3.14.4), программная библиотека 
ОТЬС версии 2.19; 

• вычислительная система № 1: 

— аппаратная часть: Іле! Соге! 7-4960Х (частота процессора 
3,6 ГГц; 6 вычислительных ядер; объем кэша [3-15 мегабайт), 
МУ!Ч!а СеРогсе СТХ 780 (частота процессора — 1071 МГц; 
объем ОЗУ - 3 гигабайта (технология Ср0ЮК5); количество 
вычислительных ядер - 192), 16 гигабайт ОЮКЗ ОЗУ; 

— программная часть: ОС СМОЛлпих (дистрибутив Аге 1пих; 

версия ядра Глпих 3.14.4), программная библиотека СБС 

версии 2.19; 

• вычислительная система № 2: 

– аппаратная часть: [пе] Соге 17-4770К (частота процессора 

3,5 ГГц; 4 вычислительных ядра; объем кэша [3-8 мегабайт), 

МУійа ОСеРогсе СТХ780 (частота процессора – 1750 МГц; 

объем ОЗУ - 3 гигабайта (технология СООК5); количество 
вычислительных ядер - 240), 16 гигабайт ОЮКЗ ОЗУ; 

— программная часть: ОС @МОЛлпих (дистрибутив Аге 1пих; 

версия ядра Піпих 3.14.4), программная библиотека СЫС 

версии 2.19. 

Каждое ядро процессоров вычислительных систем, включен- 
ных в кластер, использовалось в качестве вычислительного узла 
кластера и, таким образом, кластер состоял из 14-ти вычисли- 
тельных узлов (хостов). 
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Для уменьшения временных затрат на проведение экспери- 
ментального исследования каждый хост использовал мощности 
видеопроцессора, имевшегося в составе соответствующей системы, 
для дополнительного распараллеливания части вычислений с 
помощью технологии общих вычислений на видеокарте ((англ.) 
сепега|-ригрозе ртарћісѕ ргосеѕѕіпо ипіїѕ, СРСРО). 


4.2.3. Входные и выходные параметры 
экспериментального исследования 


Тестовые примеры 


Для проведения экспериментальных исследований было сфор- 


мировано множество Туи „и = {оа е  Т@стовых примеров 
[ 


=<В М Я и! > › ГДЕ: 


&ора/ 1е5ї 1е51 ? 1е51 ? 


° В, – спутниковый снимок, В, є Ве ; 
М... - ЦКМ; 
• Р, – множество пар реперных точек, одна из которых 


отмечена на В, ‚а другая – на М 


Ут. — Вектор, описывающий смещение В, „>, „‚ = 


ея ? тез! ? 


= оў 5е(В „У у) изображения В, 


еѕг ой 5е! ° 
На вход тестируемого метода совмещения поступает смещен- 
ный спутниковый снимок В и ЦКМ М. Задача исследуе- 


1е51 0} 5е! 
мого метода состоит в совмещении входных данных друг с другом. 
Результат совмещения оценивается с помощью среднеквадратич- 
ной меры ошибки Е (1.8) по множеству реперных точек Р 


1651 ° 


іеѕг ° 


В качестве источника данных ДЗЗ был использован сенсор 
«Тветайс таррег» (ТМ), установленный на КА «І апаѕаї 5». Для 
экспериментального исследования были использованы спектраль- 
ные каналы В 10 (синий), В 20 (зеленый), В 30 (красный), В 40 
(ближний инфракрасный) и И 50 (ближний инфракрасный) разре- 
шением 30 м на пиксель таким образом, К = 5. Размер спутни- 
кового снимка, полученного от КА «Г.ап45а{ 5», составляет 8231 на 
7511 пикселей. Для уменьшения временных затрат на проведение 
экспериментального исследования в качестве В, были использо- 


1051 
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ваны окна размером / = / =1024 пикселей, выделенные на исход- 
ных спутниковых снимках. 

Экспериментальным полигоном были выбраны Рязанская об- 
ласть, а также области и республики, граничащие с Рязанской 


областью. Множество тестовых данных ДЗЗ Вья состоит из спут- 
никовых снимков перечисленных регионов России за период с 1991 
по 2008 годы. 

В качестве источника ЦКМ был использован проект Ореп 
Ѕїгееѓ Марѕ [63]. 


Вектор у 


ой зе = {У етомай > Уконие > Узсае} ОПисывает геометрическое ис- 


кажение спутникового снимка В 


1е5ї ? 


которое необходимо скорректи- 


ровать посредством совмещения. Вектор Уот л 


состоит из следую- 
щих компонентов: 


• ү и! — СМещение в метрах (направление смещения не имеет 
значения, так как не существует зависимости информационной 
наполненности анализируемой сцены от направления ее обхода); 

6 у ик — УГОЛ поворота вокруг центра спутникового снимка 
относительно оси ОХ в радианах; 

• у – Масштабный коэффициент линейного преобразования 
координат, безразмерная величина. 


Множество Т, было составлено, исходя из следующих 


оба! 1е51 
соображений. 


Множество Т 


#їођаЇ ізі должно содержать репрезентативное 


представление всех возможных состояний окружающей среды, 
поскольку вариации спектральных яркостей и текстурных 
характеристик образов одной и той же сцены на спутниковых 
снимках зависят в наибольшей степени именно от этого. Большой 
размер /х/ спутниковых снимков, получаемых от КА «[Гапдѕаї 5» 
(полоса захвата 185 км), и выполнение дистрибьютором несколь- 
ких этапов предобработки обуславливает нормальное распределе- 
ние погодных условий на снимке с модой в точке: полуденная 
яркость солнечного света, отсутствие осадков, низкая облачность. 
Таким образом, главным параметром, определяющим специфику 
состояния окружающей среды, оказывается время года и зависящие 
от него состояние растительности и характер выпадающих осадков. 
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Для обеспечения репрезентативности представления множеством 
Това ем Генеральной совокупности всех возможных состояний 


окружающей среды, которые могут быть различимы на спут- 
никовых снимках в мере, значимой для оценки качества совмеще- 
ния, необходимо обеспечить равномерное распределение в 
множестве Т, тестовых примеров / относительно меся- 


1офа! 1еѕі &1офа! 1е51 


ца производства спутниковых снимков В 


1е$1 * 


Множество Т, должно являться репрезентативным предс- 


оба! 1е5ї 


тавлением генеральной совокупности векторов у,, „„, Которые мо- 


гут быть теоретически скорректированы предложенными систе- 
мами искусственного интеллекта. 


Мощность генеральной совокупности векторов у,,,,, 
определяется как (4.3): 
у ==:0 в 
гетоуа/ тах гетоуа/ тіп 
еее дак = 5 ! (4.1) 
У У 
у г 
ж ғгогаѓе тах ге тт ( 4 2) 
У, ае 
ж У се тах — Е тіп ( 4 3) 
К А 
ОУ се 
ГДЕ: У „тах Уо тт — ДИапазон возможных значений у, 


ду 


можных значений у 


— Шаг квантования у у у — диапазон воз- 


гетотаі! ? ’гоее тах ? гоіаіе тіп 


ду 


гоіаіе 


гетоуа! 
Шаг квантования у 
ду 


гоїіа1е ? ғо1іате р. 


у 


ѕсаіе тах ° У тіп 


— диапазон возможных значений у 


ѕсаіе ? 


— шаг 


ѕсаіе 


квантования у 


ѕсаіе * 
На основании инвариантности процесса совмещения к направ- 
лениям смещения и поворота можно принять: 


О тіп — 0; 
У оне пп — 0; 
Е тах — Л. 


На основании пространственного разрешения тестовых данных 
ДЗЗ можно принять: 
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бу 
можно рассчитать, исходя из угла поворота, 


= 30. 


гетоуа! `` 


Значение ду 


ғогаіе 
обеспечивающего (минимальное) смещение на один пиксель по 
диагонали по наиболее удаленной от центра спутникового снимка 
стороне: 


 2%302/2 
У ше = аГсЗШ-——————— = 0,06, 
(1024*30)/2 
из чего следует, что: 
У етоуа! тах — Уғетоуа! тіп нБ 52 

У ае 0, 06 
Очевидно также, что у, „ъъ = ОУ 
Таким образом, (4.3) сводится к: 

_— Унетома тах. 52. У саіе тах ый (4.4) 


о} ѕег тах 3 0 бу 


ѕсаіе 

Исходя из теоретических ограничений на качество работы 
алгоритма моделирования системы искусственного интеллекта, вы- 
полняющей совмещение данных ДЗЗ и ЦКМ путем поиска образов 
объектов, связанных с минимально допустимым количеством обра- 


зов объектов, пригодных для совмещения, 7, „„ = 20, можно рас- 


считать У, итак › ДЛЯ чего необходимо оценить по В, мак- 


симально возможное смещение, при котором количество объектов 
на соответствующей ЦКМ М, , совпадающих более чем на 50 % 


площади с соответствующими образами на В, , больше либо рав- 


е51 ? 


Но Т, пи · В ходе подготовки к проведению экспериментальных ис- 
= 1000. 


были оценены, исходя из типовых 


следований значение у 


гетоуа! тах 


было оценено как у 


гетоуа! тах 


Параметры у и ду 


ѕсаіе тах зсаіе 


значений коэффициента масштабного искажения, полученных в 
процессе обработки данных ДЗЗ, как у =3 и ду, =0,1. 


ѕсаіе тах ѕсаіе 


Окончательно № + гта» 
М аи ЕКЕ * 52 = 50267. (4.5) 
оў ѕег тах 3 0,1 


было оценено как: 
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Поскольку №,» ‚а тах >> Вен, ТО: 
| | Т ова! 1е51 |= Му 5ег тах — 50267, (4.6) 


при этом должно выполняться требование равномерного распреде- 


ления в Т 


гіоЊа! тезі месяцев производства спутниковых снимков В 


іеѕг ° 


Общие входные параметры экспериментов № 1 и № 2 


Общим параметром экспериментальных исследований алго- 
ритмов моделирования системы искусственного интеллекта, вы- 
полняющей совмещение данных ДЗЗ и ЦКМ путем поиска образов 
объектов, присутствующих на ЦКМ на спутниковом снимке, явля- 
ется множество Н =Н ОН содержащее в себе заранее 


классифицированные образы, выделенные на спутниковых 
снимках В 


іеасћ 1еѕї ігасћ 1е51 ? 


1е51 * 


Подмножества множества Н, и подмножества множества 
Н 


качества обучения простых классификаторов, входящих в состав 
древовидных классификаторов. 

Для корректного обучения классификаторов требуется, чтобы 
Н „һ и Н,„ были репрезентативными подмножествам генераль- 


ной совокупности образов объектов, пригодных для обучения клас- 
сификаторов в процессе обучения (настройки) разработанных сис- 
тем искусственного интеллекта. 

Древовидные классификаторы, обучаемые в ходе проведения 
экспериментальных исследований, должны классифицировать об- 
разы, относящиеся к образам водных объектов, лесных массивов и 
сельскохозяйственных полей, с вероятностью не менее 0,85, что, 
как показано в [17—19], является в среднем максимальной точнос- 
тью, которая может быть достигнута машиной опорных векторов и 
ИНС прямого распространения без обратных связей в процессе 
бинарной классификации. 

Для достижения заданной точности для каждого класса была 
рассчитана гистограмма многомерного распределения текстурных 
признаков Харалика образов объектов данного класса, выделенных 


применяются соответственно для обучения и тестирования 


1251 


1е51 


на спутниковых снимках Ви. Текстурные признаки Харалика бы- 


118 Экспериментальные исследования возможности использования 
систем искусственного интеллекта 


ли использованы по причине их более высокой информативности 
по сравнению с другими средствами описания текстур. 
Гистограмма распределения текстурных признаков Харалика 
для каждого из анализируемых классов сеС является унимо- 
дальной, следовательно, с помощью правила трех сигм возможно 
формирование множества Н ЕН содержащего обра- 


еасй 1е5ї с геасй 1е51 ? 
зы, относимые к классу с. 

При обучении каждого простого классификатора множества 
образов Н, и Н,, должны содержать одинаковое количество 
образов каждого из двух классов, причем подмножества, содержа- 
щие образы соответствующих классов, должны повторять пара- 
метры распределения образов данных классов в соответствующих 


множествах Н 


1е51 


іеасћ 1еѕі с‘ 


300, 
Н | = 700. Для формирования данных множеств на основе мно- 


Исходя из рекомендаций, изложенных в [17, 18], |Н 


іеасћ | Е 


жества Н в ходе проведения экспериментальных исследова- 


еасй еѕї 


ний был использован датчик псевдослучайных чисел, распреде- 
ленных нормально. 


Общие выходные параметры экспериментов 


Общими выходными параметрами экспериментов являются: 

• время выполнения совмещения (без учета времени обучения 
модели) с точностью до 1-й миллисекунды; 

• качество выполненного совмещения Е , оцененное по (1.8) и 
усредненное по 7., с точностью до 1-го метра. 


ораі 1е5! ? 


Значимость качественной эффективности совмещения превос- 
ходит значимость временной эффективности, так как теоретичес- 
кие предпосылки позволяют говорить о том, что дисперсия времен- 
ной эффективности работы каждой из исследуемых систем искусс- 
твенного интеллекта будет существенно меньше, чем дисперсия 
качественной эффективности. 


Глава 4 119 


Входные параметры, 
специфические для каждого эксперимента 


Различные комбинации фактических значений входных пара- 
метров, специфических для каждого эксперимента, образуют мно- 
жество векторов параметров Г „т = {У миа } - 


рағат 
Для каждого из рассматриваемых параметров максимально 
допустимое значение выбирается заведомо большим, чем 
соответствующие теоретические оценки. 
Эксперимент № 1. Входными параметрами-компонентами 
вектора /’„ – экспериментального определения наилучших с ка- 


ағат 
чественной и временной точек зрения параметров алгоритма моде- 
лирования системы искусственного интеллекта, выполняющей сов- 
мещение данных ДЗЗ и ЦКМ путем поиска образов объектов, при- 
сутствующих на ЦКМ на спутниковом снимке, с использованием 
древовидного классификатора и ИНС прямого распространения без 
обратных связей в качестве простого классификатора, являются: 
бит, — Размер нейронной карты Кохонена по каждой из 


размерностей. 

В ходе проведения эксперимента № 1 для простоты вычис- 
лений количество размерностей нейронной карты Кохонена будет 
принято равным 3, причем по каждой из размерностей нейронная 
карта будет иметь одинаковый размер (таким образом, в прост- 
ранстве координат нейронов используемая нейронная карта Кохо- 
нена может быть представлена как куб со стороной, равной 5, 
нейронов). 

Диапазоном возможных значений параметра был выбран диа- 
пазон 5, „„„„ © [1;51) ‚ шаг квантования параметра был принят 
равным 05 =5; 

• минимальный объем информации Т,„, , сохраняемый в клас- 
сифицируемом векторном пространстве в результате преобразова- 
ния Карунена – Лоева. 


Диапазоном возможных значений параметра был выбран диа- 
пазон Т,„ Е (0; ], шаг квантования параметра был принят равным 


ӘТ, = 0,1 > 
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• количество шагов развертки грамматики графовой генерации 
Китано Ч. 

Диапазоном возможных значений параметра был выбран 
диапазон $ Е1;10], шаг квантования параметра был принят равным 
ОФ =1; 

Максимальный уровень ошибки Т 


апп Е тах 


на тестовом мно- 


жестве векторов Н’.., при котором генетический алгоритм выбора 


1е51 ? 
структуры ИНС прямого распространения без обратных связей 
останавливается. 

Диапазоном возможных значений параметра был выбран диа- 


пазон Т; © (0;0,3] , шаг квантования параметра был принят 
равным ӧТ,,,„„ = 0,03. 
Количество векторов параметров равно: 
мы + у 03 =10000. (4) 
й 5 0,1 ] ‚ 03 


Эксперимент №2. Входными параметрами-компонентами 
вектора у экспериментального определения наилучших с 


рагат 
качественной и временной точек зрения параметров алгоритма 
моделирования системы искусственного интеллекта, выполняющей 
совмещение данных ДЗЗ и ЦКМ путем поиска образов объектов, 
присутствующих на ЦКМ, на спутниковом снимке, с использова- 
нием древовидного классификатора и машины опорных векторов в 
качестве простого классификатора, являются: 

• размер 5,„„„ нейронной карты Кохонена по каждой из раз- 
мерностей. 

В ходе проведения эксперимента № 1 для простоты вычисле- 
ний количество размерностей нейронной карты Кохонена будет 
принято равным 3, при чем по каждой из размерностей нейронная 
карта будет иметь одинаковый размер (таким образом, в прост- 
ранстве координат нейронов используемая нейронная карта Кохо- 
нена может быть представлена как куб со стороной, равной 5, 
нейронов). 
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Диапазоном возможных значений параметра был выбран 
диапазон 5,,„„„ Е [1;51), шаг квантования параметра был принят 


равным 05. = 5; 
• минимальный объем информации Т, , сохраняемый в клас- 


сифицируемом векторном пространстве в результате преобразова- 
ния Карунена – Лоева. 

Диапазоном возможных значений параметра был выбран 
диапазон Т,, Є (0;1] , шаг квантования параметра был принят 


равным ӧТ,,, = 0,1; 


• параметр $ обучения машин опорных векторов, опреде- 


5ит 


ляющий верхнюю границу поиска множителей Лагранжа в ходе 
оптимизации соответствующей функции. 

Диапазоном возможных значений параметра был выбран диа- 
пазон 5 є (0;10], шаг квантования параметра был принят равным 


5ит 
С =0,1. 
Количество векторов параметров равно: 


51-1 1 10 
=——.—_.—_ = 10000. 4.8 
| рағат | 5 0,1 0,1 ( 


Эксперимент № 3. Входными параметрами-компонентами 
вектора у экспериментального определения наилучших с 


рағат 
качественной и временной точек зрения параметров алгоритма 
моделирования системы искусственного интеллекта, выполняющей 
совмещение данных ДЗЗ и ЦКМ путем поиска корреляции между 
контурами образов, выделенных на спутниковом снимке, и грани- 
цами объектов, присутствующих на ЦКМ, являются: 

• размер 5 окна выделения контуров (5 „= М' = № 


сопіои" $12е сотоиг 5іс 


для нелинейной ячеистой нейронной сети и 5 „ге = 5 зао ли" ДЛЯ 


нейросетевого алгоритма выделения контуров, основанного на 
использовании вейвлета Габора). 

Диапазоном возможных значений параметра был выбран диа- 
пазон 5 є [3;30], шаг квантования параметра был принят 


ее 


сотоиг $12е 


равным 05 


сотоиг 5126 
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• количество $ н уі2іЫе 
стохастической машины Больцмана. 
Диапазоном возможных значений параметра был выбран диа- 


пазон 5 є [1;30], шаг квантования параметра был принят 


• количество 5, „м, Нейронов скрытого слоя ограниченной 


нейронов видимого слоя ограниченной 


гьт уе 


равным 05 


гфт ме 


стохастической машины Больцмана. 
Диапазоном возможных значений параметра был выбран диа- 
пазон 5„„ „и, Е [1:30] , шаг квантования параметра был принят 


равным 5, =1. 


Количество векторов параметров равно: 


Ра (2022: у о) [92-900 (4.9) 


Эксперимент № 4. Экспериментальное сравнение качествен- 
ных и временных показателей совмещения, выполненного с по- 
мощью разработанных систем искусственного интеллекта, и совме- 
щения, выполненного рядом существующих распространенных 
алгоритмов, не имеет специфических входных параметров. 


4.2.4. План проведения экспериментальных исследований 


План проведения экспериментов № 1 и № 2 


Эксперименты № 1 и №2, предполагающие эксперименталь- 
ную проверку качественной и временной эффективности алгорит- 
мов моделирования системы искусственного интеллекта, выпол- 
няющей совмещение данных ДЗЗ и ЦКМ путем поиска образов 
объектов, присутствующих на ЦКМ на спутниковом снимке, были 
выполнены независимо по одинаковому плану, реализующему 
полный факторный эксперимент по компонентам векторов у 


рагат ' 


План проведения экспериментов № | и № 2 состоял из следую- 
щих этапов: 


ВИ =У 


рагат сиғепі рагат ? 
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2) обучение древовидного классификатора для каждого набора 
параметров из Г с использованием множества образов 


рағат сиггет 


еасй 1еѕт 5 
3) обработка каждым алгоритмом моделирования 1 % тестовых 


примеров из Т... е › выбранных с помощью датчика псевдо- 


случайных чисел, распределенных равномерно; 

4) построение гистограммы распределения качества совмеще- 
НИЯ; 

5) поиск мод гистограммы; 

6) формирование множества И 


ағат һе ВеКТОров параметров экс- 
перимента, лежащих на модах гистограммы, такого, что дисперсия Е 


по И не должна превышать 0,1 (т.е. отбор лучших векторов 


рагат Беѕї 


параметров у алгоритмов моделирования с точки зрения ка- 


рағат 
чества совмещения для случайног о подмножества тестовых при- 


меров Т. „ьа „„› На основании предположения о равномерном расп- 


ределении оценки качества совмещения по ова 389 ); 
7) У, =У 


агат сиғғепі рагат Без! ? 


8) повторение пунктов 2-7 для Г,, последовательно с 


гат сиггепі 
использованием 10%, 25%, 50%, 75 % тестовых примеров из 


Т оова еч На этапе 3, что обеспечит постепенное повышение оценки 


точности качества совмещения для лучших векторов параметров из 


рагат * 
Для оценки временной эффективности совмещения для каж- 
дого вектора параметров у временные затраты на непосредст- 


рагат 
венное совмещение (без учета этапа обучения классификатора) 
аккумулировались и по окончанию эксперимента нормировались 
по количеству выполненных совмещений. 


План проведения эксперимента № З 


Эксперимент № 3, предполагающий экспериментальную про- 
верку качественной и временной эффективности алгоритма моде- 
лирования системы искусственного интеллекта, выполняющей сов- 
мещение данных ДЗЗ и ЦКМ путем поиска корреляции между 


124 Экспериментальные исследования возможности использования 
систем искусственного интеллекта 


контурами образов, выделенных на спутниковом снимке, и грани- 
цами объектов, присутствующих на ЦКМ, был выполнен по плану, 
реализующему полный факторный эксперимент по компонентам 


вектора ут. 


План проведения эксперимента № 3 состоял из следующих 
этапов: 
1) Р, = / 


ағат сиггепі а 
2) обработка каждым алгоритмом моделирования 1 % тестовых 
примеров из Т, выбранных с помощью датчика псевдослу- 


ођа! іеѕі ? 
чайных чисел, распределенных равномерно; 

3) построение гистограммы распределения качества совме- 
щения; 

4) поиск мод гистограммы; 

5) формирование множества Г 


азан. Вее векторов параметров 
эксперимента, лежащих на модах гистограммы, такого, что дис- 


персия Е по И не должна превышать 0,1 (т.е. отбор лучших 


рагат вез 


векторов параметров у алгоритмов моделирования с точки зре- 


рагат 
ния качества совмещения для случайного подмножества тестовых 


примеров Т„„„,„„› Исходя из предположения о равномерном 


распределении оценки качества совмещения по Той РҮ ); 
6) И = 


рагат сиғғепі рағат Беѕї ? 


7) повторение пунктов 2—6 для Г 


рағат сиггет 


последовательно с 


использованием 10 %, 25 %, 50 %, 75 % тестовых примеров из 
Т оова вы На этапе 3, что обеспечит постепенное повышение оценки 


точности качества совмещения для лучших векторов параметров из 
рағат ° 


Для оценки временной эффективности совмещения для каж- 
дого вектора параметров у временные затраты на непосредст- 


рагат 


венное совмещение аккумулировались и по окончанию экспери- 
мента нормировались по количеству выполненных совмещений. 
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План проведения эксперимента № 4 


Эксперимент № 4 предполагал сравнение результатов экспери- 
ментов № 1 и № 3 с оценками времени и качества совмещения, вы- 
полненного следующими методами и алгоритмами: 

• корреляционно-экстремальный метод совмещения, рассмот- 
ренный в [48]; 

• метод стабилизации изображений, основанный на преобразо- 
вании Фурье; 

• метод, основанный на поиске на спутниковом снимке и ЦКМ 
особых точек с помощью детектора $ІЕТ; 

• метод, основанный на поиске на спутниковом снимке и ЦКМ 
особых точек с помощью детектора 50КЕ; 

• метод, основанный на поиске на спутниковом снимке и ЦКМ 
особых точек с помощью детектора Харриса. 

План проведения эксперимента № 4 состоял в обработке мно- 


жества тестовых примеров Т, „„ Каждым из перечисленных 


алгоритмов, усреднении временных показателей и показателей 
эффективности с точки зрения качества совмещения и сравнении 
полученных результатов с результатами экспериментов № 1 и №3. 


4.3. Результаты экспериментального исследования 


4.3.1. Методика анализа результатов 
экспериментального исследования 


Методика анализа 
результатов экспериментов № 1 и № З 


Целями анализа результатов экспериментальных исследований 
качественной и временной эффективности алгоритмов моделирова- 
ния систем искусственного интеллекта, выполняющей совмещение 
данных ДЗЗ и ЦКМ, стали: 

• поиск наилучших с качественной и временной точек зрения 
наборов параметров исследуемых алгоритмов моделирования сис- 
тем искусственного интеллекта, выполняющих совмещение данных 
ДЗЗ и ЦКМ; 

• визуализация зависимостей выходных параметров моделиро- 
вания (временной оценки и оценки качества выполненного сов- 
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мещения) на наборе тестовых примеров Тин 


от параметров 
исследуемых алгоритмов моделирования. 

Алгоритм визуализации означенных зависимостей состоит из 
следующих этапов: 

1) расчет среднеквадратичного отклонения исследуемого вы- 
ходного параметра эксперимента для значения входного параметра 
для каждого из входных параметров, специфических для анали- 
зируемого эксперимента; 

2) выбор двух входных параметров эксперимента с наиболь- 
шим СКО; 

3) построение трехмерного графика зависимости исследуемого 
выходного параметра эксперимента от двух входных параметров с 
наибольшим СКО. 


Методика анализа результатов эксперимента № 4 


Целью анализа результатов экспериментального сравнения вре- 
менных показателей и оценок качества совмещения, выполненного 
с помощью разработанных систем искусственного интеллекта, и 
совмещения, выполненного рядом существующих распространен- 
ных алгоритмов, стала визуализация результатов сравнения в виде 
сводной таблицы, содержащей следующие данные: 

• для каждой из систем искусственного интеллекта: 

_ оценка качества совмещения и оценка временной эффек- 
тивности совмещения для наилучшего с точки зрения ка- 
чества совмещения набора параметров у данной систе- 


рағат 
мы искусственного интеллекта; 

— оценка качества совмещения и оценка временной эффек- 
тивности совмещения для наилучшего с временной точки 
зрения набора параметров у данной системы искусст- 


рағат 
венного интеллекта; 

• для каждого из существующих распространенных алгоритмов 
совмещения — оценка качества совмещения и оценка временной 
эффективности совмещения, усредненные по множеству тестовых 
примеров 7, 


ора! тез! ° 
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4.3.2. Результаты экспериментального исследования 
Эксперимент № 1 


В табл. 4.1 приведены значения среднеквадратичного отклонения 
(СКО) оценки качества совмещения Е (1.8) и СКО оценки времени 


совмещения для каждого компонента вектора у входных Пара- 


рағат 


метров эксперимента № 1. 
На рис. 4.1 приведена зависимость функции оценки качества 
совмещения Ё от параметров 7,„„ и Ч эксперимента № 1. 


Таблица4.1. СКО выходных параметров эксперимента 
ЕЕ СКО оценки качества СКО оценки времени 
Е совмещения ЕЁ совмещения 
0,313 0,516 


0,750 0,498 


0421 0791 


На рис. 4.2 приведена зависимость времени выполнения совме- 
эксперимента № 1. 


щения от параметров Фи Т 


апп Е тах 


0.6 


ф 0.4 


2 0.2 Ттуо 


Рис. 4.1. Экспериментальная зависимость функции оценки качества 
совмещения Е от параметров Г„„ и‘ эксперимента № | 
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В табл. 4.2 приведены векторы параметров у„„„ эксперимен- 


та, наилучшие с точки зрения временной и качественной эффектив- 
ности совмещения. 


Тель Е тов 5. 10-2 


Рис. 4.2. Экспериментальная зависимость времени выполнения 
совмещения от параметров Фи ТГ, .„, эксперимента № 1 


Таблица 4.2. Результаты эксперимента № 1 


Оценка качества | Оценка времени 
совмещения Ё совмещения 
а (метры/пиксели) (секунды) 


5 ковопен =] 1, Тир Е 0,300; Ҹ=8; 
23,546 / 0,785 
я Е тах 58 0, 120. 


Вектор параметров 
эксперимента у 


понт 0100; 
У ТУЯ іпјо 24,813 у 0,827 
8 Горы = 05120. 
$ доһопеп — 16; 125 = 0,200; 
Ф=10; я Е тах — 0,150. 26,204 / 0,873 
$ доропеп г 1; Тољ = 0,500; 
26,322 / 0,877 
Е ра 0126 
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Окончание табл. 4.2. 


29,896 / 0,997 


В табл. 4.3 приведены значения СКО оценки качества сов- 
мещения Ё (1.8) и СКО оценки времени совмещения для каждого 
компонента вектора у входных параметров эксперимента № 2. 


Т, 


в = 0,200; 


= 0,150. 

= 0,300; 

= 0,240. 

Тиь = 0,400; 

= 0,090. 

Г 8 = 0,300; 
= 0,090. 

Ти = 0,300; 
= 0,240 

Ти» = 0,400: 
= 0,300 


$ доропеп = 6; 


Ф = 7; Т 


апп Е тах 


Тир 


5 Ковонен =6; 


ф = 10: Т 


апп Е тах 


ойвиен = 6; 


= 1; Т 


апп Е тах 


ОО г ] 1; 


= 8; Т 


апп Е тах 


5 Койопеп = 2 І | 


Ҹ=6: Т 


апп Е тах 


5 ароми Б ] 1; 


= 5; Т 


апп Е тах 


Эксперимент № 2 


рагат 
На рис. 4.3 приведена зависимость функции оценки качества 
совмещения Е от параметров Т„„ и С эксперимента № 2. 


5ит 
На рис. 4.4 приведена зависимость времени выполнения совме- 
щения от параметров 5 „И С,„„ эксперимента № 2. 


5ит 


Таблица 4.3. СКО выходных параметров эксперимента 
Параметр СКО оценки качества СКО оценки времени 
совмещения Е совмещения 
0,328 0,601 


0,645 0,731 
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Сет 2 02 | Ттуо 


Рис. 4.3. Экспериментальная зависимость функции оценки качества 
совмещения Ё от параметров Т„„ и С эксперимента № 2 


Рис. 4.4. Экспериментальная зависимость времени выполнения 
совмещения Ё от параметров 5 и С,,„ эксперимента № 2 


5ит 


В табл. 4.4 приведены векторы параметров у экспери- 


рагат 


мента, наилучшие с точки зрения временной и качественной 
эффективности совмещения. 
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Таблица 4.4. Результаты эксперимента № 2 


Оценка качества 
совмещения Е 
(метры/пиксели) 


Оценка времени 
совмещения 
(секунды) 


1,139 
22,638 / 0,755 1,981 
23,118 / 0,771 1,913 
5 дойопеп =] 1; 
24,146 / 0,805 3,930 
Сеа е 9,900. 
5 ковопет | Г; Тир — 0,400; 
С, = 9,900 24,146 / 0,805 3,930 
Тир = 0,100; 
24,472 / 0,816 2,561 
Гь = 0,300; 34,512/1,150 0,755 
27,534/0,918 0,760 
28,781/0,959 0,817 
5 = 1; Г 
Койопеп іпјо 25,212/0,840 0,819 
С, = 6,800 
ии, 5 Ть= 0,200; 
Е: 25,452/0,848 0,849 
С, = 9,300 


Вектор параметров 
эксперимента у 


рағат 


Т 


о = 0,900; 


рода — | 1; 


С = 0,100 


эит 


22,088 / 0,736 


Т 


в = 1,000; 


ОБИ = ] 1; 


С. =1,800 


5ит 


Т 


в = 0,100; 


бай = 6; 


Со = 1,900 


5ит 


Т 


о = 0,400; 


СЕРЕ ра 
С = 4,400. 
5 онопеп = 46; 
С, = 9,700 


ѕит 
Гн јо 


5ит 


= 0,300; 


Вей а 1; 


С. „= 2,200 


5ит 


Т 


в = 0,300; 


5 Койопен =46; 


С = 0,100 


= 0,200; 
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Эксперимент № 3 


В табл. 4.5 приведены значения СКО оценки качества совмещения 
Е (1.8) и СКО оценки времени совмещения для каждого компонента 
вектора у входных параметров эксперимента № 3. 


рагат 
На рис. 4.5 приведена зависимость функции оценки качества 
совмещения Е от параметров 5, И 9 „и ЭКСПеримента № 3. 


фт уігіЫе 


Таблица 4.5. СКО выходных параметров эксперимента 
Параметр СКО оценки качества СКО оценки времени 
совмещения Ё совмещения 
0298 


061 0451 
0689 0,482 


5 


гьт уі2іЫе 


Рис. 4.5. Экспериментальная зависимость функции оценки качества 


совмещения Е от параметров $, И 5, ЭКсперимента № 3 


фт уігіЫе 


На рис. 4.6 приведена зависимость времени выполнения сов- 
мещения Е от параметров 5 И 5 „рае, ЭКсперимента № 3. 


сотоиг 5і2е 


экспери- 


агат 


В табл. 4.6 приведены векторы параметров у, 


мента, наилучшие с точки зрения временной и качественной 
эффективности совмещения. 
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ныть маеп 5. 5 беотбоит гќхе 


Рис. 4.6. Экспериментальная зависимость времени выполнения 
совмещения Ё от параметров 5 И Эми ЭКсперимента № 3 


сопіоиг 5і2е 


Таблица 4.6. Результаты эксперимента № 3 


Вектор параметров 
эксперимента у 


Оценка качества | Оценка времени 
совмещения Ё совмещения 
(метры/пиксели) (секунды) 


рағат 


К = 27; 
$ н тааеп = 29 


Ау =21; 


К) 


гЬт уЫе 


сотоиг ѕіге 


К) 


гЬт уігіЫе 


= 29; 


сотоиг 5і2е 


Е радеп 5 г, 


5 сотоиг ѕіге с 2 | , Эра уігіЫе Бе 
5 гьт ааенп = 2 8 


5 = 27; 5 


гфт уі2іЫе = 


сотоиг 512е 


н ааеп Е 9 


К) 


; 5 


ғБт зле 


=. 


сотоиг Я ге 


о рааеп = 


$ =30; 
5 


5 


гфт уігіЫе 


сотоиг $1те 


ғЬт шааеп г 19 
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Окончание табл. 4.6. 


В табл. 4.7 приведены результаты экспериментального сравне- 
ния временной эффективности и оценки качества совмещения, 
выполненного с помощью разработанных систем искусственного 
интеллекта, и совмещения, выполненного рядом существующих 
распространенных алгоритмов. 

В табл. 4.7 используются следующие обозначения: 

• ЕЕ-обозначает алгоритм моделирования системы искусствен- 
ного интеллекта, выполняющей совмещение данных ДЗЗ и ЦКМ 
путем поиска образов объектов, присутствующих на ЦКМ, на спут- 
никовом снимке с использованием древовидного классификатора и 
ИНС прямого распространения без обратных связей в качестве 
простого классификатора; 

• 5УМ-обозначает алгоритм моделирования системы искусст- 
венного интеллекта, выполняющей совмещение данных ДЗЗ и ЦКМ 
путем поиска образов объектов, присутствующих на ЦКМ, на спут- 
никовом снимке с использованием древовидного классификатора и 
машины опорных векторов в качестве простого классификатора; 

• КВМ-обозначает алгоритм моделирования системы искусст- 
венного интеллекта, выполняющей совмещение данных ДЗЗ и 
ЦКМ путем поиска корреляции между контурами образов, выде- 
ленных на спутниковом снимке, и границами объектов, присутст- 
вующих на ЦКМ. 

Процесс совмещения проиллюстрирован рис. 4.7—4.16. 


9 = 6; 5 


гфт ааен я 28 


= 27; 5 


сотоиг $12е 


7; 


сотоиг $12е гЬт уігіЫе 


6; 


Вт уігіре = 
ЭО 


= 30; 5 


гьт уігіЫе =? 


ЕЕ 


сотоиг 5і2е 


ғЬт шааеп = 30 


24; н уе 7. и. 


сотоиг 5і2е 


9 
9 
9 
9 
9 
9 
9 


гфт ћіадеп = 2 9 
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На рис. 4.7 приведен исходный спутниковый снимок Рязанской 
области, сделанный камерой ТМ, установленной на КА «І апаѕаї 5». 
Спутниковый снимок сделан 29 июля 2011 года. 

На рис. 4.8 приведен результат комплексирования исходного 
спутникового снимка и ЦКМ. Невязка между спутниковым сним- 
ком и ЦКМ обусловлена смещением по обеим координатным осям, 
а также поворотом спутникового снимка вокруг своего центра. 

На рис. 4.9 и 4.10 приведены результаты совмещения исход- 
ного спутникового снимка и ЦКМ, выполненного путем поиска 
образов объектов, присутствующих на ЦКМ, на спутниковом сним- 
ке с использованием древовидного классификатора. На рис. 4.9 при- 
ведены результаты совмещения, полученные с применением ИНС 
прямого распространения без обратных связей в качестве простого 
классификатора. На рис. 4.10 приведены результаты сов-мещения, 
полученные с применением машины опорных векторов в качестве 
простого классификатора. 

На рис. 4.11 приведены результаты совмещения исходного спут- 
никового снимка и ЦКМ, выполненного путем поиска корреляции 
между контурами образов, выделенных на спутниковом снимке, и 
границами объектов, присутствующих на ЦКМ. 

На рис. 4.12—4.16 приведены результаты совмещения исходно- 
го спутникового снимка и ЦКМ, выполненного с помощью кор- 
реляционно-экстремального метода совмещения (рис. 4.12), с по- 
мощью метода стабилизации изображений, основанного на преоб- 
разовании Фурье (рис. 4.13), с помощью метода, основанного на 
поиске на спутниковом снимке и ЦКМ особых точек с помощью 
детектора ЗТЕТ (рис. 4.14), с помощью метода, основанного на 
поиске на спутниковом снимке и ЦКМ особых точек с помощью 
детектора ОКЕ (рис. 4.15), с помощью метода, основанного на 
поиске на спутниковом снимке и ЦКМ особых точек с помощью 
детектора Харриса (рис. 4.16). 

Оценка качества выполненного совмещения и оценка време- 
ни, затраченного на совмещение, приведены для данного примера 
в табл. 4.8. В табл. 4.8 используются те же обозначения, что и в 
табл. 4.7. 
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Таблица 4.7. Результаты эксперимента № 4 
Метод/алгоритм совмещения Оценка качества | Оценка времени 
совмещения Ё совмещения 
(метры/пиксели (секунды) 
Е 1,546/0,052 1,212 
1,341/0,045 0,755 
1,088/0,036 2,139 


Е 
1,821/0,061 0,781 
метод 
изображений, основанный на 121,361/4,045 1,981 
реобразовании Фурье 
Метод, основанный на поиске 
на спутниковом снимке и ЦКМ 91,421/3,047 


особых точек с помощью 
75,941/2,531 14,261 


детектора МЕТ 
Метод, основанный на поиске 

Таблица 4.8. Оценка качества выполненного совмещения 
и оценка времени, затраченного на совмещение 


на спутниковом снимке и ЦКМ 
особых точек с помощью 
детектора ЗОВЕ 
Метод, основанный на поиске 
на спутниковом снимке и ЦКМ 
особых точек с помощью 
детектора Харриса 


Оценка времени 
совмещения 


(секунды) 


9187. 


Оценка качества 
совмещения ЕЁ 
(метры/пиксели) 


Метод/алгоритм совмещения 


Е 
М 27,918/0,931 1,241 


Корреляционно-экстремальный 45,321/1,511 1,971 
метод 
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Окончание табл. 4.7. 


98,89 1/3,296 1,760 
82,473/2,749 13,091 
67,537/2,251 
54,83 1/1,828 


Метод стабилизации 
изображений, основанный 
на преобразовании Фурье 
Метод, основанный на поиске 
на спутниковом снимке и ЦКМ 
особых точек с помощью 
детектора МЕТ 
Метод, основанный на поиске 
на спутниковом снимке 
и ЦКМ особых точек с 
помощью детектора 50КЕ 
Метод, основанный на поиске 
на спутниковом снимке 
и ЦКМ особых точек с 
помощью детектора Харриса 


Рис. 4.7. Исходный спутниковый снимок (зеленый, красный и 
ближний инфракрасный спектральные каналы; Рязанская 
область; камера ТМ КА «ІГ апаѕаї 5») 


138 Экспериментальные исследования возможности использования 
систем искусственного интеллекта 


| _ Окно 1 | 


Рис. 4.8. Результат комплексирования 
исходного спутникового снимка и ЦКМ 


Окно | Окно 2 Окно 3 


Рис. 4.9. Результаты совмещения спутникового снимка и ЦКМ с 
использованием древовидного классификатора и ИНС прямого 
распространения без обратных связей 
в качестве простого классификатора 
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Рис. 4.10. Результаты совмещения спутникового снимка и ЦКМ 
с использованием древовидного классификатора и машины опорных 
векторов в качестве простого классификатора 


Окно | Окно 2 Окно 3 


Рис. 4.11. Результаты совмещения спутникового снимка и ЦКМ путем 
поиска корреляции между контурами образов, выделенные на снимке 
и границами объектов, присутствующих на ЦКМ 


аве 


т ФЕ 


А, 


Окно | Окно 2 Окно 3 


Рие. 4.12. Результаты совмешения спутникового снимка и ЦКМ, 
выполненного с помощью корреляционно-экстремального метода 
совмещения 
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Рис. 4.13. Результаты совмещения спутникового снимка и ЦКМ, 
выполненного с помощью метода стабилизации изображений, 
основанного на преобразовании Фурье 


|4 


Окно І Окно 2 Окно 3 


Рис. 4.14. Результаты совмещения спутникового снимка и ЦКМ, 
выполненного с помощью метода, 
основанного на поиске на спутниковом снимке 
и ЦКМ особых точек с помощью детектора МЕТ 


Г 4 


4 


и |) 
Окно | Окно 2 Окно 3 


Рис. 4.15. Результаты совмещения спутникового снимка и ЦКМ, 
выполненного с помощью метода, 
основанного на поиске на спутниковом снимке 
и ЦКМ особых точек с помощью детектора ЗОВЕ 
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Рис. 4.16. Результаты совмещения спутникового снимка и ЦКМ, 
выполненного с помощью метода, 
основанного на поиске на спутниковом снимке 
и ЦКМ особых точек с помощью детектора Харриса 


Заключение 


Описанные алгоритмы моделирования систем искусственного 
интеллекта, численные методы обучения данных систем и комп- 
лекс программ, реализующий перечисленные алгоритмы моделиро- 
вания и численные методы, позволяют выполнять высокоточное 
оперативное автоматическое совмещение данных ДЗЗ и ЦКМ, что 
подтверждено соответствующими экспериментальными исследова- 
НИЯМИ. 
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